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BACKGROUND 


NASA has recognized that continually advancing computational capability in all its many 
facets is an absolute necessity to carry out its present and future roles and missions. As 
a premier research and development agency of the U.S. Government NASA-MSFC is 
forecasting its needs in the larger context of national policies. It is recognized that 
MSFC’s current and upcoming needs deal in large measure (but not exclusively) with 
the control and utilization of those elements of the nation's space effort that deal with 
large propulsion systems, launch vehicles and orbital or payload structures (including 
the Space Station and it appendages). Although policy at the National level is still 
evolving, it is safe to assume that constrained resources are all that are foreseeable in 
the near future. Therefore, careful planning to allot available computational resources is 
mandatory. This effort addresses itself to questions pertaining to this allotment. 
Material gathered during the course of this investigation is presented as are specific 
recommendations of actions that MSFC should consider taking. 
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APPROACH 


To form the basis of recommendations as to how MSFC should allot (and justify) its 
computational resources, a number of different activities were performed. These 
included in house (at UAH) computing projects; a trip to AMES Research Center to 
attend a workstation workshop; various demonstrations held at the UAH Visualization 
and Simulation laboratory; discussion with principals in the field and literature searches. 
The details of many of these activities are contained in the appendices. 

Specifically, these activities include the following. Discussions with and a presentation 
by the SDI sponsor and Georgia Institute of Technology (GIT) developer of the Special 
Purpose Operational Computing Kernel (SPOCK). This device is a unique parallel 
functional processor being developed specifically for SDI types of intercepter missions. 
Detail of SPOCK are contained in Appendix A. A discussion was held with Dr. Mike 
Sayers at the University of Michigan Transportation Institute. Dr. Sayers has developed 
AUTOSIM, a code generator for systems of rigid bodies AUTOSIM, is written in 
common LISP and runs on a MAC computer. It does not have provision for flexible 
mode incorporation, is copyrighted by the University of Michigan and costs 9000 dollars. 
Information about it is contained in Appendix B. A fascinating workshop on workstations 
technology was attended which was put on at AMES mainly by USRA. This workshop 
covered numerous applications of computers, but in the main to subjects far afield from 
MSFC's needs and indeed was very futuristic in many of its aims, goals, and concepts. 
A full accounting of the workshop is enclosed as Appendix C. At MSFC's suggestion, 
arrangements were made to borrow a limited time demonstration copy of SD/FAST 
SD/EXACT. Before receiving this program, the author spent a week at ADI gaining 
familiarity with the ADSIM language. This is a program which generates the equations 
of motion (in FORTRAN) of multiple connected rigid bodies. In addition, SD/FAST has 
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the necessary software "hooks" to let it be run on the Applied Dynamics AD- 100 which 
is a parallel processor digital machine optimized to the computational tasks found in 
simulating continuous systems and capable of 20-MFLOPS. UAH has one of these 
machines and MSFC has at least four, two of which are fiber optic linked and dedicated 
to real time SSME simulation. This code was exercised on the VAX front end of the AD- 
100 to generate equations of motion which were then executed along with control logic 
programmed in the ADSIM language (Version VII) on the AD-100. Dr. Li USRA/ED-12 
also ran a problem of his own. His results were transmitted to MSFC on a 5 1/4 " floppy 
disk for his use there. This was a very interesting exercise and is covered at some 
length in Appendix D. Once again at MSFC's suggestion, two demonstrations of 
CADSIs DADS (Dynamic Analysis and Design System) software were held at the UAH 
Visualization and Simulation (V&S) Laboratory. This software system includes 
capabilities for static, kinematic, dynamic analysis, inverse dynamics, inclusion of 
vibratory model information, advanced graphics and a library of feedback, control and 
hydraulic components. Both demonstrations were attended by MSFC personnel, as well 
as by some personnel employed by local industry e.g., TBE. Information on DADS is 
covered in Appendix E. MSFC's TREETOPS program has just been transferred into the 
V&S lab, but time ran out on exercising it in the AD-100 environment. In addition, of 
course, considerable general literature was read covering the field. 
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RECOMMENDATIONS 


Based upon the information emassed during the activities outlined above, the following 
near term recommendations for MSFC activity emerged. 

One; even though it did not appear that the AMES/USRA, computer activities are very 
applicable to MSFC programs at this time, they should be monitored on a continuing 
basis. It seemed clear that JPL, AMES, Stanford and others have a "good thing going" 
and have their own ideas, largely science oriented, as to how things should go. They 
also clearly have the ear of NASA Headquarters. 

Two; find the means to join the Center for Simulation and Design Optimization of 
Mechanical systems (see Appendix F for some detailed information). GSFC and LARC 
are members along with many other Government and Industrial organizations. This 
organization is advancing all the aspects of Mechanical Simulation that appear quite 
applicable to MSFC programs and projects. 

Three; start an activity at MSFC using several optically linked AD-100 computers 
currently available at MSFC modeling an LSS of some complexity. This to be done in 
real time with provision for using the AD-100's Real-Time/Input Output devices 
connected to a control computer. This would develop some hands on experience with 
real-time parallel computing (with HWIL too). Not only would this have minimum initial 
cost impact, but also if a verified model of say the LSS beams hanging in building 4619 
were programmed many more control strategies per unit time could be explored. 
Because of its possession of an AD-100 the UAH V&S Lab could be a substantial 
contributor to this activity. 
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Four; if TREETOPS is to be maintained as a contemporary tool then it should have 
continued effort on its graphics e.g., add color, perhaps expand the number of 
computing engines especially workstations on which it can readily be installed and 
continue algorithmic work to enhance its speed of execution. A somewhat more 
ambitious avenue would be to re-code TREETOPS for running on massively parallel 
machines and provide the necessary I/O so that real-time simulation with HWIL can be 
readily accomplished. 

Five; initiate an effort to become integrated into the massively parallel computer world. 
A current example of this type of effort is the 1 5 million dollar contract that Intel has 
been awarded to build a machine with 32 billion floating point operations per second 
(theoretical) using 528 of Intel's i 860 microprocessors. It is supposed to be completed 
by next spring for the Concurrent Supercomputing Consortium. This is a group of 14 
U.S. research organizations (mostly Federal) led by the California Institute of 
Technology. At the present time, the so called Delta System is to be applied to military 
and "grand challenges" (e.g. global warming, chemical reactions). Why not try to make 
LSS control one of the "grand challenges"? 

Six; hold meetings of national interest on the subject of computing and its applications to 
matters of interest here in Huntsville. UAH's Beville Center is equipped to handle all the 
arrangements if it is so desired. 

Seven; attend several well chosen symposia and national meetings each year on this 
subject. Present paper(s) if possible. 

Eight; establish and maintain meaningful relationships with the computer related parts of 
regional Universities so that bilateral enrichment of each institution occurs at an 
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intellectual level and also to assure a supply of properly trained personnel into the 
workforce. 

These recommendations are presented in WBS form in Appendix G. 
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CONCLUSIONS 


A large computing thrust is a national priority Larae scalp rnmm +■ 
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Appendix A 

Special Purpose Operational Computing Kernel (SPOCK) 


Special thanks are due to Mr. Buster E. Kelley of the U.S. Army Strategic Defense 
Command. Through his good offices the Georgia Institute of Technology developers 
(Dr. Alford) came to MSFC and made a presentation on SPOCK. 

A new innovative concept, the SPOCK, is being developed by GIT. The SPOCK 
concept, which utilizes an array of function processors connected by a solid state 
crossbar matrix, provides three major benefits over conventional processing concepts. 
These are: 1) One-to-oneness with the block diagram of the system implemented, 2) 
emulation of flight hardware, and 3) real-time performance evaluation of flight hardware 
prior to flight hardware fabrication. 

Arrangements were made with GIT to travel to Atlanta and exercise the machine on a 
problem of interest to MSFC, however this never occurred. 
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Appendix B 


AUTOSIM Code Generator 
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Through Mr. J. Munson of ADI, we were made aware of the AUTOSIM code. This code 
was developed by Dr. Mike Sayers of the University of Michigan. It deals only with 
connected rigid bodies (no vibrational modes). In a conversation with Dr. Sayers, it 
became apparent that the users for whom it was developed were the automotive 
designers. I was surprised therefore that vibration was not a consideration. For what it 
will do , 9000 dollars seems pricey (especially when TREETOPS is in the public domain) 
but it will run on a MAC. Note on page 7 of the enclosure a chart comparing various 
code efficiencies. 


B-2 



International 





6/20^90 J. Munson #2 


> 

a 



5 ° 

53 

If 

s 


O 

CO 


=r h 

o zr 
3 © 

3 0 

p 0 


V w 

0 


o 


© 3 : 
C/> o 

u W 

aa 

s > 

Q. H 

cr O 
< C/D 


7T r-+. 
0 

0 

CO Si 
0 


0 

CO 


o 

o 

5 

0 

0 

0 

0 



B-4 


AUTOSIM was developed by Dr. Mike Sayers at 
the University of Michigan Transportation Research 
Institute. 
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AUTOSiM Technical Summary 
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nTIBOOY COORDINATES 


0(1): Translation of BO ralativa to O, (nil. (a) 
g(2)i Translation of B0 ralativa to o, (n2). (a) 
0(3): Translation of B0 ralativa to o, (n3). (a) 


( rad) 
(rad) 


0(4)1 

Rot. 

0(5): 

Rot. 

0(6): 

Rot. 

0(7): 

Rot. 

0(6): 

Rot. 

0(5): 

Rot. 

;(10): 

Rot. 

‘JOCFCRDBfT : 

U(l): 

Abs. 

U(2): 

Abs . 

3 ) : 

Abs. 

U < 4 > : 

Abs . 

U( 5) : 

Abs. 

U< 6) : 

Abs. 

U(7): 

Rot. 

U(8): 

Rot. 

U(9): 

Rot. 

'J(IO): 

Rot. 


( rad) 


trans. speed of BCM, axis 1. ( a/s > 

trans. speed of BCM, axis 2. (m/s) 

trans. speed of BCM, axis 3. (n/a) 

rot. speed of B about axis #1. (rad/s) 

rot. speed of B about axis #2. (rad/s) 

rot. speed of B about axis 13. (rad/s) 

speed of C relative to B, axis 3. (rad/s) 

speed of D relative to C, axis 1. (rad/s) 

speed of E relative to B, axis 3. (rad/s) 

speed of F relative to E, axis 1. (rad/s) 


OMENTS 


Tl: boo*- torque Z; Expression - -F(l); Direction - (b3). Acts on the 
boo* fro* the bus 


T2: boo*-torque X; Expression * -F(2); Direction « (el). Acts on the 
boo* fro* the bus 

LOCKT: torque fro* clock motor? Expression - F(3); Direction ■ (b3). 
Acts on the clock fro* the bus 


AMT: torque fro* ca»era motor? Expression * r(4); Direction * (cl). 
Acts on the caaera fro* the clock 
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each danvitivi evaluation requires 544 m uitiply/dividaa f 4 72 
idd/aubtraets , and a function/aubroutme calls. 

Calculate trig (unctions of state variables 


S ( 7 ) • SIN(Q( 7) ) 

S(8 ) - S IN ( Q ( 8 ) ) 

C( 7 ) - COS ( Q { 7 ) ) 

C ( 8 ) - COS ( Q ( 8 ) ) 

Kinamatical equations 

2(1) - P( 1 ) *U( 5 ) 

2(2) - P(1)»IH4) 

2(3) « ~ ( U( 2 ) -2(2)) 

QM1) * (Q( 5 ) # U( 3 ) + U(l) 4- Z ( 1 ) 4- Q ( 6 ) *2 ( 3 ) ) 
2(4) m (UC 1 ) 4- Z(l) ) 

QM2) - -<0(4)MJ(3) -U( 2 ) ♦ 2(2) -Q(6)»2(4)) 
Q'(3) ■ -(-U(3) 4- <3 ( 4 ) *2 ( 3 ) 4* <3(5)*Z(4)) 

QM4) - -(Q(6)MJ(5) -U( 4 ) ) 

QM5) - (Q( 6 ) *U( 4 ) 4- U(5) ) 

QM6) - (U(6) -Q(5)*Q' (4) ) 

QM7) » U( 7 ) 

Q'(8) - U( 8 ) 

Q ' (9 ) - U( 9 ) 

QM10) - U(10) 


External aubroutines and extra variable* 

CLKO®, CAMO© ■ CMD ( SYSTEM TIME) 

define expression for boon-torque t 

r<l) - ( KB*Q( 9 ) + B8*U(9)) 

define expression for boon-torque x 

f{7) m ( KB*Q{ 10 ) + BB*U( 10 ) ) 

lefine expression for torque fron clock motor 

M3) * ( KCLOCK * ( CUCCMD -Q( 7) ) -BCLOCK*U( 7 ) ) 

lefine expression for torque fron camera motor 

M4) - ( KCLOCK* (CAMQtD -Q(8)) -BCL0CK»U( 8 ) ) 

<»ftne expression for thruster torque tl 

(5) - LTT1 *TKRUST ( SYSTEMJTIME , 1, fQ(4) + GY?0*U(4)}) 
efine expression for thruster torque #2 

(6) - LTr2*THRUST( SYSTEMJTIME, 2, (Q(5) + GY*0*U(5))) 
efine expression for thruster torque #3 

(7) - LTT3 *THRUST ( S YSTEMJTIME , 3, (Q{6) +GYRO*U(6))) 

/n ami cal equations 

ORIGINAL PAGE IS 
OF POOR QUALITY 
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o 

o 

II 

Qp 16) 

o 

o 

QPI7) 

= 0.8! 

QPt8> 

* -0.! 

CPI9) 

- 0.0 

opuo) 

a 

o 

o 

UP 1 D 

o 

o 

UP! 2) 

o 

o 

If 

UPI 3 ) 

- 0.0 

UP ( 4 ) 

o 

o 

* 

UP 1 5 ) 

- 0.0 

UPI6) 

9 

O 

o 

Upl 7 ) 

M 

O 

o 

UPI 8 ) 

* 0.0 

UP 1 9 ) 

o 

o 

UPI10) 

=• 0.0 

OION initial 



rad, Rot. of Bp ralativa to Bpp about axis #2. 

,! rad. Rot. of B ralativa to Bp about axis #3. 
073464102069 , : rad. Rot. of C ralativa to B about axis 13. 

, i rad, Rot. of D ralativa to C about axis 11. 

,1 rad. Rot. of E ralativa to B about axis *3. 

,l rad, Rot. of P ralativ# to E about axis #1. 

m/s, Aba. trana. spaad of BOI, axis 1. 

,: m/s, Aba. trans. spaad of BQ4, axis 2. 

,1 m/s, Aba, trana. spaad of BCM, axis 3. 

,1 rad/a, Aba. rot. spaad of 8 about axis #1. 

,! rad/s. Aba, rot. spaad of 8 about axis i2. 


Abs. 

rot. spaad 

of B about , 

axis 13. 


Rot. 

spaad of 

C 

ralativa 

to 

B, 

axis 

3 

Rot . 

spaad of 

D 

ralativa 

to 

C, 

axis 

1 

Rot. 

spaad of 

E 

ralativa 

to 

B, 

axis 

3 

Rot. 

spaad of 

r 

ralativa 

to 

E, 

axis 

1 


PRECOMPUTED CONSTANTS 


PCI) - L2*MC/(MB * MCI 
PI 2 ) « IL7 + LB) 

PI 3) - ( ID11 -ID33 ) 

PI4) • < ID22 -ID33 ) 

PIS) - L5 *MD 
PI 6 ) • L6*KD 
PI 7) - (L7 + LB ) *MT 
PI I ) - LB *MT 

PI9)-IMB^MC+MD+MP) 

PI 10) » I rc + MB* | L2 *MC ) * * 2/ 1 MB 4- MC)**2 4- MC*(L2M1 -MC/IMB 4 - 
MC) ) >**2 -f IP1 4- IB11) 

PI1D - IIC + MB* I L2*MC) **2/|MB 4* MC)**2 4- MC*(L2M1 -MC/IMB 4- 
MC) ) ) • * 2 4- IP2 + IB22) 

PU2) • IMP* ( L7 + LB ) **2 4- in 4* IB33) 

PI 13 ) - (IL5*L5 + L6 *L6 ) *MD + ID11) 

P 1 14 ) - (L8ML7 4* LB ) # MP + Ifl) 

PI15) ■ I MT*L8 *L8 + IP1 ) 

PI 16 ) - 1 . 0/PI 1 5 ) 

PI 17 ) ■ PI B ) *P 1 16 ) 

PUB) - PI 8 ) *P| 17 ) 

PI19) - (P(9) -P( 18 ) ) 

PI 20 ) - 1 . 0/P 119) 

PUD - P{8 ) *P( 20 ) 

PI 22 ) - PI 14 ) *P( 16 ) 

PI 23 ) - P( 8 ) *P( 22 ) 

Pf 74 ) • PM 4 ) *Pf 17) 

Pi 25) • Pi 14 ) *P| 22 J 
PI 26 ) - P( 1 ) *P(8 ) 

PI 27) - L8 # P( 8 ) 

PI 28 ) • L7 * P < 8 ) 

PI 29 ) - IIP1 + P( 27 ) ) 

PI 30 ) ■ L7*MF 

PI 31 ) - PI 20 ) *P( 23 } 

PI 32) - IPI12) -PI 25 ) ) 

:D REGION 


NAMXC continuous 
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3(62) 

2(63) 

2(64) 

Zf65> 

2(66) 

2(67) 

2 ( 68 ) 

2(69) 

2(70) 

2(71) 


- -< -Ll * (-Z( 5 ) *Z( 6 ) * 2(401 * Z < 41 1 ) + L5*(Z(9)*Z(14) -Z(43>) 
♦ L6*(Z(S)*Z(14) + Z f 44 ) ) ) 

i ;r;,;, o,9, ‘ o,s ’ + 3 -°* u,io "* z < i7 > ^.io^ms,, 

- rD12*Z<9> 

- ID23 *Z ( 8 ) 

- rD12*Z(14> 

- ; ID ”-2'39, * ID33 *Z ( 4 3 ) 4 1023^44, -Z , 9 , • , -ID22 M , U , + 
Z<64) 4- Z,65)l 4 Z ( 1 4 ) * ( -1011 *Z ( 9 ) 4 Z(66> 4 z, 67M) 

- ( ID23 *Z ( 9 ) + ID1 3*2(14 ) ) 

" zILT*'” 1 * ID33 ‘ Z ' 431 + * Z(I)*(P(3)*Z(14) ♦ 

Z(64) + 2(65) ) -Z(14)*Z(69)> 

’ z!«T Z(39> * I013 ‘ ZM3 > * “>»**««> -Z( 8 ) * ( P( 4 ) *Z ( 9 ) ♦ 
Z(66) + 2(67)) + Z(9 ) *Z{ 69) ) 


2(72) - MD*Z ( 25 ) 
2(73) - MD*Z ( 23 ) 
2(74) - MD*Z ( 21 ) 
2(75) - ID13 *Z ( 10 ) 
2(76) • ID12*Z(U) 
2(77) - 1011 *C( 7) 
2(78) » MD*Z ( 28 ) 
2(79) - MD*Z ( 27 ) 


Z ( 80 ) » MD*Z ( 26 ) 

Z(81| . f-ID23*Z(12) 4 ID22*Z( 13 ) + ID12*S<7)) 

Z(82l - <-IDJ3*Z( 12 ) + ID23*Z ( 13 ) 4 ID13*S(7)I 
2(83) • 1013*2(12) 

Z(84) - ID12*Z( 13 ) 

2(85) - ID11*S(7) 

Z(86) • (-2(83) + 2(84) + 2(85)) 

2(67) - MD*Z { 29 ) 

Z(8I) - ( ID33 *C( 8 ) +• 1023*5(8)) 

2(89) - ( 1023 *C( 8 ) ♦ ID22*S(8)) 

2(90) » ( ID13 *C( 8 ) + 1012*5(8)) 

2(91) • MF*Z ( 32 ) 

2(92) » MT*Z ( 31 ) 

2(93) - MT*Z(30) 

2(94) - MF*Z ( 34 ) 

2(95) - MF*Z(33) 

2(96) - iri*Q(10) 

2(97) - IF1*Q(9) 

2(98) - P( 30 ) *Q(9 ) 

2(99) - XP2 *Q( 10 ) 

ZU00) - -HD*(ZI10)*Z(60) -Z(U)*Z(61) + Z(62)*C(7)J 
Z(101| - MD* (-Z( 12 1 *Z( 60 ) 4 - Zl 13 ) *Z(61 ) 4 - Z(«2)‘S(7)) 

2(102) » MO* ( Z ( 60 ) *C( 8 ) ♦ Z(61)*S(8)) 

Z(lOl) - <r(5) + xc* (U( 5 ) *Z( 5 ) -Z(36)) -Z<10)*Z<68> 4 Z,ll,. 2( 70, 
— Z(62 ) *Z( 72 ) -Z ( SI ) *Z ( 73 ) 4- Z(60)*Z(74) -Z(71)*C(7)) 
ZU04I - -( Z ( 11 ) *Z( 73 ) 4- Z ( 10 } *Z ( 74 ) -0<9)*Z<91) + Z (93> 

— Z( 72 ) *C( 7 ) ) 

21105) - (-Z(13)*Z(73) Z (13) *2(74 ) ♦ Z(91) -<>(10)*Z(92) 

— Z( 72 ) *S( 7 ) ) 

1,1061 - (Q( 10 ) *Z (91 > 4 . Z(92) 4 Z{ 74 ) *C( 8 ) -Z(73)*S(8)» 

2,1071 - (-Z ( 11 ) *Z ( 79 ) -Z( 10 ) *Z ( 80 ) 4- Z(94) -Z(7B)‘C(7)) 

ZI108) . -( Z( 13 ) *Z ( 79 ) 4- Z( 12 ) *Z ( 80 ) -Z(93) 4- Q(9)*Z(94) 

-Z( 78 ) *S( 7 ) ) 


Z(109) - ( Z ( 95 ) 4- Zf 80 ) *C( 8 ) -Z(79)*S(8)) 

Z(llO) . (IB12 -IF2*Q(9) -Z(25)*Z(78) 4- Z(23)*Z(79> 4 Z(21).z(80) 
-Z(11)*Z(81) 4 - Z( 10 ) *Z( 82 ) 4 - Z(32)*Z(93) -Z(30)*Z(94) 4 - 
Z(31 ) *2(95 ) + Z ( 97) + 2(86) *C( 7) ) 

2(111) - 2 ( 62 ) *Z ( 87 ) 
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Z(51 - <U(6) ♦ U ( 7 ) ) 

Z(6) - <0(4>*C<7) ♦ U< 5 ) *S ( 7) ) 

Z ( 7 ) - (U( 5 ) *C( 7) -0( 4 ) *S ( 7 ) > 

Z ( 8 ) - <Z(5)*C(8> -Z ( 7 ) *3 ( 8 ) ) 

2(9) - (Z(7)*C(8> + Z { 5 ) *S { 8 ) ) 

2(10) - S { 7 ) *S ( 8 ) 

Z( 11 ) - C( 8 ) *S( 7) 

Z ( 12 ) - C I 7 ) *S f 8 ) 

2(13) • C ( 7 ) *C ( 8 ) 

2(14) - ( U( 8 ) ♦ 2(6) ) 

2(15) • (U(6) + U(9)t 
2(16) - (Q(9)*0(4) -U { 5 ) ) 

2(17) • (Q( 10 ) *Z ( 15 ) -Z ( 16 ) ) 

Z ( 18 ) • P ( 1 ) *S ( 7 ) 

Z ( 19 ) - Pf 1 ) *Cf 7) 

2 ( 20 > - <L5 -L3*C( 8 ) ♦ U*S{8)) 

2(21) - < 2 ( 20 ) *C( 7 ) -Z(19)*S(B) J 
2(22) - L3*S ( 8 ) 

2(23) - (C(7)*tL6 ♦ 2(22) + L1*C(8)) -2 ( 19 ) *C( 8 ) ) 

2(24) - LI *5 ( 7 ) 

2(25) =* ( L5 *2 ( 10 ) + L6 * Z ( 1 1 ) -2(18) + 2(24) ) 

3(261 - I 2( 20 ) *S ( 7 ) -2 ( 18 ) *5 ( 8 ) } 

2(27) » ( (-2(16) + Z ( 24 } ) *C( 8 ) + ( L6 4- Z(22))*S<7>) 

2(28) - ( L5 * 2 ( 1 2 ) ♦ L6 *2 ( 13 ) -2< 19 ) + Li‘C(7)) 

Z ( 29 ) - (U -L5*C( 8 ) 4- L6 *3 < 8 ) ) 

2(30) - P( l ) *Q( 9 ) 

2(31) - ( PC 2 ) -P( 1 ) *0( 10 ) ) 

2(32) - (PCD + L7*Q(10J) 

2(33) - L8*Q<9) 

2(34) - (PCD -L8*Q(10) ) 

2(35) - 0(4)*0(7) 

2(36) - U( 5) *0(7} 

2(37) - 2(36) *C( 7 ) 

2(38) - Z<35)*3(7) 

2(39) - (2(37) —2 (36)) 

2(40) - 2 ( 35 ) *C( 7 ) 

2(41) - Z< 36 > *3 ( 7 ) 

2(42) - (2(40) 4- Z(41)J 

2(43) « -( 0( 8 ) *Z(9 ) —2(421*3(8)) 

2(44) - ( U( 8 ) *2( 8 ) -2(42) *C(8 ) ) 

2(45) - U( 4 ) *U< 9 ) 

2(46) - -( 0( 5) *0(9) -Q{9 ) *Z(45) ) 

2(47) - U( 6 ) *0( 4 ) 

2(48) - U( 4 ) *0( 5 ) 

2(49) - U( 5) *U(6 ) 

2(50 ) - U( 4 ) *U( 4 ) 

2(51 ) - ( Z ( 50 ) 4* 0(6) *0(6 ) ) 

2(52) - -(0(6) *0(2 ) -0 ( 5 ) *0 ( 3 ) ) 

2(53) - -(0(4) *0( 3 ) -0(6) *0( 1 ) ) 

2(54) - Z(5)*2(7) 

Z( 55) - (2(37) -2(38) 4- Z(54)) 

2(56) - 2 ( 8 ) *2 ( 9 ) 

2(57) - 2 ( 14 ) *Z( 14 ) 

Z(58) - -LI *(-2(37) +• 2(38) + 2(54)) 

2(59) - L3 *2 ( 5 ) *2 ( 5 ) 

2(60) - -< L5* (Z( 39 ) 4 Z ( 56 ) ) -L6*(2(57) 4- Z(9)*2(9)) -L3 *2( 55 ) *C( 8 ) 
4 (2(58) -2(59) ) *S<8) ) 

2(61) - (L6*<2(39) -Z { 56 ) > 4- L5*(2(57) + 2(8)*Z(8)) + (Z(58) 

-2(59)) *C( 8 ) + L3*2(55)*S(8) ) 
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:ii67i =. mini + z<ii5>»2<i6i n/2ti6«i 
11168) - (ZOO) -Zdl5>*2(162) I/ZI166) 
2(169) - ( Z( 135) * 2(11502(163) )/Z(166) 
ZI170) - (2(133) 4 2(11502(164) ) /Z ( 166 ) 
2(171) - ( z ( 134 ) 4 2(115)*Z(165) > /Z (166) 
2(172) » P( 1 7 ) ft Z ( 14 7 ) 


Z( 173 ) • P ( 1 7 ) *Z ( 14 3 ) 


Zfl74) - (Z(117) + Z(160)*Z(173) ) 

Z(175) - (P(9> -Z(167) *2(174) ) 

2(176) - ( Z 1 138 ) -Z ( 154 ) *Z< 172 ) -Z( 162 ) *Z( 173 ) 4 
2(164 02(174 ) )/Z( 175) 

Z(177) - ( Z ( 118 ) * P< 22 02 ( 143 ) -2(16302(173) 
-2(169 02(174) )/Z(175) 


2(178) - ( Z(10S) * Z( 147 02(151 ) -z (155 ) *2 , 172 ) -Z ( 164 ) *2 ( 173 ) 

-Z( 170 ) *2(174 ) )/Z( 175) 

2(179) - (2(108 ) * 2(143 02 ( 159 ) 4 Z ( 136 ) *2 , 172 ) -2,165 0 2 0 73) 
-3 (171 OZ ,174)) /2 ,175) 

Z ( 180 1 - (ZOO) -2(1371*2(160)) 

Z(181> - (2(138) -Z(139J*Z(153) -Z ( 137 ) *2 ( 161 ) ♦ Z( 167) *2(180) ) 
Z(182) - (P(13) -2( 139 ) *Z ( 154 I -2(1371*2(162) -2 ( 168 ) *2 ( !80 ) 

-Z( 176 ) *Z( 181 ) ) 

Z(183) - (2(90) 4 Z ( 137) *2(163 ) -Z ( 169 ) *2 ( 180 ) + 

Z ( 177 ) *Z ( 181 ) )/Z (182) 

2(184) - (2(140) 4 2(1391*2(155. 4 2(1371*2,164) -2 ( 170 ) *2, l 80 ) * 
2(1781*2(181) )/Z( 182 ) 


2(185) - ,2(141) -2(139. *2, 136) 4 Z( 137 , *2, 165 ) -2, 171 ) *2( 180 ) * 

3(179 ) *2(181 ) )/Z (182 ) 

2(186) - ( — P ( 24 3 4 Z( 116 ) ) 

2(187) m (2(135) 4 Z( 160 1*2(186) ) 

2(188. - ,2,118) ♦ P,14 ) *2, 157 ) -2 , 161 ) *2 ( 186 . -2,1671*2(187), 
2(189, - ,200) 4 2(162, *2, 186, -2 ( 168 ) *2, 187 ) ♦ 2(176 , *2,188 , , 
2(190) - (P(32) 4 ZI129, 4 2(130) 4 2(131, -Z( 163 ) *z ( 186 I 
-Z(169)*2(187, -Z ( 177 > *Z( 188 ) -2(1831*2(189)) 

2(191) - (2,123) 4 P( 14) *2(158 ) -Z( 164 ) *2(186 ) -2,170 ) *2, 187 ) 
-2(178) *2(188) -2 (184) *2 (189) 1/2(190) 

2(192, - (2(128) -P( 14 ) *2(159 ) -2 ( 165 ) *2 ( 186 ) -2(1711*2(187) 

-Z( 179 ) *2(188 ) -Z(185)*Z(189))/Z(190) 

2(193) ■ ( Z ( 106 ) -P{ 17) *Z(148 ) ) 

2(194) - ( Z ( 104 ) + P( 17) *Z( 144 ) ) 

2(195) - ( Z ( 133 ) 4* Z ( 160 ) *Z ( 194 ) ) 

2(196) - (2(105) 4 ZC148 1*2(150) -Z ( 153 ) *2 ( 193 ) -Z ( 161 ) *2 ( 194 , 

-2( 167) »2( 195 ) ) 


2(197) - (2(140) 4 Z(154) *2(193) 4 2(162) *2(194) -2(168) *2(195) 4 
2(176)»2(196) ) 


2(198) - (Zd23) 4 P< 22) *2(144 ) -2(163) *2(194 ) -2,1691*2,195) 

-2(177 ) *2(196 ) -Z(183)*2(197)) 

2(199) - (P( 10) 4 2( 25 ) *2( 72 ) 4 Z(23)*2(73) 4 2(21)*2(74) 4 

2,321*2(91) + Z, 311*2(92) -2(1481*2(151) -21 155) *2, 193 ) 
-2(1641*2(194) -2(1701*2(195) -Z ( 178 ) *2 ( 196 ) -2 ( 184 ) *2 ( 197 ) 
-2(1911*2,198, * (2(75) -2(76) 4 Z(77)|‘C(7) 4 2(10) 
*(1033*2(10) -ID23*2( 11 ) 4 1D13*C(7)) -2 ( 11 ) * ( ID23 *2 ( 10 ) 
-rD22*Z(ll) ♦ IDI2*C( 7 ) ) ) 

2(200) - (ZI110) -2(1481*2(152) 4 Z ( 144 ) *Z( 159 ) 4 2 ( 156 ) *2 ( 193 ) 

-Z( 165 ) *2, 194 ) -2(171 ) *2(195) -2(179 ) *2,196 ) -2(1851*2(197) 
-ZU92)*Z(198))/Z{199) 

Z{201) • (-2(109) 4* P ( 17 ) *Z(149 ) ) 

2(202) - ( Z ( 107 ) -P( 17 ) *Z( 145 ) ) 

Z(203) - ( Z ( 134 ) + Z (160 ) *2(202 ) ) 

2(204) « (ZOOS) 4 2 ( 145 ) *2 ( 157 ) 4 2 (153 ) *2 ( 201 ) -2(16102(202) 
-Z(167)*2(203) ) 
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:U12) a 3(68 ) *C(8 ) 

2(113) - 2 ( 70 ) *S ( 8 ) 

2(114) - (F<7) 2(111) -Z ( 11 2 ) -Zf 113) ) 

2(115) - Z ( 8 7 ) *C( 7 ) 

2(116) - (P(7) -Z( 115 1 ) 

2(117) » 2 ( 87 ) *S ( 7 ) 

2(118) ■ ( -P ( 7 ) •<) ( 9 ) + Z ( 98 ) + 3(117} ) 

2(119) - Z(25)*Z(07) 

2(120) - Z ( 10 ) *2 ( 68 ) 

2(121) - 2(11) *2(69) 

2(122) • 2 ( 90 ) *C( 7 ) 

2,1231 - IIB13 -P(7) *2(30) + Z(32CZ(98) -Z( 119 ) ♦ Z(120S -ZU21) ♦ 
Z(122 ) ) 

2(124) - Z ( 28 ) *2 ( 87 ) 

2(125) - Z ( 1 2 ) *Z( 88 ) 

2(126) « 2 ( 1 3 ) *Z( 09 ) 

2 1 127) - Z ( 90 ) *S ( 7 ) 

2(128) - (IB23 + P( 7) *2( 34 ) -Z(96) 4- Z{99) 4* 2(124) -2(125) 4- 2(126) 
♦ 2(127) ) 

2(129) - Z f 29 ) *2 ( 67 ) 

2(130) - 2 ( 88 ) *C{ 8 ) 

2(131) « Z ( 89 ) *S ( 8 ) 

2(132) - <F(3) ♦ Z ( 111 ) -2(112) -2(113) ) 

3(133} * < -Z ( 119 ) + 3(120) -Z(121) + Z(122)) 

ZI134) - ( Z( 124 ) -Z ( 125 ) + 2(126) + Z(127)) 

2(135) - ( Z ( 129 ) ♦ 2(130) ♦ Z < 1 3 1 ) ) 

2(136) • (F(4) 4- P( 5 ) *Z( 60 ) -P(6)*2(61) -2(71)) 

2(137) • ( P( 5 ) *Z{ 10 ) ♦ P( 6 ) *2 ( 11 ) ) 

2(138) - ( P( 5 ) *Z( 12 ) + P( 6 ) *2 ( 13 ) ) 

Z(139) « -( P( 5 ) *C( 8 ) -P(6 ) *S(0 ) ) 

2(140} - f P( 5 ) *2 ( 21 J + P(6 ) *Z( 23 ) 4. 2(75) -2(76) 4 2(77)) 

2(141) - < Pf 5 ) *Zf 26 ) + P(6)*2(27) -Z(83) + 2(84) 4- 2(85)) 

2(142) • -(-f(l) -P( 26 ) *2(47) + P(28)*Z(40) -P(8)*(Z(52) + 

0(9 ) *2( 33 ) ) + P( 29 ) *2 ( 63 ) ) 

2(143 ) - P< 0 ) *0(9 ) 

2(144) - ( P(8 ) *Z( 30 > 0(9) *2(99 ) ) 
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U * f 5 ) 

* 2(236) 

U' (A ) 

- 2(237) 

UM6) 

- 2(238) 

LT(0) 

- 2(239) 

UM2) 

- Z ( 240 ) 

U' (7) 

- 2(241) 

U'(l> 

* 2(242) 

U' (9) 

- - ( P ( 22 ) *2 ( 238 ) 


2( 144 ) *Z( 237) 

U'(3) 

- 2(243) 

U' (10) 

■ -( P( 16 ) * ( Z ( 1 46 ; 


2 ( 147 ) *Z ( 240 ) ) 


(237) 


>>) 


‘TO DYNAMIC 
<ECUTE 'CONTINUOUS' 

•••* cod. .bov. thi. an, g.n.r.t.d ^ Atm5SIM .... 

JNCTION cmd (cikcmd, camcmd ■ t) 

TP t .it. 1 THEN 
cikcmd m 4 - .pi, 
camcmd ■ -0.5 
ELSEir t .It. II THEN 

cikcmd - 4 - 0.025 * (t - 1 ) - Pr# 
camcmd * -0.5 + 0.01 * (t-l) 

ELSE 

cikcmd - 3.75 - .pi. 
camcmd ■ - 0.4 
ENDIF 

ND FUNCTION 


UNCTION thrust (thr - t, sxis, tsrror) 

ECTOR flr,( 3) , toff(J) 

\TA dband • 0.0025 
VTA tmin * 0.02 

\TA flra - (0,0,0), toff * (0,0,0) 

rr tarror .It. -dband THEN 
firafaxis) - l 
toff (axis) - t + tmin 

ELSEIf tarror .gt. dband THEN 
fira(axis) ■ -l 
toff (axis) » t + tmin 

ELSEIF t . ga . toff (axis) THEN 
f ira( axis ) » o 

ENDIF 

thr m fira(axis) 

0 FUNCTION 
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The "Workshop on NASA Workstation Technology" held in Palo Alto, California, during 
the week of March 1 2, 1990 was attended. This work is sponsored by NASA Code R 
(Lee Holcomb). It is centered in the Research Institute for Advanced Computational 
Science (RIACS). This is a NASA funded organization formed by USRA "to provide 
leadership in developing computers” and is headquartered at NASA AMES Research 
Center. The Workshop was emceed by Dr. Bob Brown, RIACS’ technical director, 
acting under the overall direction of RIACS’ assistant director, Dr. Berry Leiner (Dr. 
Brown's boss). The meetings were actually held at the Hyatt Rickeys. The meeting 
room was certainly adequate, the noon meals, served buffet styles in an adjoining dining 
room were outstanding as were the break snacks and the "heavy hors d' oeuvres" 
served on the second evening during the vendor show. 


The Workshop, while exceedingly interesting from the viewpoint of being exposed to 
new vistas of thought, was in no wise along the lines we had discussed before my 
going. And at no time did anyone define a Workstation. I will offer my own definition for 
discussion later. Even at a distance of two weeks it is still not entirely clear what 
electicism was used to choose the variegated program. I will include a brief summary of 
my impressions of most of the talks but perhaps a synopsis of the keynote address will 
provide a general "sense of the meeting." 


The keynote address was given by Dr. William Bricken, Chief Scientist, Human Interface 
Technology Center, University of Washington. In his talk Dr. Bricken presented a 
"Vision of Virtual Reality." In his address he suggested the use of the computer (and an 
appropriate human sensory environment) as a replacement for reality. One would, in 
this scheme, seem to be "wearing a computer." In this scenario anything imaginable 
would be sensorally possible, e.g., the known laws of physis would be irrelevant. It 
would be possible as an example to sense being on Mars with a population of penguins 
moving mountains with the flick of a (sensed) finger. A pale subset of this direction of 
sensing is seen in the very realistic scene generation, flying maneuvers etc currently 
available in USAF training commands. Of course, he pointed out that carrying this work 
forward will require new paradigms of computing and boundary mathematics. Some of 
the audiences' comments and discussions at the break regarding this address included 
a few catch phrases, e.g., post STARSHIP, sanity, realism, "electronic crack." My own 
reaction is that with time much of this could be possible and we will then pass through 
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another door into yet another Brave New World. Stopping short of the ultimate there 

seems to be considerable use for this sort of thing in education, training, entertainment, 
and repair work. 


Mr. Lee Holcomb gave the charge to the group. A group of us also talked to him during 
one of the breaks. The subject during the break was the upcoming national High 
Performance Computing Initiative. Mr. Holcomb had attended a meeting not long ago 
chaired by the President's science advisor (Dr. Bromley) and attended by OMB, NSF, 
NASA, and DOD representatives. The information was that the advisor had been 
roughly handled in Congress because of the lack of an Administration's High 
Performance Computing Initiative plan and the subsequent lack of any budget for it. 
This will be rectified in the next ('92) budget and this was agreed to by the OMB 
representative present. In his charge to the Workshop Mr. Holcomb listed three 
applications of advanced computing. The first was its use in operations. This included, 
as examples, both space and aeronautical operations to include real time control 
(telescience), anomaly identification and correction and scheduling and planning of 
operations. 


His second application was that of science analysis. There are almost unimaginably 
large volumes of data currently available. These data are being added to daily and as 
projects such as the Space Station come on line an even more rapid stream of data will 
occur to fill the data banks. The challenges then include finding ways to allow 
appropriate people access to these huge data banks, how to couple various disciplines 
(e.g. plankson biology and oceanography) to produce coupled phenomenalogical 
models and just in general how to facilitate the collaborations of multidisciplinary 
scientists on a global scale. His third and last charge (and as far as I could see the only 
explicit mention of the totality suggested) was to address how to build complex 
aerospace systems. He mentioned that computational fluid dynamics, structural design 
etc are being pursued but that the overall challenge is to pull together or model the 
interaction of fluid flow, aerodynamics, guidance and control, structures, stealth and so 
forth so as to produce an integrated, optimal design - or, in the case of something like 
NASP, perhaps any viable design at all. He concluded expressing hope that the 
Workshop will encourage an exchange of general information, technology, applications, 
capabilities vendors and software. He also asked rhetorically what standards should 
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NASA support? For instance networking and the UNIX versus X-WINDOX 
versus. ..discussion. He also expressed concern over high maintenance cost e.g. 25% 
of DEC'S income is from maintenance work. 


There is supposed to be a mailing of hard copy of the presentation materials. A video 
camera was run during the presentations but as we left we were told that the quality of 
the video was very poor and only the audio was useable. In all events the media 
presentation materials were generally either in viewgraph or VCR form. Therefore, I 
shall set forth below some of my impressions of the various speakers' messages. 


Mr. Jay Costenbader of GSFC, mission operations, picked out three workstations 
usages at GSFC. These were in mission operations, data capture and guidance and 
control. My impression is that although there are many automation aids via the 
computers involved there is much human intervention and decision making. He 
indicated that this is because of the very conservative approach taken by the various 
system (e.g. Hubble ST) managers. They have looked at the ADA language and found 
it too restrictive, they (as are a number of others) are moving on to C++. 


Mr. Jim Jeletic of GSFC, Flight Dynamics Division, reinforced a number of Mr. 
Castenbader's points and added a few. To wit, they have to support long term ongoing 
programs with current (i.e. may be very old) inflight protocalls (may go on for years). 
This makes for a mix of operational systems to support. He sees more ethernet 
"bridges" and more distributed computing in the future. He pointed out that there is an 
ADA champion at GSFC thus ADA is used whenever possible. He sees concerns 
including resource consumption (memory, performance), system compatibility, system 
complexity (which workstations to use, what operating system) and distributed systems 
making for harder troubleshooting. 


Mr. Mike Wiskerchen of Stanford University had a most interesting story to tell. Under 
Code R support he undertook to study the application of MIS to NASA operational 
needs. The vehicle for this study was the process of making sure the Orbiter tiles were 
ready for flight. Apparently he had worked for NASA at one time. Anyway, he went to 


C-4 



KSC and found that for the tiles only 360,000 sets of initials or signatures are required to 
release an Orbiter for flight. He found an almost absolute lack of MIS use in the system 
and as he progressed through the system (he actually performed himself all the jobs in 
the chain, it took a number of months) he developed some interesting statistics. These 
included his 85-12-3 rule. This rule is that 85% of the people favor the status quo, 12% 
actively oppose change, and the 3% support change in the modus operandi. He found 
that 90% of the people at KSC do not use computers. Thus, technology is 1 0% of the 
problems and people 90% of the problem. Only the direct involvement of top 
management (who do not know about computers now) and the establishment of 
rewards for technology change (such rewards do not now exist) will cause a change. 
On the technical side there is a problem. The technology of MIS changes on a 6 to 9 
month lifecycle. Programs evolve over 10-15 years . Can you design to replace 
technology every year or so? 


Mr. Randy Davis of the University of Colorado which apparently has a mission 
operations center for some scientific missions gave a presentation on OASIS. OASIS is 
their software package (Operations and Science Instrument Support Teleoperations 
support software). Overall they propose moving control of payloads to the Pi's own 
home premises (versus current gather-everyone-at-a-NASA-center). Their software is 
written in ADA and used on DEC/VAX, Sun Ill's, (S IV upcoming). The general theme 
was to work on a coherent architecture for payload control. 


Mr. Contenbader of GSFC returned and talked about TAE + (Transportable Application 
Environment Plus). This environment separates the users interface from applications 
programmers. It now has spread, after a slow start, to 500 users and in available from 
COSMIC. Some technical requirements were given. 


Mr. Andrew Potter of MSFC (I cannot find him in the October 1989 MSFC phone book) 
gave a talk about the goals to be obtained in designing a workstation for the Space 
Station. It emphasized instruments on a screen ('a la the "glass cockpit") and talked 
about search strategies. His presentations was not on MSFC viewgraphs. 


C-5 



Mr. Tim Castellano of AMES RC talked about the workstation aspects of the SHOOT 
project. This project is supposed to fly in the fall of 1992 on STS61. It is a superfluid 
helium transfer flight demonstration. It includes doing some control and observations of 
the experiment on MAC I Is and a laptop from the POCC at GSFC. It is an expert based 
system and has 1 00 to 1 25 rules so far. 


Ms. Patricia Liggett of JPL aggressively headed a trio of presenters detailing the 
VANESSA project. It was a system rather hurriedly put together from existing 
equipment for the Voyager Neptune Encounter Science Support Activity. In a nutshell it 
allowed the delivery of both engineering and encounter science data in near real time. 
Of course the Voyager data stream bit rate is quite low so equipment and software 
speed was not a problem. The problem was to fan out the data from the data stream 
(not the data base, once in there retrieval is arcane and slow) to a multiplicity of user 
data formats. JPL furnished some of the on site (at JPL) workstations other PI'S 
brought their own. This was a one time effort but showed what type of data display etc 
will be required in the future. 


Mr. Graius Martin of JPL’s Al group (including Cal Tech's Doyle) described the SHARP 
system development status. SHARP stands for Spacecraft Health Automated 
Reasoning Prototype. The objective of this development is to shrink the "army on the 
ground" which is currently seen to be proliferating with the ever increasing number of 
spacecraft borne experiments. 


The goals include an 80% staff reduction, increased safety because of more data being 
under surveillance, more reliability because of system wide monitoring and vastly 
increased productivity because one operator can be spread across a number of flight 
systems through the use of Al. They have put together a Symbolics based workstation 
and later port to a Sun Workstation to support the Magellan mission spacecraft. Various 
aspects are being checked out with the archived Voyager data. 

Mr. Tom Engler of MSFC (can not find him in the phone book either) gave an overview 
of what MSFC hopes to do in automating the HOSC. HOSC is currently actually four 
centers. One each supports Shuttle, HST, Spacelab (POCC) and Space Station. The 
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goal is to multiplex people, use emerging technologies, use Al and possibly neural 
networks etc. Very preliminary design or conceptual goal type of paper. Questions to 
him were mainly had he looked at the other existing control centers. The answer was 
yes, at GSFC, JPL, and JSC. 


Mr. Jerry Blackburn of JSC covered their work toward updating the shuttle through the 
development of MPAC (Multipurpose Applications Console). This would represent an 
update to 386 technology, use ADA and a 1553B bus. They are also developing Space 
Station DMS and have MacDac and IBM "on the team." 


Mr. Eric Hibbard of AMES RC reported on the letting of 4 seven year contracts to DEC, 
Silicon Graphics, Sun, and Convex concerning workstations and midrange computers. 
They are to develop graphics standards in software and hardware. They are looking 
toward animated stereographies, distributed graphics (i.e. created from various non co- 
located sites) and various rendering tools. They now are using X-windows system 
(release 4) and advocate CGI extensions. 


Mr. Val Watson of AMES RC had a number of recommendations for improvements in 
Scientific Workstations. See his handout, a copy is enclosed. 


Mr. Keith Lane "nationally well known for multimedia" who just left Olivetti Computer 
(80% force reduction) talked on a portable work environment (which is mu just a 
portable computer). The gist is that your portable computer is a window to your 
resources (which may be anywhere). He suggested audio at the 768 Kbyte to 2 M Byte 
level, video at 90-240 or 1 0-60 (compressed) and data up to video rates. It seemed to 
me that all this was very much in the hand wavy hobby shop stage. 


Mr. Barry Leiner, Assistant Director of RIACS, described the National Collaboratory 
(Collaborative + Laboratory). Basically this idea grew from a conference held at 
Rockefeller University in 1989. A report of that conference is available. It is basically 
another attempt to deal with world wide knowledge/data explosion. Its tools would 
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include a digital library (Knowbots, knowledge robots, to include digitizing the world's 
knowledge), interoperable data, information fusion, smart/remote instruments, joint 
authority tools and symbolic assistants (the theorist's workbench e.g. mathematical 
objects versus numbers). NSF has a new initiative in this direction; 3 AOs are for this 
sort of thing. 


Mr. Joe Hale of MSFC (EC23) talked about the development surrounding the use of the 
Private Eye. This is a static data display instrument which is worn over (or slightly 
under) one eye and produces a virtual image focused from 10 inches to infinity. It may 
or is being combined with a Dataglove to provide Space Station documentation which 
goes along with an operator (versus using a "floating in zero -G" manual). While in early 
state of development seems promising to me. 


Mr. Michael McGreevy of AMES RC put on a razzle dazzle picture show of various 
geological data. His conclusion was that desktop data presentations are not adequate. 
Something like the head mounted virtual reality techniques are needed (he buttressed 
this thesis by showing mosiacs from the Moon and Mars). 


Because no one else offered their notion of what constitutes or defines a workstation I 
will, in the spirit of polemics, offer one. I propose to define one as a single person 
communication node that has access to outside information resources in addition to its 
own capabilities and storage means. 


My conclusion as far as this NASA effort is concerned is that it is a "going concern" what 
with RIACS off and running. AMES clearly has a large interest and influence in the 
effort e.g. RIACS is housed there. There is a very strong science flavor to the effort. 
This science flavor manifests itself in several ways. One way is in the pure development 
of the various leading edge technologies involved e.g. A.I. Another is in the research 
and application interests of the present group. This includes planetary exploration, the 
processing and display of very large scientific data bases e.g. the portrayal of 3D 
turbulent flow fields created by CFD. The other type of use growing out of the programs 
(e.g. Space Station) is in operations. Examples include mission operations (reduce the 
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Army on the Ground) and training about and repair of orbital hardware. The direct 
application of all this to current known ED-12 programs is not obvious. In the long run 
faster running generally more capable workstation technology will benefit LSS control, of 
course. Based on this conference little or no direct hardware in the loop control 
problems are being addressed specifically (at least of a type applicable to say CASES). 
Probably monitoring progress reports would be worth the effort especially in real time 
data collection disciplines. To try and do any serious steering to the direction of the 
effort would be difficult and require dealing with the entrenched interests of AMES, JPL, 
and RIACS. 
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Workshop on NASA Workstation Technology 

Agenda 


Dav One - March 13, 1990 

8:30 AM Registration desk opens. 

Coffee and pastries served in the Ballroom lobby. 

9:00 AM Welcoming and opening remarks. 

Barry Lciner, Peter Denning, and Robert Brown 

9: 15 AM Charge and goals for the workshop. 

Lee Holcomb 


9:30 AM 


10:10 AM 


12:00 PM 
1:00 PM 


3:00 PM 
3:15 PM 


5:00 PM 
6:30 PM 


Keynote speaker. „ 

William Bricken, Chief Scientist, Human Interface Technology -enter, univer- 
sity of Washington 
Mission Operations Workstations. 

Jay Costenbadcr, GSFC, Three Mission Operations Workstation Projects 
T im Jeletic, GSFC, Workstation Technology Used for Flight Dynamics Mission 
Support 

Mike Wiskerchen, Stanford University, Application of Advanced Workstation 
Technology to Shuttle Operations 

Break for lunch. 

Location: Ballroom C 

Application Development Environments 

Eric Hardy, CMU, The SEI User Interface Project 

Randy Davis, University of Colorado, Oasis : Present and Future 

Jay Costenbader, GFSC, TAE+ 

Coffee break. 

Andrew Potter, MSFC, The SoftPanel Prototype 

Mission Science Workstations. 

Tim Castellano, ARC/RIA, SHOOT 
Patricia Liggett, JPL, The VNESSA System 

Technical sessions end. 

Vendor show, poster and demonstration session. 

Hors d' oeuvres and no -host bar 
Ballroom C &.D 


9:30 PM Day one ends. 
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Dav Two 
8:00 AM 

8:30 AM 

10:30 AM 
10:45 PM 

12:00 PM 
1:00 PM 

3:00 PM 
3:15 PM 

5:00 PM 
5:15 PM 


March 14. 1990 

Coffee and pastries. 

Ballroom lobby 

Mission Operations (continued) 

Gaius Martin, JPL, Sharp 

T om Engl er, MSFC, MSFQ Workstation Lab:SSF and AXF 

Gregory Blackburn, JSC, SSF Multipurpose. Applications Console (MPAC) 

Coffee break. 

Science Data Visualization. 

Eric Hibbard, ARC, Code RCD Graphics 

Leo Plume, JPL, Linkwinds — A Prototype Scientific Visualiztation System 

Break for lunch. 

Location: Ballroom A 

Val Watson, ARC, Workstation Applications Office Projects 
Stephen Coles, JPL, EASE: An Engineering Analysis Subsystem Environment 
for Real-Time Spacecraft Control 

Productivity and Collaboration Tools. 

Keith Lantz, Consultant, Multimedia Workstations for Collaboration 
Coffee break. 

Barry Leiner, RIACS, The National Collaborator 

Other Workstation Projects. 

Larry Foster, Boeing, TMIS, Micro! hostllntegration 
Michael McGreevy, ARC, Ames Virtual Reality Projects 

Workshop wrap-up. 

Robert Brown 

Break for the day. 


c-u 



-3- 


nav Three - March 15. 1990 

The third day of the workshop consists of working sessions at NASA Ames 
Research Center. We are interested in having 15-25 people work in a small number 
of groups to capture and document the important parts of the workshop and to begin 
the process of creating the workshop report. 

Key elements of the report could be: 

The state of workstation technology inside NASA. 

The priorities for future activities. 

A structure will be created in which the third-day participants will work. This 
structure will be reviewed in a plenary session at the beginning of the day. Then, 
small groups will form and work on a detailed outline of the report. 

Who should participate? Speakers, managers, those with opinions or an interest in 
impacting the report These people will become reviewers for the report as it is being 
integrated and edited over the course of the following two months. 


9:00 AM 

Description of the report-writing process. 
Robert Brown 

10:00 AM 

Break into groups by session. 

12:00 PM 

Break for lunch. 

2:00 PM 

Plenary session for presentation of results. 

5:00 PM 

Workshop ends. 
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Support Personnel 

-- Bob Brown, Program, Demonstrations 

- Dee Doyle, General Scheduling & Management 

- Lorraine Fisher, Local arrangements, messages 

- Others (Green Badges) 

Posterboard Session 

-- Posters hold 9 (portrait) or 6 (landscape) letter-sized pages 
-- Give materials to Lauri Howard or Bob Brown 

- Posters will be presented at the vendor show 

Demonstrations 

-- Limited demonstration capability 
-- Sun 3/160 color & monochrome 
-- SGI IRIS (personal & 4D/220) 

- NeXT 

-- Macintosh (MacOS, MacX) 

- IBM RS6000 

- Sun SPARCStation & Sun 4/330 

- DEC DECStation 3100 

Vendor Session 

-- 6:30-9:30PM, adjoining ballroom 
-- Hors d’oeurves served in lobby 
-- No-host bar 

- Demonstrations mostly on vendor equipment 
-- Video tapes running in the corners 
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Vtl Wtoon 

NASA A me* Rc*e«ch Center 


Outline 

Workstation features important for scientific 
applications 

Illustrations of scientific analysis using 
visualization 

Current capabilities for analysis 
Comparison of current capabilities with "ideal" 
Recommendations for the future 
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Benniremen * t QL Rendering the Representations 

All currently popular scenes can be created with 

• 3D surfaces 

• 3D lines 
- points 

Surfaces composed of 4 sided n on-planar polygons 
(flat shading sometimes preferred) 

Typical scene represented with 10,000 of these 
10 "frames/sec" needed to understand dynamics 

Therefore, the derived requirement is 
100 000 polv eons/ sec with hidden surfaces removed 


Current N AS Workstation Procurement 

300 workstations over 3 year period 
Graphics performance is most important factor 
Threshold of pain for cost approximately S100K 
Most weight given to graphics performance 

Performance measured with our benchmarks 

- surface and line graphics benchmark 

- fluid dynamic code floating point benchmark 


NAS firanhics Benchmark 

consists of 3D scene of Shuttle with particle traces 
and function mapped surface 

measures 3D coordinate transformation rate 

measures polygon and line rendering rate 

- with hidden surface removed 

- with mostly 4 sided non -coplan ar polygons 

- with and without smooth shading 

- with and without transparency 

- with and without antialiasing 


Distribution of Computing Power 

after the Procurement 


supercomputers 

□ 



Current Work on Software 

Combination of old modules into a single program 
Changes to take advantage of new workstations 
Changes to make visualization more interactive 


Visualization Software 


Simulation 

Visualization and Recording 


Scene 

Creation 

Scene 

Viewing 

Animation 

Sequence 

Creation 

Recording 
on Him 
and Video 

0 — PLO 

*— su 

Rl 

T3D — l 
RF — 1 
P — 1 

> 

- HAQ _ 


— - Ctl 

UAj 



L 
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Ad ditional Hardwar e Improvements 
Recommended 

An increase in the field of view of die display to 
more nearly match the eye's field of view 

Effective six-degree of freedom controls 
Sound 

Voice recognition 

- discrete word or phrase now 

- continuous voice later 
Specialized hardware to support the tasks 

described in the next viewgraph 


Additions! Software Improvements. 

Recommended 

Software for new data transformations, filters and 
comparators that will be useful in fluid dynamics 
Software to extract the essence of the fluid flows 
Software to clearly display the essence of the fluid 
flows (to create displays highly tuned to the 
human's cognitive capabilities) 

Software to generate pictures for monitoring 
accuracy and stability of the computer 
sinrm tarinns using the human's pattern 
recognition capa bi l itie s 


r.irrpnt Visualiza tion Research 

Scott Fisher’s "Virtual Environment" project 

- Helmet display fills the field of view 

- Display corresponds to head position 

- Voice recognition 

- Data gloves for manual control of objects 
(scientist feels like he is inside the simulated field) 

Professor Hesse link's research project 

- feature extraction using expert system 

(based on flow field topology rules) 

- display that Highlights critical matures and 

suppresses "visual clutter" 


The Scienti st’s Environment of the Future 

A sub stantial portion of the computing power will 
be distrib uted to the scientist's office 

Most scene rendering will occur in "real time" 

Most simulations will be performed interactively 
with die scientist observing and steering the 
simulation in "real time" with highly effective 
controls 


The Scientist 1 * Environment of the FutlUS 



Tn the scientists office 

Interface 

lane field of view stereo vision 
stereo sound 

6 degree^rf- freed om controls 
voice recognition 

Tailored fractions 
personality module 
Amtm management module 
visualization module 
special data analysis mndnle 
special data analysis module 
special data analysis module 


Conclusions 

Breakthroughs in workstations have greatly 
enhanced our capability to experience and 
understand computer simulations 

There is still room for major improvements 
in w orkstatio ns for scientific analysis. 

The most critical task facing scientists doing 
computer simulations is extraction and display of 
the key features of die simulations — more 
research is required on this task. 
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INTRODUCTION 


This report is the outcome of a project initiated to evaluation the use of SD/FAST, an 
equation generating program, in the specialized environment of the AD 100 simulation 
engine of the Visualization and Simulation Lab of The Research Institute. A program 
was written to simulate a slewing maneuver of a four body spacecraft model. The 
model is described on pages 13 to 36 of NASA JPL Technical Report 32-1592, "Attitude 
Dynamics Simulation Subroutines for Systems of Hinge-Connected Rigid Bodies" by 
G.E. Fleischer and P.W. Likins. This simulation was converted to ADSIM from the 
sample program, written in FORTRAN, found in Appendix B of Symbolic Dynamics' 
SD/FAST User's Manual. Plots were then generated in order to evaluate the 
performance and accuracy of the ADSIM simulation in this simple problem. 


IQQL 

SD/FAST is a program used for the generation of dynamic equations from a simple user 
input file describing the system to be modeled. This program can handle any system 
that can be described by a series of rigid bodies connected with hinges. The input file 
consists of descriptions of each body in the system (mass, inertia matrix, reference 
orientation), joint descriptions, and any constraints added to the system. From this file, 
the SD/FAST program generates two other files. The first is an initialization file that is 
used to initialize state variables and set up the environment in which the simulation is to 
run. The other file generated by SD/FAST is the dynamics file. This file contains all of 
the equations of motion needed to model the system. These two files generated by 
SD/FAST do not constitute a complete simulation themselves. The user must add all of 
the necessary controllers and sensors in his own file. The user must also add any data 
concerning prescribed motion in any part of the system to be modeled. This information 
is then used to drive the SD/FAST generated portions of the simulation (particularly the 
dynamics file). 
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CONSIDERATIONS 


There were many considerations that had to be addressed when converting the sample 
program from FORTRAN to ADSIM. The first thing that had to be done to realize that 
the two languages use SD/FAST differently. ADSIM uses the initialization file generated 
by SD/FAST as a place in which to incorporate its user-written code. The dynamics file 
generated by SD/FAST is incorporated into this framework. FORTRAN handles the task 
quite differently. The user includes calls to the initialization and dynamics files. The 
incorporation of the user's ADSIM code into the SD/FAST generated files makes the 
ADSIM much smoother and easier to follow. The addition of a sophisticated, variable 
time-step integrator within ADSIM was another major difference between the two 
languages. This integrator had to be written into the FORTRAN code of the sample 
program, but was not needed in ADSIM version. Other considerations that had to be 
addressed were the differences in the two program languages mechanics. The 
renaming of variables within the different subroutines of the FORTRAN example had to 
be traced back through the code in order to decipher which ones needed to be ported to 
ADSIM and which ones were only needed by the FORTRAN. The global scope of 
ADSIM variables dismissed the need for renaming these variables within different parts 
of the ADSIM code. Also, the logic that is used in the controller model of the user- 
written portion of the FORTRAN simulation had to be incorporated into different areas of 
ADSIM code since ADSIM deals with procedural statements differently. The 
mathematics of the controller model also had to be deciphered in order to better 
understand the dynamics involved with the slew maneuver as commanded by the user- 
written code. The ADSIM code that was obtained by this conversion still retains most of 
the flow of the FORTRAN code and should be easy to follow. 


MODEL 

As seen in the accompanying figure the spacecraft model contains a number of 
elements (four bodies with ten degrees of freedom). The spacecraft model used during 
the simulation includes a bus fixed to a boom by means of a u-joint. Atop the bus is a 
camera attached by means of shaft (or clock). The shaft is connected to the bus by 
means of a pin joint while the camera is connected to the shaft by another pin joint. The 
positions of these joints are Q-4 and Q-5, and their rates U-4 and U-5 respectively. 
These are the joints that have commanded rates during the simulation. The boom is 
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flexible; and to simulate this, spring and damping torques are applied to the U-joint 
between it and the bus. The craft is also equipped with three gas jet thrusters and rate 
sensing gyros to correct bus attitude should it exceed certain limits. The clock and 
camera can be slewed by means of variable torque motors at each hinge. A controller 
model in the program drives these motors. 

The reference configuration of the clock and the camera is pointing away from the boom 
with no inclination. At the beginning of the simulation, Q-4, the clock position, is 4 
radians from the reference configuration; and Q-5, the camera position, is given to be - 
0.5 radians. One second into the simulation, the slew will begin and will last for ten 
seconds. During this time, the clock will cover a 0.25 radian rotation from 4 to 3.75 
radians and the camera will cover a 0.10 radian decrease in the declination of the 
camera from -0.5 to -0.4. The plots of this can be seen in Appendix B. 

As the bus drifts, the controller fires the thrusters as necessary in order to correct its 
attitude. There is a deadband of .0025 radians before this occurs. The effects of these 
thrusters on the bus and the boom can be seen in the plots in Appendix B. 


PROGRAM 

Appendix A contains all of the simulation code and all input files associated with the 
simulation. The first item is the SD/FAST input file that was provided by Symbolic 
Dynamics. The next two items are the initialization file and dynamics file. Finally, there 
are several plots that were obtained from the simulation running on the AD 100. 

The simulation program is written in ADSIM with much of the initialization occurring in 
the FORTRAN sections at the beginning of the code. The ADSIM consists of the main 
user-written body of code and five functions that control each of the three thrusters and 
the two hinges that follow prescribed motion. The dynamics file is then incorporated into 
the code with an INCLUDE statement after the main body of code. The blocks are 
defined below. 
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The Example Spacecraft 
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Initialization File: 


FORTRAN blocks - These blocks set up the variables in FORTRAN. These will 

be reassigned as ADSIM variables later in the dynamics file. 
There is also a description of the system from the input file. 

REGION initial - This block initializes the stat variables and all initial values of 

all other variables used to drive the simulation. 


FUNCTION thrust 


- These three blocks each control the three thrusters that 
correct the attitude shifts of the bus. 


FUNCTION camrat 


- These two function control the commanded rate of the clock 
and camera as the slew maneuver is simulated. 


DYNAMIC continuous - This block of code contains both the equations for computing 

the hinge torques and the equations for computing the force 
applied to the system by each of the three thrusters. 

Dynamics File: 

REGION sd-initial - This section converts the FORTRAN variables in the 

initialization file into ADSIM variables. 


REGION sd-dynamics 


This block computes the state derivatives. 


REGION sd-states 


- This block performs the integration routine. 


NOTES ON THE SIMULATION 

The simulation environment consisted of a VAX 3100, 3400 and an Applied Dynamics 
AD-100. The AD-100 features the ADSIM language. The version of ADSIM used is 
version 7. This version has the ability to import FORTRAN code generated outside of 
the AD-100 environment. In this instance the VAX 3100 hosted the SD/FAST software. 
The spacecraft problem was input to the SD/FAST program in the 3100 and it took 
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27.06 sec of CPU time to create the FORTRAN coded equations of motion. The ADSIM 
portion of the code was input to the 3100 and then the ADSIM compiled. This 
compilation took 50.1 1 seconds of CPU time. The problem was then executed on the 
AD-100. One run (the one shown in the graphics) was made at 10 times real time. It 
was then found by experiment that the problem could be run in 96 times real time i.e. it 
took 0.3125 (30/96) seconds for a total problem time that would take 30 seconds in real 
time. The default integration scheme called the Adams-Bashford 2 was used in all 
instances. This is a single pass integration scheme suitable for real time or hardware in 
the loop simulations. It is one of over a dozen integration schemes selectable in the 
ADSIM environment (you may, of course, write your own if you wish). 

CONCLUSION 

After becoming familiar with the mechanics of SD/FAST in regards to ADSIM and 
converting the sample program from FORTRAN to ADSIM, the worth of a program such 
as SD/FAST becomes apparent. It gives the programmer more time to concentrate on 
the design and control of a system without having to worry about producing the long, 
complicated equations of motion of the system by hand. The time spent learning how to 
use SD/FAST is more than made up by the savings in man hours and computer time 
that it will save later. The sample situation, although far from complex, gave a solid 
performance when running under ADSIM. The relative strengths of the AD 100 and the 
ADSIM programming language are more than capable of running this program's 
generated files with almost no loss in performance except in the largest of simulations 
where runtime performance is of optimum importance. In these cases, SD/FAST is an 
option, although it may not yield as good results as other methods of deriving equations 
of motion. 
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Input File 
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# 

# This is an SD/FAST Input File describing the four-body spacecraft 

# from Appendix B of the SD/FAST User's Manual- Note that no external 

# forces are allowed and that only the base body can have any external 

# torques. 

# 

language = adsim 
body = bus 
mass = 410 

inertia = 115 -14 14 

^ Jib ~ 3 4 * £ 

Y 440 

force = 0 0 0 

body - clock i n b = pus i o i n x - om 
mass = 6.8 

i ner 1 1 a = 0 . 35 0 . 35 .> 

b od v T o J O' l n t = 0 0 - 0 • 7 5 

inbTo Joint = 0.0 -1.5 

pin = 0 0 1 

force = 0 0 0 
torque = 0 0 0 

body = camera inb 
mass = 57.5 
inertia = 4.85 


bodyTo Joint = 0 
inbTo Joint = 0 
pin = -1 00 
force = 0 0 0 

torque = 0 0 0 

body = boom mb - bus loint - u -joint 
mass = 10.7 

inertia = 27.2 0.2 27.2 

body To Joint =0 3.3 0 

inbTo Joint - 0 -1.2 0 

pin = 0 0 1 

pin = 100 

force = 0 0 0 

torque = 0 0 0 

# @(#)slew.dat 1.1 (Symbolic Dynamics, Inc. .) 88/05/02 20*29:46 


= clock joint = pin 

0.41 - 0.07 

2.2 0.54 

'? ^ 5.5 

- 0.1 - 0.75 
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4 

4 This is an SD/FAST Input File describinq the four -body spacecraft 
4 from Appendix B of the SD/FAST User's Manual. Note that no external 
4 forces are allowed and that only the base body can have any external 
if tor dues. 

4 

language = adsim 
body = bus 
mass = 410 

inertia = 115 -14 14 

? 31b -34.6 

> ? 440 

force = 0 0 0 

bodv = clock mb - bus ioint - om 
mass =6.8 

l ne r 1 1 a = 0 . 35 U . 35 0 

bodvToJomt = 0 0 -0.75 

inbTc'Juint 
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ORIGINAL PAGE iL 

OF POOR QUALITY 



Initialization File 
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n n n 


-ORTRAN SECTION declare 
-ORTRAN 


3D/FAST INITIALIZATION 


Z INITIALIZATION ROUTINE 

C GENERATED 16 -AUG -1990 12:38:31.66 BY SD/FAST, ORDER ( N :» FORMULATION 
(SDFAST A. 2. 14 #20112.) ON MACHINE ID 0A000005VS3100 
0 COPYRIGHT (C) 1988 BY SYMBOLIC DYNAMICS, INC. 

C DESCRIPTION OF SYSTEM: 



BODIES 

INB 


COOPDS 1 

:\ 

Q 

# 

NAME 

BODY 

HINGE TYPE 

P =PRESCRIBED 

i w 

1 

bus 


ATTITUDE 

1 2 

■—> 
■ j 

- 


doc k 

1 

TRANSLATION 
P I N <1 D > 

3 9 

4 

10 

c 

n 

earner a 


PIN <. ID) 

5 



4 

Doom 

1 

U -JOINT <’2D) 

6 7 



C THIS ROUTINE SHOULD BE EDITED TO SUPPLY VALUES FOR ALL VARIABLES 
C WHICH WERE GIVEN AS LiUESI’lON MARKS IN THE INPUT. THESE VARIABLES 
C ARE SET TO BEl-OW AND SHOULD BE CHANGED AS APPROPRIATE. 

C 

C MEANING OF VARIABLES: 

C 

C IN THE INDICES BELOW, V IS A BODY NUMBER, I AND J SELECT 

C ELEMENTS OF A VECTOR OR 3X3 MATRIX, AND Q SELECTS A PARTICULAR 

C PIN OR SLIDER AT A HINGE COUNTING UP FROM INBOARD TO OUTBOARD. 

C 

C ME no -- MASS OF BODY K 

C I ('. E , I , J ) -- IJ-TH ELEMENT OF INERTIA MATRIX FOR BODY E 

C LCQ.I) -- I-TH ELEMENT OF PIN VECTOR AT COORDINATE 0. SINCE 

- Q<4> IS THE 1ST HINGE COORDINATE. L<4,I? SELECTS THE 

C 1ST PIN OR SLIDER. 

C RE < E , I > -- I-TH ELEMENT OF VECTOR FROM MASS CENTER OF BODY E TO 

C ITS INBOARD HINGE FOINT 

C R I ( E , I :> -- I-TH ELEMENT OF VECTOR FROM MASS CENTER OF THE INBOARD 

C BODY OF BODY E TO THE INBOARD HINGE POINT OF BODY E 

C 

DOUBLE PRECISION ME (4) , I <4, 3, 3> 

DOUBLE PRECISION I. (7, 3) . RE <4, 3) 

DOUBLE PRECISION RI(4,3) 

COMMON BLOCE FOR SYSTEM PARAMETERS 

COMMON /SDCQMN/ ME, I , L . RE , R I 
C 

END FORTRAN 



FORTRAN SECTION start 
FORTRAN 
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c 

C VALUES FROM INPUT FILE 
C 

C BODY 1 >bus) 

C 


MK ( 1 

.) = 

410 

. ODO 

Id. 

1. 

1 

') = 

115 

. ODO 

Id. 

1 , 


) = 

-14 

. ODO 

Kl. 

1 9 


) = 

14. 

ODO 

Id. 

9 

1 

;> = 

-14 

- ODO 

Kl. 

.3 f 

*2 

;> = 

31b 

. ODO 

1(1, 

jil v 

wJ 

) = 

-34 

. bDO 

1(1. 

,wm t 
vj 9 

1 

> = 

14. 

ODO 


I ( 1 , 3, 2) = -34. EDO 

I Cl, 3 , 3 )= 440.000 
C 

C BODY 2 (clock) 

C 


MK(2) = fa 

. SDO 

I (2, 1, 1) 

= 0.35DO 

1(2, 1,2) 

= O.ODO 

I (2, 1 .3) 

= O.ODO 

1(2,2, 1) 

= O.ODO 

I (2, 2, 2) 

= 0.35D0 

1(2, 2, 3) 

= O.ODO 

I (2,3. 1) 

= O.ODO 

I (2,3,2) 

= O.ODO 

I (2,3,3) 

= O.ODO 

RKC2, 1) = 

0. ODO 

RK (2, 2) = 

0 . ODO 

RK ( 2, 3) = 

-0. 75D0 

RI (2. 1) = 

0, ODO 

RI(2.2) = 

0. ODO 

RI (2.3) = 

- 1 . SDO 

L (4. 1 ) = ' 

0. ODO 

L ( 4 , 2 ) = 

0 . ODO 

L (4,3) = 

1 . ODO 


c 

C BODY 3 (camera) 

C 

MK ( 3 ) = 57 . SDO 
I ( 3 . 1 , 1 ) = 4 . 85 D 0 
1 ( 3 , 1 , 2 )= 0.41 DO 
1 ( 3 , 1 , 3 )= - 0 . 07 DO 
I ( 3 , 2 , 1 ) = 0 . 41 DO 
I ( 3 , 2 , 2 ) = 2 . 2 DO 
I ( 3 , 2 , 3 )= 0 . 54 D 0 
I ( 3 , 3 , 1 )= - 0 . 07 D 0 
I ( 3 , 3 , 2 )= 0 . 54 D 0 
I ( 3 , 3 , 3 )= 5 . 5 DO 
RK( 3 , 1 )= 0 . 0 D 0 
RK ( 3 , 2 ) ■ - 0 . 22 D 0 
RK( 3 . 3 ) = 0 . 2 D 0 
RIO, 1 ) = 0 . 0 D 0 
RIO. 2 )- - 0 . IDO 
R I ( 3 , 3 ) = “(i * / SDo 
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L >. 5 , 1 :> = -l.ODO 
L <! 5 , 2 ) = O.ODO 
L C 3 , 3 ) = O.ODO 

c 

C BODY 4 boom ) 

C 

MK ( 4 > = 10.7D0 


I (4. 1, 

1 ) 

= 2 7 . 2D< 

I (4, 1, 

2 ) 

- O.ODO 

I (4. 1, 

3) 

= O.ODO 

I <4,2. 

1) 

= O.ODO 

I 14, 2 , 

2 :> 

= 0.2D0 

I I 4,2. 


= 0.0 DU 

I <4,3, 

1 ) 

= O.ODO 

I ‘4,3. 

2) 

— 0 • o D 0 

I ‘ 4 . 3 , 

3 ) 

= 27. 2D' 

R K < 4 , 1 

) - 

U.ODO 

R K i4, ; 

) = 

3.3 DO 

Rfc. < 4 , 3 

) s 

0. ODO 

R I C 4 . 1 

) = 

U.ODO 

R I < 4 . 2 

) = 

- 1 . 2D0 

R I C 4 , 3 

) = 

0 . ODO 

L <6, 1 ) 

= 

0 . ODO 

L (. 6 , 2 :> 

= 

O.ODO 

L<6,3:> 

= 

1 . ODO 

l ( 7, i:> 

= 

1 . ODO 

L ( 7 , 2 ) 

= 

0. ODO 

L < 7 , 3 > 


O.ODO 


c 

END FORTRAN 
FORTRAN_SECT I ON * 

REGION initial 

1 SET THIS FLAG NON -ZERO WHENEVER YOU WANT f'O CONFUTE 
! ANGULAR MOMENTUM, KINETIC ENERGY. ETC. 

sdchek flaa = u.ODO 

1 INITIAL CONDITIONS FUR STALE VARIABLES 

□10 = O.OdO 
Q 20 = O.OdO 
q_ 30 = O.OdO 
G_4* = 4 . OdO 

Q_5@ =-0.5d0 
Q_E@ = O.OdO 
Q_7S = O.OdO 
0_3<3 = O.OdO 
□_9@ = O.OdO 
Q _ 1 0@ = O.OdO 
U _ 1 £ = O.OdO 
U_2Q = O.OdO 
U_3@ = O.OdO 
u 4^ = O.OdO 
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u 5<§ = O.OdO 
U_6<S = OdO 
U„73 = O.OdO 
u_se = O.OdO 
U _9<s! = O.OdO 
U_10@ = O.OdO 

! EXTERNAL FORCE AND TORQUE INITIALIZATIONS 


FK 

_ 1 . 

_i = 

0 . 0D0 

FK 

_1 

— 

0 , ODO 

FK 

_1 

0 — 

0. ODO 

TK 

_1 

” 1 = 

0 • OdO 

TK 

.1 

2 = 

0. OdO 

TK 

..1 

_3 = 

O.OdO 

FK 

2 

.1 = 

( i . ( ') D 0 

FK 


Z — 

0 . ODO 

FK 


<— \ _ 
■+J — 

0 . ODO 

TK 


_ 1 = 

0 . ODO 

TK 


•L — 

0 . ODO 

TK 

-2 

w ~ 

0 . ODO 

FK 

_ y-l 

_ 1 = 

0 . ODO 

FK 

— 


0 . ODO 

FK 


~3 = 

O.ODO 

TK 



_1 = 

O.ODO 

TK 

o 

o 

2 = 

0 . ODO 

TK 

3.3 

J3 = 

O.ODO 

FK 

_4 

_1 = 

0. ODO 

FK 

_4 

■7 - 

0 . ODO 

FK 

_4 

^ — 
_ *J — 

0. ODO 

TK 

_4 

" i - 

0 . ODO 

TK 

_4 

-* _ 

0. ODO 

TK 

4 

1 _ 
0 — 

0 . ODO 


! HINGE FORCE OR TORQUE INITIALIZATIONS 

TAU 4 = O.OdO 
TALI 5 = O.OdO 
TAU S = O.OdO 
TAU 7 - = O.OdO 

! a-frame initialization 


WA__ 

1 

= 

0. OdO 

WA _ 

2 

= 

0 . OdO 

WA _ 

■_> 

= 

O. OdO 

VA^ 

1 

= 

O.OdO 

VA_ 

A 

3 

0. OdO 

VA 

•— 1 

= 

O.OdO 


! USER ADDED VARIABLE ASSIGNMENTS 
CLKCMD = 4.0D0 
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CAMCMD = -0.5D0 
TIME = o. UdU 
CLKCMD = 4.0D0 
CAMCMD = -0.5D0 
qyro = 2. OdO 
1_1 = 0. 23d 0 

I ^2 = 0 . 21 dO 
1_3 = 0.31d0 
k_l = 3500. OdO 
k_2 = 3500. OdO 
k_3 = 2000. OdO 
k_4 = 2000. OdO 
b_l = 20. OdO 
b_2 = 20. OdO 
b_3 = 10. OdO 
b_4 = 10. OdO 
END REGION initial 

function thrustl(trl = time, err _ 1 ) 

dband = 0.0025d0 
tmin = 0.02d0 
fire_l = O.OdO 
toft ~ 1 = O.OdO 

if ( er r _ 1 . 1 1 . -d band ) ■chan 
f i r e _ 1 = 1 . Od 0 

toff 1 * time + tmin 

end if 

if (err 1 . at. dband) then 
fire_l = -l.OdO 

toff_l = time +■ tmin 

end i f 

if C time. ge. toff _1 > then 

fire_l = O.OdO 

end if 

trl = f ir e_ 1 
end function 

function thrust2(tr2 * time, er >-_2) 

dband = 0.0025d0 
tmin = 0.o2d0 
fire 2 = O.OdO 
toff_2 = O.OdO 

if (err 2. It. -dband) then 
fire. 2 = l.OdO 

toff_2 = time + tmin 

end i f 

if ' err 2. at . dband ) then 


D-17 



fire_2 = -l.OdO 

toff_2 = time + tmin 

end 1 f 

if ( time. ae. toff ) then 
fire 2 = O.OdO 

end 1 f 

t r 2 - f ire 2 
end function 

■function thrust3(tr3 = time, err 3) 

dband = 0.0025d0 
tmin = 0.02d0 
fire_3 = O.OdO 
toff _3 = O.OdO 

if (er r_3. 1 t . -dband) then 
fire_3 = 1 . OdO 

t of f _3 = time *+• tmin 

end i f 

if (er r_3. gt . dband > then 
fire__3 = -l.OdO 

toff_3 = time + tmin 

end i f 

if ( time. ae. toff .3 ) then 

fire 3 = O.OdO 

end if 

t r 3 = fire 3 
end function 

function cmdrat 1 ( c 1 k cmd = time.) 

STRTSL = I. OdO 
STOPSL = 11. OdO 
elk rat = -0.025d0 
if C t ime. 1 t . strtsl > then 
c 1 k cmd = 4. OdO 
end i f 

if ( ( ti me. qt . strtsl). and. (time. It. stops 1 ) ) then 
elk cmd = 4. OdO + (elkrat * ( time- 1 . OdO) ) 
end if 

if ( t ime. qt . stops 1 ) then 
elk cmd = 3.75d0 
end i f 

end function 

function cmd r at2 ( cam cmd - time.) 
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STRT5L = l.OdO 
STOPSL = ll.OdO 
camrat = O.OidO 
if ( time. 1 t . strtsl ) then 
camcfTid = -O.SdO 
end i f 

if < (time. gt. strtsl). and. c time. 1 1 . stoosl 
camcmd = -O.SdO + f . camrat * (time-i. 
end i f 

if C t ime. at ■ stoosl !> then 
camcmd = -0.4d0 
end if 

end tunc t ion 


DYNAM I C continuous 

clkcmd = cmdratl (time.) 

camcmd = cmdrat2< time) 

time = system^time 

tau_4 = -k_l*(q_4 - CLKCMD) 

tau_5 = -k_2*<qJ5 - CAMCMD) 

tau_S = ~k_3*q_6 

tau_7 = -k_4*q _7 

tri = thrust 1 < time, err_l > 
err _ 1 = qyro * u_l + q_ 1 
t k _ 1 _ 1 * 1_1 * tri 

tr2 = t h rust 2<! time, err j2> 
err_2 = qyro * u 2 + a 

tkj,2 = 1..2 * t r 2 

t r 3 = thrust 3c! time* err 3) 
err 3 - ayro * u 3 ■+* a 

tk.i.3 =1,3 * tr3 

enddynamic 

INCLUDE f SLEW . DYN f 

r unspecs 

endtime = 30, 
speedup = 10 

EXECUTE ' SDCONTINUGUS' 


- b_l*u_4 

- b_2*u 

- b_3*u_S 

- b 4*u 7 


) ) then 
OdO) ) 
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Dynamics File 
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REGION sd initial 


' VARIABLES FROM l'HE SD/FASi INPUT FILE 


MK 

1 


= C M K (1)3 




I .. 

1 . 

i . 

_1 

= C I ( 1 , 

1, 

1 : 

> 3 

I _ 

1 _ 

i 

2 

= C I ( 1 , 

1 

2 

) 3 

I_ 

1_ 

i_ 

'-y 
_ >-j 

= LKl. 

1 , 

0 . 

) 3 

I _ 

1_ 

2 

_1 

= C I ( 1 , 

42 

1 

) 3 

I_ 

1_ 

2 

2 

= C I ( 1 , 

2 1 

4 - , 

) 3 

I 

1_ 

2 

“3 

_ 

= [1(1, 

2 

■—v 

O 

) 3 

I_ 

1. 


^1 

= ci a. 


1 

) 3 

I_ 

1_ 

’ 3 ' 


= C I ( 1 . 

O 

2 

) 3 

I_ 

1_ 


_ N_> 

= c i c 1 , 


wl 

) 3 

MK 

2. 

= 

c m k: c 2 ) 3 




I_ 


.1. 

_1 

= [1(2, 

1 

1 

) 3 

I _ 


.1 

2 

= Cl (2, 

1 

2 

) 3 

I_ 


,1. 

__ wl 

= Cl (2. 

1 

w 

) 3 

I . 

4 ^ 

4 ^ 

_1 

= CI(2, 

2 

1 

) J 

I_ 

2 

2 


= CKO. 

2 

2 

) 3 

I 

2 

2 


= Cl (2. 

2 

*-y 

O 

) 3 

I _ 


'*3 

_1 

= Cl (2. 

2 

1 

) 3 

I . 

*L. _ 

••J 

2 

= Cl (2. 

0 

2 

) 3 

I_ 

*T« 


O 

= Cl <2. 

0 

0 

) 3 

RK 

2 

1 

= CRKC2. 

1 

J ] 


RK 


42 

= CRKC2. 

2 

.) J 


RK 

ju! 

3 

= CRKC2, 

yJj 

> ] 


RI 

2 

1 

= C R I ( 2 , 

1 

) 1 


RI 

2 

2 

= C R I ( 2 , 

2 

:> j 


RI 

2 


= CRIC2. 

0 

:> -j 


L 

4 

_1 

= 

CL (4, 1) 

3 



L_ 

4 

■ ' 

= 

CL (4,2) 

3 



L. 

4 

O 

= 

C L C 4 . 3 ) 

3 



MK 

. O 

= 

CMKC3) 3 




I_ 

O 

_1 

_1 

= C I C 3 , 

1 

• i 

) 3 

I_ 

3^ 

’l 

-fc_ 

= CIO, 

1 

* *- 

) 3 

I _ 

<— 1 

• ^ - 

.1 


= CIO, 

1 

•— 1 
f 0 

) 3 

I_ 

O . 

~ 1 


= CIO. 

2 

, 1 

) 3 

I 


2 

2 

= CIO, 

-> 

v 2 

) 3 

I_ 



O 



= CIO, 

2 

V O 

) 3 

I, 


■“1 

1 

= CIO, 

■y 

O 

1 

) 3 

I_ 

*—> 

'—I 

_ O 

_2 

= CIO, 

O 

v_; 

f <*- 

) 3 

I. 

_3 

1 3 

: 3 

= CIO, 

3 

1 ^ 

1 f ^ 

!) 3 


RK 3 1 = CRKO, 1)3 
RK_3 2 = CRKO, 2) 3 
RK 3 3 = CRKO, 3) 3 
R I 3 1 = CRIO, 1 ) J 
R I 3_2 = CRU3.2;] 

R I ” 3 3 = CRIC3.313 
L_5 _ 1 = CL (5.1)3 
L 5 2 = CLC5.2) J 
L 5 3 = C L ( 5 , 3 ) 3 
MK 4 = C M K ( 4 3 3 

I 4 _ 1 . 1 = C I (4. 1 , 1 ) j 

I ",4„ 1 .2 = Cl (4, 1.2) 3 
I ~4ll._3 = Cl (4. 1,3) J 
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14 2 1 

= CI <4, 2, 1)1 

I _4_2 .2 

= C I C 4 , 2 . 2 ) ] 

I 4_2_3 

= Cl (4,2,3) ] 

I 4 3 _ 1 

= : i <: 4 , 3 . i :> i 

14 3 2 

= C I (4,3. 2)] 

14 3 3 

= C I < 4 , 3 , 3 ) ] 

R K. 4 1 

= CRKC4, 1.) ] 

RK 4 2 

= L R K ( 4 . 2 ) I 

RK 4 3 

= C R K ( 4 . 3 ) j 

RI 4 1 

= C R I ( 4 , 1 ) ] 

RI 4 2 

= C R I ( 4 . 2 ) ] 

PI 4 3 

= C R I i. 4 . 3 ’) 1 

1-6 1- 

C L <: 6 . 1 ) J 

L .6 .2 = 

C L < 6 , 2 ) 1 

L _fa _3 - 

CLCS.3) I 

L_7 _ 1 = 

C L ( 7 , 1 .) 1 

L_7_2 = 

CLC7.2) ] 

L_7 ~3 = 

L L ( 7 , 3 ) ] 

» REGION =d 

_ initial 

ION 3d_dynamic5 


SUBROUTINE SDNSIMCQ, U, FK, TAU, TK, WA. VA, QDOT, UDOT) 


1 

I 

I 


THIS IS THE DERIVATIVE SECTION FOR A 4 -BODY 
SYSTEM WITH 4 HINGE DEGREE'S) OF FREEDOM, 
PLUS 3 ATTITUDE AND 3 TRANSLATIONAL DEGREES 

GENERATED 16 -AUG -1990 12:38:30.91 BY bD/FAST 
(SDF.AST A. 2. 14 #20112) ON MACHINE ID 0 AO 0000 


COPYRIGHT i.C) 1988 BY 1 

DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 

DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 


YMBGLIC DYNAMICS, INC 


0 

<; 1 0 ) 

u 

C 1 0 > 

FK 

(4, 3) 

TAU 

(7) 

TK 

(4,3) 

WA 

(3) 

VA 

(3) 

QDOT 

< 1 0 ) 

UDOT 

( 1 0 ) 

RK 

(4,3) 

RI 

(4,3) 

L 

(7,3) 

I 

(4,3,3) 

MK 

c 4 :> 

I AD J 

( S o 'i 

Cl K. 

<•5,3,3 :> 

WADOTC 

f- 3 .) 

RSK 

< 5 , 3 .> 

MRSK 

<5.3 !) 

K I K 

■- 3 . '-j ) 


FREE-FLYING 
OF FREEDOM. 


, ORDER (N> FORMULATION 
5VS3100 
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DOUBLE 

PRECISION 

WPK 

(7,5,3) 

DOUBLE 

PRECISION 

VPK 

(7,5,3) 

DOUBLE 

PRECISION 

WK 

(5,3) 

DOUBLE 

PRECISION 

WI 

(5,3 

DOUBLE 

PRECISION 

TTK 

(5,3) 

DOUBLE 

PRECISION 

ftk: 

(5,3 ) 

DOUBLE 

PRECISION 

OK 

( 5 , 3 > 

DOUBLE 

PRECISION 

AK 

, o ) 

DOUBLE 

PRECISION 

AKI 

( J , O ) 

DOUBLE 

PRECISION 

GMK1 

( *2 , a. , O , O ) 

DOUBLE 

PRECISION 

OGMK 1 

• -L- , •*— , O , O ' 

DOUBLE 

PRECISION 

GFk 1 

(2,3 ) 

DOUBLE 

PRECISION 

OGFK 1 

(2.3) 

DOUBLE 

PRECISION 

ACC1 

( 3 ) 

DOUBLE 

PRECISION 

ALPH1 

( 3 .? 

DOUBLE 

PRECISION 

r< I f . 1 

( _j , ) 

DOUBLE 

PREC I S I ON 

WKH I 1 

( 3 ) 

DOUBLE 

PRECISION 

N R S K 1 

{ 3 .> 

DOUBLE 

PRECISION 

GflK 2 

' • — f ^ ’ 

DOUBLE 

PRECISION 

QCMK 2 

( 2 , 2 , 3 , o 

DOUBLE 

PRECISION 

G F K — 

( jL. , w ) 

DOUBLE 

PRECISION 

OGFK 2 

(2,3 ) 

DOUBLE 

PRECISION 

ACC2 

(3) 

DOUBLE 

PRECISION 

ALPH2 

(3) 

DOUBLE 

PRECISION 

RIKT2 

(3,3) 

DOUBLE 

PRECISION 

WKRI2 

<3) 

DOUBLE 

PRECISION 

WRSK2 

(3) 

DOUBLE 

PRECISION 

GMK3 

jL. , f S-J } 

DOUBLE 

PRECISION 

OGMK 3 

y ^ , o.' 

DOUBLE 

PRECISION 

GFK 3 

(2,3) 

DOUBLE 

PRECISION 

OGFK 3 

(2, 3) 

DOUBLE 

PRECISION 

ACC3 

(3) 

DOUBLE 

PRECISION 

ALPH3 

(3) 

DOUBLE 

PRECISION 

RIKT3 

(3.3) 

DOUBLE 

PRECISION 

WKRI3 

( 3 ) 

DOUBLE 

PRECISION 

WPSK3 

(3) 

DOUBLE 

PRECISION 

GNK4 

, W « 

DOUBLE 

PRECISION 

OGMK 4 

(2,2, 3.3) 

DOUBLE 

PRECISION 

GFK. 4 

(2,3) 

DOUBLE 

PRECISION 

OGFK 4 

( 2 , 3 > 

DOUBLE 

PRECISION 

ACC4 

(3) 

DOUBLE 

PRECISION 

ALPH4 

(3) 

DOUBLE 

PRECISION 

RIKT4 

(3,3 ) 

DOUBLE 

PRECISION 

WKRI4 

(3) 

DOUBLE 

PRECISION 

WRSK4 

(3) 

DOUBLE 

PRECISION 

GMK5 

(2, 2, 3, 3) 

DOUBLE 

PRECISION 

0GNK5 

( 2 , 2 , 3 , 3) 

DOUBLE 

PRECISION 

GFK 5 

(2,3 ) 

DOUBLE 

PRECISION 

0GFK5 


DOUBLE 

PRECISION 

ACCS 

(3) 

DOUBLE 

PRECISION 

ALPH5 

(3) 

DOUBLE 

PRECISION 

RIK.T5 

(3,3) 

DOUBLE 

PRECISION 

WKRI5 

(3) 

DOUBLE 

PRECISION 

WRSK5 

(3) 

DOUBLE 

PRECISION 

LO 

r •"* -I* r> "D \ 

, U y 

DOUBLE 

PRECISION 

D I 

f ” ”• T* 'T* \ 

■ ♦ V-/, W 
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DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 


V I M 1 

TAUP1 

MPPl 

(2.2) 

VIP1 

(2,3 ) 

V I M2 

TAUP2 

MPP2 

C 2.3) 

V I P2 

C 2 , 3 ) 

VI M3 

TAUP3 

MPP3 

(2,3) 

VIPS 

(. 2 « 3 .) 

VIM4 

TAUP4 

MPP4 

(2, 3) 

VIP4 

(.2.3) 

VIMS 

TAUP5 

npP5 

< 2 . 3 ) 

v [ PS 

(2.3 ) 

V [Mb 

TAUPb 

MPP6 

( 2 . 3 ) 

VI Pb 

(.2,3) 

V I M 7 

TAUP7 

MPP7 

(2, 3) 

VIP7 

(2,3) 

SOLNb 

(6) 

WE 

(6) 

V6 

(6) 

LU 

(3,3) 

L21 

(3,3) 

L22 



w , w J 

D 1 

(3,3 '.) 

D2 

{ <3 r> 'i 
1 , O J 

D1INV 

(3,3 ) 

D2INV 

(3,3) 

L11D1 


D1L21 

( ry \ 

L22D2 

CO \ 

y w , o ) 

CL 1 1 

(. 3 , o .) 

CL 22 

(3,3) 

RCL 

( 3.3) 

CL 1 1 D I 

(3,3) 

CL22D2 

(.3,2) 

RCLD1 


CF 

TEMP < 100.) 
SI , Cl 
32, C2 
3 , C O 

54 , C4 

55, C5 
Sb, CE 
S7 , C7 

i O 1 

■ ^ j J 
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COMPUTE SINES AND COSINES jP j 

5 1 . C 1 = sin.cos(Q.l) 

52. C2= sin_cos(Q 2) 

53. C3= sin .cos (0 3) 

34 . C4 = sin_cos(Q _4) 

35, C5= sin cos'.'Q 5) 

56 , Cb = sin .cosdl.b) 

57, C7= sin_cos(Q_7) 

COMPUTE DIRECTION COSINES 

( C2*C3 ) 

-( 33*C2) 

c csi*< S2*C3) ) +(S3*ci ) :> 

( (Cl *C3 ) -'SI* ( S2*S3 ) ) ) 

-<si*c 2 ) 

< (SI *S3 :> - <C 1 * < S2*C3 ) ) > 

( (SI *C3 ) + ( C 1 ■* ( 32*33 J ) > 

<ci*C2) 

COMPUTE QDQT CKINEMATICAL SONS ) 

WADOTC 1= ((Cl K 1 3 1 *WA 3) + C (CIK 1 1 . 1*WA_1 ) +(CI K_ 1 _2_ 1*WA_ 

WADOTC _2» ( (ClK 1 3 2*WA 3 ) + ( ( C I K 1 1.2*WA 1)+(CIK 1 _2_2*WA_ 

WADOTC _3 = ( (CIK 1 J 3*WA 3) + ( <CIk 1..2_.3*WA.2) -KWA_1*S2) ) > 

QDOT_ 1 = ( ( (U 1 -WADOTC 1 ) *(03/02) > + ( (U .2 -WADOTC 2>*<-S3/C2) )> 

QDOT _2 * ( (S3*dJ 1 -WADOTC 1 .> ) + (C3* (U..2-WADQTC_2) ) ) 

QDOT~3* ( (U.. 3 -WADOTC 3)+(((U 1 -WADOTC 1 )*(-<. S2*C3) /C2> > + (( U _ 
WADOTC _2) *( (32*33) /C2> > > ) 

QD0T_4= U_4 
QD0T_5 = U_5 
QDOT_b = U_S 
QD0T_7 = U _7 

□DOT_8 = ( ( (U_10*S2> + ( ( C I K_ 1 _ 1 _ 1 *U _8 ) + ( C I K _ 1 _ 1 _2*U _9 ) ) ) + < ( (Q_ 
WA_3) - ( Q_ 1 0*WA _2 ) ) -VA_1 ) ) 

QD0T_9= ( ( (Cl K_ 1 _2_3*U ..10) + ( (Cl K _ 1 _2_ 1 *U _3) + (Cl K_1 _2_2*U _9) ) 
( (7 q_10*WA_1) - ( C3_8*WA_3 ) ) -VA_2) ) 

QD0T_ 1 0 = ( ( ( C I K_ l _3_3*U _ 1 0 ) + C ( C I K _ 1 _3 _ l *U _S )+‘CIK_l_3 _2*U _9 ) 
( ( (Q_3*WA_2) - ( Q_9*WA_ 1 ) ) -VA_3) ) 

COMPUTE WPK (PARTIAL ANGULAR VELOCITIES) 


CIK_1_1_1 = 
C I K_ 1 _ 1 _2 = 
C I K_ 1 _2_ 1 = 
C I K_ 1 _2_2 = 
CIK_1_2_3= 
CIK_1_3 1= 
C I K _ 1 _3_2 = 
CIK: 1 3 3 = 


COMPUTE WK (ANGULAR VELOCITIES) 

AND WI (U-JOINT PARTIAL ANGULAR VELOCITIES* 


wk..2_i = 

WK_2 _2 = 
WK_2_3= 
WK_3 _1 = 
WK „3_2 = 
WK..3 _3 = 

WK_4_ 1 = 
WK 4_2= 


( ( U 1 *C4 )■*■'. U 2 *S4 ) ) 

( (U_2*C4) U _ 1 *34) ) 

(U..3+U 4) 

(WK 2_1 -U 5) 

( (WK 2 2 *05 ) -( WK 2 3*S5) > 
<. t. WK_2._2#S5) + '■ M- 2 5) ) 

i (. U _ 1 ^CG J + UJ 2^'db) 7 
C < U _2*C6 ) - < U _ i *ofc ') ) 


2 ) > ) 
25 ) ) 


3* 

) + 

) ) + 
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WK_4_3= 

CU._3+U_6) 

WK .5 _1 = 

( U..7+WK..4 1 ) 

WK_5_2 = 

( CWK_4_2*C7) + < WK .4 _3*S7) ) 

WK_5_3= 

( ( WK..4. 3*07) -< WK 4 2*S7'> ) 

COMPUTE OK 

(ANGULAR ACCEL. < 

OK .2.1 = 

CU 4*WK 2 2) 

OK .2.2 = 

- ( U 4*WK 2 1 .» 

QK_3_2= 

( (OK. .2 2*C5 ) - ■ U 5*WK 3 .3) 

OK 3 3 = 

( C OK .2 _2*S5 ) + < U 5*WK 3 2) 

0K_4 1 = 

(U 6*WK 4 2) 

OK 4 2 = 

-CU.£*WK 4 1) 

0K„5_2= 

( (QK_4 _2*C7) + t IJ _7*WK_5 .3) 

OK 5 3 = 

-( (OK _4_2*37 > +(U _7*WK_5_2 


COMPUTE VPK (PARTIAL VELOCITIES) 


COMPUTE RSK/MRSK 


COMPUTE RIK AND IADJ 


COMPUTE AK C L I NEAR ACCELERATIONS) 

AK_1..1 = C(U ,2*U .10) -<U „3*U 9) ) 

A K _ 1 _ 2 = C CU 3*U_8)-CU. l*U 10)) 

AK 1 3 = C CU 1*U .9) -(. U 2*U B) .) 

AKI 2 1= (AK 1 1-C1.dD0*CU 1*U .3) ) ) 

AK I .2 .2 = ( AK .1 2 - 1 . SDO* C U 2*U 3 ) ) ) 

AK I .2 3 = ( AK . 1 3+ ( 1 . SDO* < <U..1*U 1) + ( U 2*U_2) ) » ) 

AK JZ 1 = C (AKI .2 _1*C4>-KAKI 2 _2*S4) > 

AK _2 2= ((.AKI .2 „2*C4) -CAKI_2 1*S4) ) 

AKI 3 1= (AK 2 1 -CO. IDO* (WK . 2 1*WK 2 2))) 

AKI ,3.2 s (AK 2 2+ ( 0 . 1 DO* C ( WK _2 . 1 *WK .2 1 ) + ( WK_2 3*WK _2 .3 
AKI '3.3= (AKI _2 3 - ( 0 . 1 DO* ( 0K_2 l + CWK 2_2*WK_.2_3) ) ) ) " 
AK_.3_2 = ( (AKI .3 _2*C5) - ( AK I _3_3*S5 ) ) 

AK~3 _3= ( (AKI_3 .2*S5) + (AKl”3_3*C5) ) 

AKI _4_. 1 = ( AK_ 1 _ 1 - ( 1 . 2D0* < U~ 1 *U_2) ) ) 

AKI _4^2= ( AK_l_2+( 1 . 2D0*( (U_1*U_1 ) +(U_3*U_3) ) ) ) 

AKI _4_3= ( AK~ 1 _3 - ( 1 . 2D0* ( U_2*U_3 ) ) ) 

AK_4 1 * ( <AKI_4_1*CS) + (AKI~4_2*S6) ) 

AK_4 ~2= ( (AKl”4_2*C£) -(AKI _4_1*SS) ) 

AK_5 _2= ( ( AK_4__2*C7 ) + ( AK I _4_3*S7 ) ) 

AK_5~3= ( ( AkT_4_3*C7) -(AK~4~2*S7) ) 

COMPUTE FTK AND TTK (INERTIA FORCES AND TORQUES) 

WRSK3 1= 

FTK 2... 1 = 

FTK 22= 

FTK 23= 

.) .) ) .) 

FTK 3 1 = 


- ( ( 0 . 2D0*WK 3 2 ) + ( 0 . 22D0*WK 3 3 ) ) 

(6. SDO* ( AK 2 1 + ( O. 75DO* ( OK 2 2+('WK 2 1*WK 2 3) 
(6, SDO* ( AK 2 2*( 0. 75D0*( (WK 2 2*WK’ 2 3) -OK _2 1 
(£.SDO*(AKI 2 3 - ( 0. 75DO* ( ( WK 2 1 *WK 2 1)+(WK 2 

( 57 . 5D0* (AK t 3 1 + ( < ( 0 . 2 2 DO* ( wK 3 l *WK 3 2 


.) ) ) ') 


) ) ) ) 

) ) .) ) 

. 2*WK_ 2 .2) 
. 2D0* ( 
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WK_3_1*WK...3_3> ) > -< <0. 2D0*0K_3_2!> + (0'. 22D0*0K_3_3) ) .) ;> ) 

FTK_3_2= < 57. 5D0* ( AK_3_2+ C < 0. 2D0*0K_2_ 1 ) + < ( WK_3_3*WRSK3_ 1 > - ( 

0 . 22D0* (' WK_3_ 1 *WK_3_ 1 ) !> !> ) ) ) 

FTK_3_3= (57. 5D0* C AK_3_3-t- ( < 0. 22D0*0K_2_1 > + ( (0. 2D0* ( WK_3_ 1*WK_3_ 1 > 

) -( WK._3_2*WRSK3_1 !> J) !> > 

FTK_5_ 1 * ( 10. 7D0* C AK_4_ 1 + <3. 3D0 * '■ 0K_5_3 - < ! WK_5_ l *WK_5_2!> ) ;> > ) 

FTK _5_2 = C 10. 7D0* ( AK_5_2+ ( 3 . 3D0* ( < WK._5_l*WK_5_l ) + < WK~5_3*WK_5_3> ) 

> ) ) 

FTK _5_3* <10.7DCMUAlt .5 ,3-( 3. 3D0*< OK.. .4 . l+(WK_5_2*WK .5_3!>:> !> ) !) 

TTK_ 1 _ 1 = ( ( (U_2* ( ( 440 . 0D0*U .3 ) + ( ( 1 4 . 0D0*U _ 1 ) - C 34 . 6D0*U_2 .) ) ') ) - 
C U 3* ( C < 3 16 . 0DO*U 2) -(14. 0D0*U 1 > > - ( 34. 6D0*U .3 ) .) > ) -TK. 1 _ 1 > 

TTK 1 2= ( C < L) 3+ < ( 14. oL0*U 5 > + • ■ 115. 0D0*U 1 ) -( 14. 0D0*U_2 ) ) ') ) - 
< U 1*( (440. 0DO*U 3 j + ( '■ 1 4 . UiJO+O 1 > - > 34. 6D0*U 2) ) .» !> ) -TK 1 2> 

TTK. 1 3= ( ( (U 1*( 1 (316. ODO*U 2 .) -< 14. oD0*U_ 1 ) ) -(34. 6D0*U _3 ) !• ) - 

(U 2* ( ( 14. ODO*U 3.) + * < 1 15. 0UO*U l < - < 1 4 . OD0*U 2) ) ) > ) -TK_1 ..3) 

TTK 2 1= ( <4. 17500* (OK 2 l - « WK 2 2 +WK 2 3) ) .) -(5. 1D0*AK_2..2) ) 

TTK_2 ,2= ( (4. 1 75DU4 < UK 2 2+< WK .. 1 *WK 2 3) ) .) + < 5. 1 D0*AK.._2_ 1 ) ') 

TTK '3.' 1 = (((11. 5Do*AK 3 2 > + ( 1 2 . 65D0*AK .3 3 ') ') + C ( ( WK _3„2* ( ( 9 . 2S3D0* 
WK _3_3.) + ( (3. o7D0*wk 3 2) - <U.07DO*WK 3 1 .1 .) .» .1 - (WK 3 3* C (3.07D0* 

WK 3..3- 1 + ( (0. 41D0*WK 3 1 > + (4. 5D0*WK. _3_2) ) ) ') ) + < ( <0. 4 1 D0*0K_3_2 ') + 
(9. 933D0*0K:„.2 _1 !>.»-» 0. o7D0*0K_3_3 > ' J 

TTK_3_2 = ( ( ( ( 3 . 07D0*0K J3 _ 3 > + ( ( 0 . 4 1 D0*GK_2_ 1 .* + ( 4 . 5D0*0K_3_2 ) ) ) +■ 

( ( WK_3_3*( < (0. 4 1 DO*Wk _3 _2) +(9. 333D0*WK_3_1 > ) -(0. 07D0*WK_3_3) ) ) - 
C WK_3_ 1 <9. 283D0*WK _3_3 ) + ( <3. 07D0*WK_3_2> -(0. 07D0*WK_3_1 > ) ) :> > ) 

-(11. 5D0*AK I _3_ 1 ) ) 

TTK_3_3 = ( ( < (9. 283D0*0K_3_3) + C ( 3. U7D0*0K_3_2) - ( 0 . 07D0*0K_2_ 1 ) ) ) + 
("CWK._3_1*( ( 3 . 07D0*WK_3_3 ) + ( CO . 4 1 D0*WK_3_I ) + (4. 5D0*WK_3~ 2) > ) ) - 
<WK:_3_2*C ( (0.41 D0*WK_3_2.) + < 9. 933D0*WK_3_1 ') ) -(0.07D0*WK_3_3) ) ) :> ) 
-( 12. 65D0*AKI_3_1 > ) 

TTK_5_1 = ( < < 143. 523D0*( WK_5_2*WK_5_3> ) + ( 143. 723D0*0K_4_1 ) ') -( 

35. 31D0*AK_5_3) :> 

TTK_5._3= ( (35. 31D0*AK_4_1 ) + ( ( 143. 723D0*0K_5_3> -( 143. 523D0* ( 

WK. 5_1*WK_5_2) :> ) ) 

COMPUTE GMK AND GFK (GENERALIZED MASS 2ND FORCE) 

T AUP7 = (0. 0 (!) 69579286 0 0 8 5 D 0 * ( TAU..7 -TTK .5 1 ) !> 

GFK 5 _ 1 _3 = (FTK 5 . 3 - ( 35. 3 1 D0*TAUP7 .) 

GFK5 _2 1= (TTK. 5. l+( 143. 723D0*TAUP7) ) 

0GMK4 1 . 1 2 .2* <2. 02500643599 146D0 -MS. b/4993564O0855D0* (C7*C7) .) 

.) 

0GMK4 _ 1 _ 1 _2_3 3 (8. 67499356400855DO* ( S7*C7 ) > 

GGMK4 _1 _1 _3_2* 0GMK4_ 1 _ 1 _2 3 

0GMK.4 ~l ”l _3_3» (2. 02500643599 1 46D0+ ( 3 . 67499356400655D0* ( S7*S7 ) ) 

;» 

0GMK4 _ 1 _2_ 1 _2 * -<35.31D0*S7> 

0GMK4 _ 1 _2 _ 1 _3 = (35. 31D0*C7 .> 

0GMK4~2~1 _2_ 1 = -(35.31 D0*S7 > 

0GMK4_2_ 1 _3_ 1 = (35. 3100*07!) 

0GMK4~2_2_2_2= (0.2D0+C 143. 523D0* ( S7*S7) ) ) 

0GMK4~2~2_2_3= -( 143. 523D0*(S7*C7) !> 

□GMK4~2_2_3_2 - 0GMK4 _2_2 _2_3 

0GMK4 2_2_3_3= (0.2D0+C 143. 523D0* CC7*C7!» !> ) 

CF_2= ( CFTK_5_2*C7:> -(GFK5_1 _3*S7) ) 

CF_3= ( ( FTK_5_2*S7!) + ( GFK5... 1 _3*C7 ) ) 

QGFK4_ 1 _ 1 = FTK_5_ 1 
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0GFK4 .1 2= Cf 2 
GGFK.4 _1 _2 = CF 
GGFK4_2_1 = GFK5_2_1 

GGFK4_2_2= ( CO. 2D0*CQK. 5 .2*C7'i .. -< TTK .5 ,3*S7'i » 

OGFK 4 _2_3 = < CO. 2D0* (OK 5 2*37 > ) + ». TTK_5._3 *C7 ') > 

MPPfc = C T. 0DG/0GMK4_2_2 3._3:> 

V I Mb _ 1 _ 1 = MPPS*GGMK4 _ { 2_1_3) 

V IM6_2~2 = c MPP6*0GMK.4 *2 ~2_2J3> 

T AUP6 = C MPP6* C T AU_6 -QGFK4_2_3‘.> > 

GFK4_1_1 = C 0GFK.4_ 1 _ 1 + CQGMK4_ 1 _2_1_3*TAUP6> :* 

GFK4_2_2= C0GFK4_2_.2+CQGMK4_2_2_2_3*TAUP6.> ) 

GFK4_2._3 = CQGFK4 _2_3+< QGMK4_2_2_3 _2*TAUP6> ) 

GMK4 __ 1 1 1 _ 1= c 10. 7D0-cQGMK4._t_2._i _3*V I MS_ 1 _ 1.) > 

GMK4_ 1 _2 .1_2* cQGMK.4 l_2._l _2 - C GGMK 4 .1...2 1 „3*VIM6_2_2) ) 
GMK4_2_2_2..2 = '■ 0GMK4 2"'2 _2 . 2 - CGGMK4 2_2..2_3*VIM6_2_2) ) 

TAUP5 3 C 0 . 1 006 745 1 3279 1 7 1 DO* C TaU £+P7K .3 1 ) .) 

GFK3 1 .2 = C FTK. _3 2 - * 1 1 . SDO* TAUP5 ) ) 

GFK3 1 3 = < F T K 3 3 ~ ■ i 2 . 8 “ l )‘ ! '* T nUi*s ; 

GFK.3. 2 1= '. fTK 3 1 - C 8 • 88300* i ALIK5 .) .1 

GFK3 2 2= < TT K 3 2 - U — 1 U'J* 7 AUF5 ’ J 

GFK.3. 2 3= (FTK 3 3+ C u. • ' 71 mj* • ACRS ' 

CL 1 1 2 2 = < 25 ~ < 0 • 33 1 4 555220 75200*35 » • 

CL 11 3 2= -C 35-m O. 33 1-553 -20/52D0>C5 > > 

CL22 2 2= * LS+CO. 0521 :»2a3720Oa9l)O*S5 > i 
CL22 _3 2= ( (.0. 2521 l628372OO8SD0*C5) -S5.» 

RCL _2_ 1 = - C C 0 . 2D0*C5 .) + t. ;..s . 22D0*S5 1 ) 

RCL _2_2 = C CO. 00 1 334 1 4 1 22337500*35 > -< 0 . 0 1 0742827 1 97694D0*CS ) > 
RCL_2_3 = C C 0 . 003 1 75257 73 1 959D0*S5 ) - • 0 . 0 1 8597933 1 4433B0*C5 ) ) 
RCL _3_..l = C 0 . 1 D0+ C (. 0 . 2D0*S5 ■> - < 0 . 22D0*C5 ) > ) 

RCL_3_2= c CO. 00 1 334 141 223875D0*C5 > + C 0 . 0 1 0742327 1 97634D0*S3 ) .) 
RCL_3 3 = <. C 0 . 003 1 7525773 1 959D0*C5 > + C 0 . 0 1 3597933 1 4433D0*35 > ) 
0GMK.2_1_1_2 .2= (6.SD0+C < 36. 53543822S7436D0*CS5*S5'.> ) + C 
44. 1 35794325329700* < CL1 1_2_2*CL 1 1_2._2> > ) > 

QGMK2_ 1 _ 1 _2_3 = C (36. 535433223743600* c S5*C5'.> > + < 

44 . 1 357 948 25 32 3 7 DO > CL 1 1 .. 2._2*CL 11.3 . 2 .» > > 

0GMK2_ 1 _ 1 ...3 _2= GGMK2.1 _1 „2_3 

QGMK2 1 _ 1 _3 .3 = C fa . 800+ C < 36 . 535433228743600* C C5*C5 > ) + C 
44. 1 857343253297D0* C CL 1 1 3.2*CL11 3 .2) .> > ) 

0GMK2 1 2. .12= C5. 1D0+C57. 5D0*RCL ,„2, l ) ) 

0GMK2 1 2 1.3 s • 57. 5D0*RCL _3 1 ) 

QGMK2 1 _ 2 2 1= - C 5 . 1 DO+ c < -3 . 6535433228 743600* c 35*C5 ) ) + < 

4.41 857348253297D0* C CL 1 1 2 2*CL1 1 3 2)))) 

QGMK.3 1 2 2 _ 2 = C C 36 . 5354 53238 74 36D(..4 ( RCL 2 3 *35 .<.)+( 

44. 1857948253237D0* < CC. 1 1 2 2*RCL . 2 2) ) > 

OGMK 2 1 .2 _2_3 = C C 36 . 53S-332237436D0* C RCL . 3 3*65 > ') * ' 

44 . : 35784:32532370* >* * C L 1 1 2 2 * RCL 3 2 .) .» ) 

OGMK 2 1 2 .3 1= - C *. 3. 6~.-354 3822 J7436D0* * C5*C5 .j ) + C 
4. 4 185 734825329 7DO* <. UL 1 1 3 2*C,Lll 8 2) ) ) 

0GMK2, 1 _2 _3 2= *. (36. 535433228743600* i. RCL 2..3*C5) !> + ( 

44. 185794825323700* (CL 11 3 2*KCL,.2_2.) ) > 

OGM K 2 _ 1 '2 _ 3 _3 = C C36. 535438223743600* *. r.CL _3_.3*C5> .>+( 

44. 135794325329 7D0* C CL 1 1 .8 _ 2* RCL 3 .2 ) ) » 

0GMK2 _2_ 1 _1 _2= OGMK 2 _ 1 _2_.2 _ 1 
0GMK2_2__ 1 _ 1 _3 = OGMK 2 .. 1 _2 .3 1 
0GMK2_2 ... 1 _2._ 1 = 0GMK.2 ... 1 _2._ 1 J2 
□GMK2_2_ 1 _2 .2 = OGMK 2 _ 1 .2 .2.2 
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0GMK.2 2 1 2 .;■= ■ _ j - 

0GMK.2 2 1 J 1= ' 5/. iDu*HCL 3 1 i 
0GMK.2 2.1 _ 3 _. 2 = 0GMK2 1 ^23 
0GMK2 .2.1 3 .3= 0GMK.2 *1 2 3 3 

0GMK2 .2 .2 , 1 1 = • 4. 1750o+>. >.<_•. 36535433223 7 436D0* (. C5*C5 > J + (. 

O. 44 1857948253297 DO* -LL 1 1 3 2 * C L 1 1 8 2 ) ) j i 
0GMK2 2 2 1 2= - <. (3. 65354 38228743600* > RCL _ 2 3*C5> > + (. 

4. 4 165794325323 /DO* ' CL 1 1 3 2*RCL 2 2 ) > > 

0GMK2. 2.2 1.3= - ( (.3. 6535433228 7436D0* >. PCL .3 2*C5) ) + ( 

4. 4165794825329700* (. CL 1 1 .3 ,2*RCL .3._2> > > 

QGMK2 .2 _2_2_ 1 = QGMK2_2_2_1 *2 

0GMK2_2_2..2_2 = ( 4 . 1 75D0+ < t ( 2. 1 65340206 1 655700* (. CL22 2 2*CL22 2 25 
> + ( 5 . 36 1 3549929U605D0* (. 35*65 > * .)+•>. '. 86. b8543622S7436D0* ( 
RCL_2_3*RCL„2_.3> > + ( > 44. 1 65794825329700* ( RCl_2_2*RCL_2_2> > + ( 

57 . 5D0* > RCL _2_ 1 *RCL _2_ 1 > > > > > > 

0GMK2_2_2_2_3 = ( ( <2. 1 65340206 1 855700* ( CL22_2_2*CL22 _3_2) > + ( 

5 . 36 1 35499290605D0* < 35*C5 .) > > + ( (. 3S • 535433223743600* ( RC L _2_3* 
RCL_3_3) > + <■. '■ 44. 185794625329700* <PCL_2_2*RCL_3_2:> > + (57. 5D0*>. 

RCL .2.. 1 *RCL „3„ 1 > > > > > 

0GMK2 ,~2 ~2 .3 1= l )GMK2 2...2 1.3 
0 G M k 2 2 J. 2 .1 3 2= 0GMK2 2 2 .2.3 

0GMK2 _2_ 2.3 .3 = ( (. (2. 1 65340206 1 3557D0* < CL22 .3 ,2*Cl22_3 „2> > + (. 

5. 36 1 35499290605D0* ( C5+C5 > > > + ( > 36. 53543622S7436D0* ( RCL _3_3* 

RCL 3 3> >+< <44. 135794625329700* (RCL 3 2*RCL 3 2> > + ( 57 . 5D0* ( 

RCL 3 1 * R C L 3 1 > > > > > 

CF_2= ( ( GFk3 1 2*C5> + >. uFk.3 l 3*85 > ’ 

CF 3 3 > GFKo 1 3*Lb.j - G+ko 1 2*35 .> < 

GGFK 2.1.1 s > FTK 2 1 +F T K 3 1 > 
ijGFK'8 1 2= ( LI - 2+P ! < 2 2> 

(JGFK2 .1.3s (CF o + FTK .2 3> 

0GFK2_2_ 1 3 ( TT K „2_ 1 + ( GFk3"_2_ 1 - ( 0 . 1 D0*CF .3 > ) > 

0GFK2_2_2 S < TTk .2 2+( >. GFK.3 .2_2*C5> + CGFK3 .2_3*S5> > ) 

0GFK.2_2_3 = >. (0. 1D0*FTK_3_1 > + ( ( GFK3_2_3*C5> -(GFK3_2_2*S5> > > 

MPP4 = >. 1 . ODO/OGMK'2_2_2_3_3> 

0 1 M4_ 1 _ 1 = < MPP4*0GMk2_1...2_l_3> 

VIM4_1_2= ( MPP4*0GMk2_ 1 _2_2_3 > 

V I M4~ 1 _3 = ( MPP4*0GMK2_ 1 _2_3_3> 

0IM4 ~2~1 = CMPP4*0GMk2_2_2_l_3> 

V I M 4 ” 2 ” 2 3 < M P P 4 * 0 G M K 2 _ 2 . . 2 _ 2 _ 3 ) 

TAUP4= ~>MPP4* > T AU_4 -QGFk 2._2_3> > 

GFK2.1 _1 3 <. QGFk2_ 1_1 + ( 0GMK.2~ 1 _2_ 1 _3*T AUP4 ) > 

GFk2 3l_2= ( 0GFK2_ 1_2+ ( 0GMK2_ 1 _2_2_3*T AUP4 > > 

GFK2 " 1*3= (aGFk2_l_3+t 0GMK2_1_27.3 - 3*TAUP4> > 

GFk2.2_l 3 (0GFK2_2_1+(QGMK2_2_2„ 1 _3*TAUP4> > 

GFk2*2_2= ( 0GFK2 .2 2+ ( 0GMK2 .2 2 . 2 3*TAL)P4> ) 

GFK2„2_3= (0GFK2 .2_3+(uGMK2_2_2.,3 3*TAUP4> > 


GMK 2 

X 

i 

I 

1 = 

i r^,»4 . 32^0 

- •: OGMK 

2 1 2 1 3*VIM4 

1 1 .) » 




GMK2 

.1 

.1 

1 


- ( QGriK2 

1 2 1 

3*V]>14 i 2> 





Gm 2 

1 

l 

1 

**> — 

-J — 

- i Y- x. 

.1 2 1 

3*V IM4 i _ 3.) 





GMk2 

.1 


jL 

2 - 

>. Li GMK 2 

: 1.2 .. 

2-CQGMK2 122 

3*V IM4 

1 

2 ) 


GMK2 


i 


■ -> _ 

— 

( UGr-1 K 2 

i 2 

3-( i.jGMK2 1 2 2 

3*VIM4 

__1 

__o) 

> 

GMK2 

1 

...l 

.3 

• *‘t _ 

-Zi — 

(.□GMK 2 

t L .j 

3 - ( UGMK 2 1 2 3 

3*V I M4 

.1 

.3) 


GMK 2 

.1 


,_1 

1 = 

- >. OGM K 2 

1 2 1 

3* VIM- 2.1.1 





GMK 2 

1 


I 


(□GMK 2 

i L . 1 . 

2 ~ ( UGMK2 1.2.1. 

3*0 I M4 

2 

2) 

;) 

GHK2 

1 



1 = 

(UGMK2 

^ 

1-(.GGMK2 12 2 

3*0 I M4 

4w 

' 1) 

.) 

GMK2 

.A 

_2 

_2 

J2 = 

< 0GMK.2 

1 ' j 

X X- 48. 

2-> DGMK2 .1 .2 2 

.3*0 1 M4 


_2 :> 
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GMK2.1.2 

w? 

_ 1 

= 

• QGMK2 

1 3 . 

1 -(0GMK2..1.2 

n 

_ 1 

,3*VIM4 .2 

1 ) > 

GMK2 . 1 .2 

-n 

— w 

__2 

= 

( qgmk;2 _ 

1 .2 3 . 

2-( 0GMK2.1 .2 

Jj 

.3*VIM4l2~ 

2) > 

GMK2_2_2 

_1 

_1 

= 

(0GMK2 ~ 

2.2 1 ... 

1 -CGGMK2 .2. 2 

. 1 

3*VIM4 2.” 

1 > .) 

6MK2 2 2 

" 1 

2! 

3 

( 0GMK2 " 

2 _ l 

2 - •: OGNK 2 2 2 

l 

3*V I ri4 2 

2.) ) 

GMH: 2 2.2 


_ 2 

= 

(0GMK2. 

f ^ 

2-C0GMK.2 2.2 


3*V I M4 2 

2.) ) 

D 1 1 NV 1 .. 

1 = 

(’.) 

• 

ODO 






if (GMK.2 

1 

1 


1 1 . ne. 

< • . oDO 

.) then 




D 1 

i N 

V 

1 

1 = (1.0 

OU/bNK 

2 1 1 1 1 1 





endi f 

L 1 1 2 _ 1 = cDlINV 1. 1112* 

Lll 3 1= (DllNV.l 1*GMK2 l .1.1. 3) 

L11D1 2 2= < GNK 2 1 1 2 .--'GMKl 11 1 _*Lll _ 1 .) .) 

LI ID 1.3 2= ( GM K 2 1 1 2. 3-1 UMK2. i 1 1 .3 *L 1 1 _ 2_ 1 .) .) 

L 1 1D1 3 .3= (6iNK2_l .1 _3 _3-Cbfli<:.2 1.1 1 J*Li 1 ,3.1)j 

D1INV 2_2= O.ODO 

it ( LI ID 1.2.2 . ne. O.ODO.) then 

D 1 INV_2_2= C 1 . ODO/L1 lDl _2_2.) 

endi f 

Lll _3 .2 = ( D 1 INV_2_2*L11D1 .3.2) 

L 1 1 D 1 .3.3 = < L 1 1D1 _3_3-CL 1 1 3 2*L11D1 3 2).) 

Dll NV.3.3 = O.ODO 

if (L11D1_3_3 .ne. O.ODO) then 

D1 INV_3_3= ( 1 . 000/ LI 1D1 _3_3) 

endi f 

D1L21_2_1= ( GMK2. 1 „2_2_ 1 - < GMK2 _ 1 _2_ 1 _ 1 *L 1 1 2.1)) 

D1L21 _2_2 = ( GMK2 .. 1 ...2 .2 2 - ( GMK2 " 1 ~ 2~ 1 ~2*L 1 1~ 2~ 1 ) :> 

D1L21_3 .1 = c.GMK2 .1.2..3 1 - ( CD1L21 2 1*L11 3 2)+<GMK2 1 "• 1 1* 
Lll 3 D)i " ' ~ - - “ ~ 

D1L21 3 2= c GMK2 „ 1 _ 2 3 2-CCD1L21 2 2*1.11 3 2)+«GMK2 1 1 "'* 

Lll 3..1))) " “ 

L21 1.1= (Dll NV 1 1*GMK2 1 2 1 1 > 

L2 1 __ 1 .2 = (Dll N V 2 2*0lu2i 2 1.) 

L21 1 .3= »Dll NV 3 3*»>1L21 1) 

L21 2 1= (Dll N V 1 1*GMK2..1 2.. 1.2.) 

L21 .2 2= <DlI NV 2 2*UIL21 2 2.) 

L2l 2 3 = (Ui I NV .o 3*1)1 L21 3 2 ) 

L22D2 1 1= CbMK2 2 2 1 1 - ( (D1L21 3,1 1 Jit! t D1L21 2.1* 

L2 1 . 1 _2 .» + ( GMK.2 1.2 1, DL21.1 1))).) 

L22D2 .2.1 = CGMK2 .2.2.. 1 . 2-( (D1L21 .3 ,1*L21 2 .3> + ( (D1L21 .2.1* 

L21 .2 .2 .' + ( GMK2 _ 1.2..1 UL21.2..1))) ) 

L22D2..2_2= CGMK2.2. 2.2 _2-( (D1L21 .3 .2*1.21 .2.3) + ( ' D1L21 2_2* 
L21 _2_2) + CGMK2. 1.2.1 _2*L21 1 ) ) ) .) 

D2INV .1.1 = O.ODO 

if <L22D2.1_1 .ne. O.ODO.) then 

D2 I N V. 1 _ 1 = Cl. ODO / L 22D2 _ 1 _ 1 ) 

end l f 

L22 _2_ 1 = CD2INV _ 1 _‘l *L22D2..2_1 ) 

L22D2.2 .2 = C L22D2.2.2 - ( L22..2. 1 *L22D2 2.1 ) > 

D2INV_2_2= O.ODO 

if (L22D2.2.2 .ne. O.ODO) then 

D2 1 NV .2 _2* ( 1 . ODO / L 22D2 2 2 ) 

end i f 

LG 1 1 2 1 = L 1 1 2 1 

lU 1 1.3 1= Lx i l 

LG 1 1 3 2= Lll 3 2 

L_ U 1 1 i = (_ L x 1 1 
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LO 2 _ 1 _1 . - 
LO_2_ 1 _ 1 _3 
LO_2_ 1 _2_ 1 
LCL2_1_2_2 
LO_2 _1_2_2 
L0_2_2_2_l 
D I _ 1 1 l_l 
DI_i .1 .2 2 
D I 1 1 3 i 
DI_2..2_.l 1 
01.2 2 2 2 
CL 111 1 = 

CL 1 i 2 1 = 
CL 22 1.1 = 
CL 22 .2 1 = 
RCL 11= C 
RCL 1 2=i 
RCL .1.2= v 
RCL 2 1= •: 
RCL ..2 .2= < 
RCL_2..3 = < 
CLlloi 1 1 
CL 1 1 D 1 _ 1 _2 
CL 1 1 D 1 _2_ 1 
CLUD1_2_: 
CL 1 1 D 1 _3_ 1 
CL11D1_3_2 
CL22D2..1 
CL22D2_ 1 J 
CL22D2_2_3 
CL22D2 2 •• 
RCLD 1 _ T_ 1 : 
RCLD1 1 .2 : 
RCLD 1 . 1 ...3 = 
RCLD 1 2 1 : 
RCLD 1 2 2 - 
RCLD1 .2 3= 
OGMK 11 1 . 

GGMK1 .1, 1 
OGMK 1.1 1 

OGMK 1 _ 1 l 
OGMK 1 . 1 .. 1 . 
OGMK 1 _ 1 _ 1 
OGMK 1 .1_1. 
OGMK 1 _ 1 _ 1 . 
OGMK 1 _ 1 1 

cclIioi. 

OGMK 1 _ 1 _2. 
OGMK 1 _ 1 _2. 
0GMK1_1_2. 
OGMK 1 _ 1 _2 
OGMK 1 _ 1 _2. 

RCLD 1 _ 1 
OGMK. 1 _ 1 _2. 

RCLD 1 „2 
OGMK 1 2. 1 


:= L21 .. 1 _2 
= L21._l_3 
= L2 1 _2_ 1 
= L21_2_2 
:= L21 _2_3 
= L22 _2_ 1 
= GMK2.1 ..1.1 .1 

= LllL'l 2 2 

!= L11D1 3 _3 
= L22D2 1 l 
:= L22D2 2 2 


CC4-CLU 1 1 

C 34+ C LU 1 1 

1 L 4 - 1 L Li ._ 2 

< 34+ ' L U 2 2 

< 1 . w> L • • *■ C 2 1 1 

< 1 . SD«»*C4 ) -*■ 


2 1 *34 i 
2 1 ,i > 

2 1 *34 .i i 

... . -<-04 i j 

2 1 > i- C ' L i J 

i ( LU 2 1 1 


I'LL! 2 1 1 342-;. -• LO 2 1 


< LU 2 


1 


i 


1 ') + i 

L 

u 

2 

1 . 5 DO 

*• 

54 

) 


( U.C) 2 

1 

1 

2 

LO . 2 . 

1 


i 



> + f . L 

U 

2 

... 1 

< CL1 

1 


i 


1 

*D I , I ... 1 


1_ 

1 .) 

-< DI 


i 


1 


2_2*S4:> 




<CL1 

r 


2 


1 

*DI_1 _i 


1 

n 

< DI_ 

i 


i 



_2*C4 ) 




< DI_ 

i 


i 


1 

.1*L0_1 


1_ 

n 

■_j _ 

(01. 

i 


i 


2 

_2*U0_1 


1, 

_ 

‘CL2 

2 


i 


1 

*Q I _2J2 


1_ 

1 .) 

-CD I 




2 


2.2*54) 




< CL2 

2 


2 


1 

*D I _2_2 


1 _ 

i :> 

> < D I 

2 


2 


2 

:_2*C4> 




CD I 1 


r 


r 


l*PCL w l 


1 ) 


CD I 1 


i 


2 


2*RCL 1 


2) 


< DI . 1 


i 


o 


3* RCL . 1 


3.) 


< 0 i 1 


i 


1 


1 *-HCL _2 


1 ) 


< DI 1 


i 


2 


2*HCL 2 


2 ) 


(Oi l 


i 


• *< 
o 


o'*KLl _2 


3) 


1 = c 

a 

i 

0 

■ 

o 

D0+ ( (.CL 

L 

1 

i 


1 2= : < CL 1 1 2 1 *CL i 1 D 1 

1 3= ••iCLllDl 1 1 * L U 1 

2 1 = OGMK 1111 _2 
2__2 = < 410. 0D0 + < t CL 1 1 2 

_2~3 - C CCLl 1D1 . 2 . 1*LG .1 
~3 ~1 = OGMl-: 1.1 1.1 ..3 

3_2 st OGMK. 1 _ 1 . 1 _2 _3 


1 1 1*04.) LO 2 12 1*S4 :•.).) 

:*C4 i - ' LO .2 1 2 2*S4.) > .1 

2 3*S4).i 

1 2 1*C4) » -vl.5D0*CLU_l 1 ) ) 

*34.) + -. LO 2 1 .2 _2*C4 l ) > 

2 3*04 .) « 


1 :> 
2 !) 


1 *CL 1 ID 1 1 1.)-<CL11DI 1 2*S4-.i :> .) 

i i :> + c c l i i d i _ i 2*C4 < .» 

1 3 1 .< + < o L i 1 D 1 l 2*1.0 ... 1 .1 3..2) > 

1 *CL 1 1 D 1 2 1 > * '.CL 1 1 D 1 2 2*C4.» .> ) 
1 _3 . 1 >+CCLl 1D1 .2_2 *lQ_ 1 .1 ...3 _.2 .i .' 


3_3 = <41 0 • 0 DO * <. D I _1_1_3 .3+ <. C CL 1 1D1 _3_ 1 *L0._ 1 _1 _3 _ 1 ) + 


3 ..2*L0_ l .. I _3_ 2) ) :* :> 

V'l = <. < CL1 1 _ l _ 1 *RCLD 1 1 _ 1 ;> - C RCLD 1 _ 1 _2*S4 > > 

i _2 = <: < cli i_i_i*rcldi_2_i :> -<:rcldi_2_2*S4:i :» 
2_ i = <: (CL11_2 _i*rcldi_i_1) + <:rcldi_i_2*C4:) 
2_2= <: C CL 1 1 _2_ 1 *RCLD 1 _2_ 1 > + C RCLD 1 _2_2*C4 ) ) 
3_ 1 = < RCLD 1..1 _3+C <: L0_r_l_3_l*RCLDr_l_l :) + CL0 
2) > .) 

3_2 = cRCLDl ..2,.3+<: CL0_l_l._3_l*RCLDl_2_.l :i + CLO 

2 ) '.) ) 

’l_l= 0GMK1 .1 .2 .1 _1 
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OGMK l 2.1 12= OGMK 1 I 2 2_1 
OGMK 1 _2 .1.1 3 = UGMK 1 . l 2.3 1 
OGMK 1 _2_ 1 _2_ 1 = OGMK 1 _1 .2 1 _2 
OGMK 1 _2_1 _2_2 S OGMK 1 ”l 2 22 
OGMK 1 _2_ 1 _2_3 = OGMK 1 _ 1 2_3 _2 

OGMKl_2_2_l_l = ( 115.0D0+C C < CL22 _ 1 1 *CL22D2... l _ 1 ) - ( CL22D2_ 1 _2*S4 ) 
+•( C RCL_1_3*RCLD.1_1 _3) + C < .RCL _ 1 _ 1 *RCLD 1 _ 1 _ 1 ) + ( RCL_1_2*RCLD 1 _ 1 _2 

i ) ) 

OGMK 1_2_2_ 1_2 = C < ( (CL22_2_1 *CL22D2_ 1 _ 1 ) + < CL22D2_ 1 _2*C4 ) ) + C ( 

RCL_2_3*RCLD 1 _1_3)+C ( RCL ,_2._ 1 *RCLD 1 _ 1 1 .» + C RCL_2 ,2*RCLD 1 _ 1 _2) ) ) 
14. ODO ) 

OGMK 1 _2 _2_2_ 1 = OGMK 1 _2 .2 . 1 _2 

OGMK 1 _2_2_2_2 = ( 3 16 . ODO+ < c <CL22_2_ 1 *CL22D2._2_ 1 ) + ' CL22D2„2 _2*C4 ) 
+ < (RCL 2 3*RCLD1 _2_3) + < < RCL _2_ 1 *RCLD1 _2_I ) + < RCL ,.2_2*RCLD 1 _2_2 
) ) ) 


CF _ 1 - 

c CGFfc 

2 ... 1 _ 1 *C4 -• 6FK2 .1 

2^-SA ) '' 


CF 2 = 

y y b F K 

.2 1 1 ★S4 > + *•. 13 F K.2 1 

2*C4 .) .) 


QBFK 1 

.1.1 = 

t CF 1 + > 4 I U « • Db * 'AK 

1 l .> ) 


Ub FK 1 

1 2 - 

f L'i- ,-i+i A ; :.i , .. ;i J.J frA. 

1 2 ) ) 


GbrV 1 

L 2j = 

6 F K 2 I 3 + ■ .a 1 O . a y 

.j * A i*. 1 


IJ 3 F K’ 1 

1 i 

2 1 = 

' T 7 K 1 1 r i (. 1 . TjYJ 0 fri 

~ F 'J. i i“ » 

•■orK 2 2 1 ^C4 ' - < bF K2 _ 2.„ 2+S4 ) ) 

) J 

OGFKl 

2 2 = 

<. IT K l 2+1. (brK2 4 

l^S4.J 

T ( bp K 2 2 2*C4 5D0*CF , 1 ) 

OGFK 1 
D1INV 

2 _ 2 = 

- 1 1 1 = 

(GFK2.2 3+TTK I 3. 
0 . ODO 

f 


if CGMK4 1 

1 11 . ne. O#0D0> 

then 



Dll N V _ 1 _ 1 * Cl. 000/ GMK4_1 _ 1 _ 1 _ 1 > 

endif 

D1INV_2_2= 0. ODO 

if (QGMK4_1_1_2_2 . ne. O.ODO) then 

D1 INV_2~2= C 1 . 0D0/0GMK4_ 1 _ 1 ._2_2) 

endif 

L 1 1_3_2 = CD1INV _2_2*0GMK4_1_1_2._3) 

L 1 ID 1 _3_3 = ■ OGMK 4 _ 1 _ 1 __3_3 - ( L 1 1 „3 _2*GGMK4_ 1 _ 1 _2_3 ) .> 

D 1 INV3 _3 = O.ODO 

if (L11D1_3„ 3 . ne. O.ODO) then 

Dll NV 3 3 = ( 1 . ODO/ L 1 1 D 1 3 3 ) 

end i f 

L. 2 1 _ 2 _ 1 = (Dll NV 1 1 *GMK4 .1 .2 1.2) 

L22D2 2 .2= (GMK4_2_2 .2 2-CGMK.4_l 2 1 2*L2 1 ...2... 1 ) ) 

D2INV 2.2= O.ODO 

if (L22D2_2..2 . ne. 0.0DO> then 

D2 1 NV _2_2 = Cl. ODU / •... 2202 2 2 ) 

sna l f 

LO .1 1 .3 2= L 1 1 3 2 

LC 2 12 1 = L21 .2 1 
Dll 1 1 1 = GMK. 4 1111 

01 1 1.2.2= 0GMK4.1 1 2 2 

D I _1 ... 1 .3 3= L 1 ID 1 _3 3 
DI_2 .2.2.2= L22D2..2 2 
RCL _ 1 _ 1 = - ( LO _2... 1 _ 2. 1 ) 

RCL 21= (LO .2 .1 _2_.1*CS > 

CL1 1D1 11= ( D I .. 1 _ 1 1 l*Cb) 

CL 1 1 D 1 1 2= -CD I .1.. 1 ,2. 2*SS) 

CL11D1 2 1= ( D I 1_1 1 1 *S6) 
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CL i 101 2= 'L>1 l l . _*Co ' 

C L 1 1 D 1 :• _' = '■ D i l l i.' ... + L U l l fa J .» 

CL 2202 1 2= -<.DI 2.2 2 2*36.) 

CL2_U2_2 .2= <DI 2 2 2 2 ■*!..**> 
n C L 0 1 11= (01 1 1 1 i I 1! 

RCL01 1 2= -i 1> 2D0* < 0 1 : i _*LU 1 1 5 2 ') 

RCLi.'i 2 1= CD I l 1 1 l *hLL 2 1 ) 

RCLL1 3 1= t 1 . 2DO** ( 0 1 L 1 1 1*C6.>.) 

RClD1_3 2= ~(1.2D0*(DI 1_1.2 2*S6) ) 

0GMK1.1 1_1 1= ( OGM K i l l 1 1+C(CL11D1 1 ... 1 *Cfa ) - c CL 1 1 D 1 _ 1 2*So> ) 

) 

OGMK1 .1 .1 1_2- (OGMK 1 i .. 1 . 1 _2+ ( < CL 1 1 D 1 .. 1 . 1 *36) + (Cl 1 1 Dl_l_2*Co) > 

) 

OGMKI _l_l .1_3 = (OGMK 1 _1 ..1 _1 _3-KCLl 1D1 _1_2*LU 1_1 3 _ 2 !> ) 

OGMK 1 ~ 1~ 1 2~ 1 = QGMK 1..1 .1 _1 ..2 

□GMK 1 _ 1 _1_2._2 = ( QGMK 1 .1 . 1 _2 .2+ ( '■ CL 1 1D1 2 „l*S6)+< CL 1 ID 1 _2_2*C6'J ) 

;> 

□GMK 1 _1_1 _2_3= (QGMK 1 _1 .1 .2 3+(CL 1 1D1 ...2_2*L0_1 1 3 2')') 

QGMK 1 _ 1 _ 1 .3 .. 1 = 0GMK1.1 1 .1.3 

OGMK 1 .1. 1 .3 2- OGMK 1 l 1 .2.3 

QGMK 1 1 l .3 3= (OGMK i 1.1.3 3+< DI 1 _1_3 3+ > CL 1 ID 1 _3_2*L0_ 1 _ 1 _3 ..2 .' 

.) > 

OGMK i_l 2_1 .. 1 = (OGMK 1 1 ..2 l .! + ( < RCLD1 .1 1 *C6 J - ( RCLD 1 _ 1 _2*S6> ') ) 

OGMKi 12 1 2= (OGMK. 1 L - 1 2i-(RCLDl 2. l*Cfa.>.) 

OGMK 1 1 _2 1 .3= ((RCLD. J 1 *Cfa > - • RCLD 1 3 2*36) > 

OGMKI 1 2 2 1= (OGMk 1 i. 2.2 li-( ( RCL01 1 1 *S6) + ( RCLD 1 _ 1 _2*C6 > .) ) 

OGMKI 1 2 2 2* (OGMKi l 2 2 2+ (RCLD l 2.1*36:).) 

OGMKI 1 2 2 3= ( ( RCLD 1 1*36) f( RCLD I 3 2*C6> ) 

OGMKI 1 .2 3 1= (OGMKI i. 2 3 1 + ( u_0 1 1 3 2*RCL01 .1..2) -U.2D0* 

D I 11 3 3.))) 

□GMK 1 1 _2 3 2= OGMKI l 2 3 2 

OGMKI 1.2 3.3= ( L0_1 1 3.2*RCL0l .3.2) 

OGMK 1 ,.2_ 1 _ 1 _1 = OGMK 1 .1. 2 .1 . 1 
OGMK 1 _2_. 1 ._ 1 _2 = OGMKI .. 1..2_2_1 
OGMK 1 _2_. 1 _ 1 . 3 = OGMKI _i_2_3_l 
OGMK 1 _2_. 1 _ 2 .. 1 = OGMK 1 ... 1 _2. 1 _2 
OGMK 1 ~2_ 1 _2_2= OGMKI _1_2_2_2 
OGMK 1 _2_ 1 _2_3 = OGMK 1 _ 1 _2 _3_2 
OGMK 1 _'2._ 1 _3 _1 = OGMK 1 _ l _ 2 _ 1 __3 
OGMK 1~2" 1~3_2= OGMKI _1_2_2_3 
OGMK 1 _2_ 1 _3 _3 = OGMK 1 _ 1 _2 ...3_3 

OGMK 1 ~2_2_ 1 _ 1 = ( OGMK 1 _2_2_ 1 _ 1 + ( (Cl. 44D0*D I ._ 1 _1 _3_3 !> + ( (RCL_1_1* 
RCLD1 _ 1 _ 1 5 - ( 1 . 2D0* ( LG . 1 . 1 .. 3 _2 *RClD 1 _1 _2) ) > ) - ( CL22D2 , 1 _2*S6 > ) .» 


OGMKI _2, 

2 

1 2= (OGMKI 2 2 

1 _2 + ( 1 

(CL22D2.. 

1 2*C6)+(RCL 2 

I* 

RCLD 1 . 

1 

_i:> ) ) 





CGMk 1 _2 

Z 

1 .3 = (14. 0D0+ ' 1 

• 2D0* < *. 

RCLD 1 _1 

. 1*C6) -(RCLD1..1 

2+56 ) ) } > 

OGMKI . 2 

Z 

.2. 1= OGMKi 2 

l z 




OGMKI .2 

z 

2 2- <' OGM s*. i z z 

2 2+(< 

' CL22i.>2 

2 2*Cfa.) + ( RCL. 2 

1* 

RCLD 1 

z 

1 •> - J > 





□GMK 1 2 

z 

2 3= ( ( 1 . 2Dd* < « 

Cl. D 1 

1 ') '} - 

34. fa DO ) 


OGMKI 2 

z 

3 1= OGMKI 2 2 

1 o 




OGMK 1 2 

z 

3 2= OGMKi 

Z Z 




OGMKI 2 

2 

3 3= (440. 0 00+' 

1 . 2D0* 1 

l *■ KCLDl 

3 1 *Cfa ) - ( RCLD 1 

3_2*S6 ') ) ) 


.) 


CF . 1 = ( ( GFK4 1 .1*06) -UJbFK.4 .1 .2*36) .) 

CF 2= (» GFK4_1 1 *S6) +(GGFK4 1 2*C6 > j 
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OGFK 1_ 1.1= < CF ,1-HJGFKl 1 1 > 

OGFK 1 .. 1 _2= ' CF..2+0GFC, l 12.) 

OGFK 1*1 _3 = (OGFK l _ 1 _3+OGFK: 4_ 1 .3) 

OGFK 1 ~2_. 1 3 <.QGFK1 <. <.QGFK4 2 1 *C6) - • 1 GFK4 _ 2*So> .* - ( 1 . 200* 

UGF K4~_ 1 3 > > > 

UGFc 1 _2..2 = C0GFK1 2 „2-«- • ' jrK4 2 2*C6) + CUGrK. 4 .. l*Sb.).i, 

OGFK 1 2 3= 1 OGFK 1 2 0 + 1 \it~ K-. 4 — ' i « ^.DU+Cr * - ■' ) -) 


SOLVE F-'uR UDUT (ACCELERATIONS > 


Dll NV 1 1= < 1 . ODO/OGMK 1 1111) 

L 1 1 2 1= (Dll NV 1 1*UGMK1 1 1 2 D 

Lll 3 1= 'Dll NV 1 1*QGMK1 1 1 _3„ 1) 

LI 1D1 _2_.2» (GGMKl .1.1 .2 2-0.11 2 . i*ObMKl _l_l _2_1.) ) 
L 1 ID 1 _3 2= ( OGMK 1 .1 1 _3 _ 2 - « . L 1 1 _2 1 *GGMK 1 1 _ 1 _3. 1 ) ) 

L 1 1 D 1 ~3_3 = t; OGMK 1 _ 1 _ 1 .3 _3 - ( L 1 1 3,1 *GGMK 1 .. 1 .1 ,.3_ 1 ) ) 
D1 INV_2 .2= C 1 . ODO/L 1 1D1 _2_2> 

Lll 3 2= (D1 INV..2 „2*L1 1D1_3 .2) 

L 1 1 5 1~3 _3 = CL 1 1D1 „3_3-< LI 1 „3_2*L 1 1D1 _3_2) ) 

D 1 INV 3_3 = C 1 . ODO/L 1 1D1 „3„3> 


D 1 L2 1 _2._ 1 = c OGMK 1 _ l _2 _2_ 1 - - Lll _2_ 1 *GGMK 1 _ 1 _2„ 1 _ 1) ) 

D1L21 2 ~2 = (0GMK1 .1 _2 .2 _,2-< LI 1 _2_1*GGMK 1 1 2 .1 _.2i> ) 

D1L21 ~2~3= < OGMK 1 * 1 _2 .2 ,3 - ( L 1 1 ,.2_ 1 *OGMK 1 _ 1 _2 _ 1 _3) ) 
DlL2ll3_l» C0GMK1 _1_2_3 _1 -( * D1L21 .2._1*L1 1 _3_2)+<Ll 1 _3._1* 
QGMK1..1 _2_1_1 ) > ) 

D1L21_3_2* C OGMK. 1 .1 „2 3 2-C (D1L21 2 /2*L11 _3 2) + C L 1 1 _3 _ 1 * 
ogmk 1 1 _2 _ 1 .. 2 :> :> :> 

D1L21 .3 3 = C0GMK1 1 2 3 3-C<DlL21 2.3*1.11.3 2)+<Lll 3.1* 
OGMK 1 1 2 ._ 1 _3 ) ) ) 

L21 1 1= (Dll NV 1 l*LlGMf 1 12 1 1) 

L21 1 2= CD 1 INV ... 2*LK L21 0 1) 

L21 1 3= 'Dll NV 3 3*DiL2l 3 1) 

L21 2 1= 'Dll NV 1 1 *OGMK 1 . I 2 1 2) 

L21 2.2= (DlINV 2 2*DlL21 2 2) 

L2 1 2 3= i.DUNV 3 3*D 1 L2 1 3 2) 

L21 3 1= 'Dll NV 1 i*U6MKl 1. 2 1 3) 

L21 3 2= (.Dll NV . 2 2 * D 1 L 2 1 2 3 » 

L21 3 3= (Dll IW 3 3*D 1 L21 3 3) 

L22D2.11* (0GMK1 _2„2.. 1 _ 1 -< CD1L21 3 ,1*L21 1 .3) + ( CDl L2 1 _2 
L21 " 1 ..2) + C L21 1 .. 1 *OGMK 1 1 ...2 . 1 _ 1 ) ) ) ) 

L22D2l2_l = (0GMK.1_2_2_.2_1 -< CD1L21 _3_1 *L2 1 _2_3 ) + ( <D1L21_Z 
L21_2 2) + ( L21 _2_ 1 *QGMK 1 _1_2._1 _1 ) ) ) ) 

L22D2_2~2= C OGMK 1 _2_2_2 _2 - C CD1L21 _3...2*L21_2 .3) + C CD1L21 Ji 
L21 _2~2) + < L21 _2_ 1 *OGMK 1 _ 1 _2 _ 1 _2) ) ) ) 

L22D2_3_ 1 = C OGMK 1 _2_2_3 _ 1 - (7 D 1 L2 1 _3_ 1 *L2 1 _3_3 ) + C CD1L21 
L2ll3~2)+CL21 _3_1*0GMK1_1_2_1 _1 >) ) ) 

L22D2_3„2 * C OGMK 1 ...2_2._3_2 - (7 D 1 L2 1 _3_2*L2 1 _3_ 3 ) + C C D 1 L2 1 
L21 ~3_2) + ( L21_3_1*0GMK 1 .1 _2_1_2) ) ) > 

L22D2_3~3 = C OGMK 1*2.2 .3 ,.3-C (D1I.21 _3_3*L21 _3_3> + ( (D1L21 
L21 _3_2) + ( L21 .3 ,1*0GMK 1 . 1 _2 1 .3) ) ) ) 

D2INV .1.1= Cl. 000/ L22D2 .1 . 1 ) 

L22..2 1= CD2INV 1 1*l22D2 2 1) 

L22 .3 1= C 02 INV 1 1*L22U2 3 l > 

L22D2 2 2= C L22D2 2 2-< L22 2 l*l.22U2 2 1 .» ) 

L22D2 3 2= (L22D2 3 2--L22 2 l*L22D2 3 1 ) ) 

1..22D2 3 3= ' L22D2 3 3-0.22 3 1*i.2wD2 3 1 1 ) 


. 1 * 
. 1 * 
2 * 


l* 


2 * 


3* 
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D2INV .2 _2» <. I.ODO/L22D2 2_2) 

L22 _3 2= (D2INV_2_2*L22D2 3_2:> 

L22D2~3_3 = ( L22D2_3_3 - < L 22_3_2*L22D2 _3_2> .> 

L0_l_l_2_l= 111 _2_ 1 
L0_ 1 ,l_3_i = LI 1_3_1 
LQ _ 1 _ 1 _3_2= L 1 1 _3_2 
LO _2_ 1 1 _1 = L 2 1 _ 1 1 

LQ .2 1 1.2= L21 1 2 
LQ .2.1. 1 3 = L21 . 1 3 
LO 2_1_2 1= L21 2 1 
LO 2 12 2= L2 1 2 2 
LO 2 1 2 3= L21 2 3 

LQ 2 1 3 1= L 2 1 3 1 

LU 2 1 32= L21 3 _ 

LO 2 1 3 3= L21 3 w 

L0_2 221= L 22 2 1 
LO 2 2 3 1= L 22 3 l 

LO 2 2 3 2= L22 3 2 

D I l 11 1 = GGMK 1.1 1 1 1 

D I _ 1 _ 1 . 2 2 = L 1 1 D 1 _ - 
D I 1.1 .3 .3= Ll 1D1 3 3 
D I 2_2_ 1 _ 1 = L22D2 1 1 

DI~2 2 2 2= L22D2 _2 _2 
DI_2_2_3_3= L22D2_3_3 

W6~2= t: LO_l _l„2_i*0GFKl_,i ..1 ) -OGFk. 1 _1 _2i> 

W6_3= C ( CL0_r_r_3 _l*OGFKl_l _1> -OGFK 1_ 1 _3) - < LO_l _ 1 _3_2*W6_2) ) 

Wb 4 = I C (. ( LO ,2_ 1 1 1 *OGFK 1 _I _ 1 1 -OGFK { 1 ) - ( LQ_2_ 1 _ i _2*W6_2 ) ) - 

(LO 2_l_l_3*W6_3:i) 

W6_5 = ~<: <: ( ( (L0_2_1_2._1*0GFK1 _1_1 ) -OGFK 1 _2_2) - ( LO_2_ 1 _2j2*W6 _2') ') - 
(LO 2_ 1 2_3*Wb _3 > 1 - ( LO _2__2 _2._ 1 *W6_4 .) .) 

W6_S=~< ( ( < ( (LO 2 _1_3 t*OGFKl _1._1 !> -OGFK 1 _2_3) - ( L0_2_ 1 _3_2*W6 _2> ) - 
<; l 0_2_ 1 ,.3_3*W6~3> ) -< L0_2.. 2~3 l*Wb_4) ) - (. L0_2_2_3_2*W6_5 ) > 

V6_ 1 = ( -OGFKl_l_l/DI _ 1 „ 1 _ 1 _ ll 
V6 .2 = t Wb .2/ oT . . 1 .. 1 .2 2.) 

Vb _3 = i Wb . 3 / D I 1 _1 .3 . '3 .) 

Vb 4= f ,*Jb 4/ 0 1 2 2 1.1.) 

Vb 5= 1 Wb b / D i 2 2 2 2.' 

Vb b = ( Wb b / 0 1 3 3 .) 

SOLNb b= Vb b 

SOLNb 5= v Vb 5-aU 2 2 3 2*3ULNb o .> ) 

SOLNb * = t. (Vb 4-1LO 2 2 2_l*SULNb 5.3-aD 2 2 3 1 *3QLNb . 6 ) ) 

SOLNb 3= (((Vb 3-i.LO 2 1 1 3*30 L No 4.)') - (LU 2 . 1 2 . 3*30LNS_5.> ) - 

( LO ,.2_ 1 _3__3*SQLNb _6 ) < 

SOLNb 2= ( ( ( (Vb_2-( LG_l 1. 3_2*SuLNb .3) .» -(L0.2..1..1 _2*SQLN6_4> .) - 

(L0l2_l_2_2*S0LN6_5j , -« loI2_1 _3 2*S0LN6_6> > 

S0LN6_ 1 = ( ( ( ( (VS_1 -(LO .1 .1 j2_ 1 *S0LNb _2> ') -(LO _1_1 .3 ,l*S0LNb_3) ) - 
(LO 2 1_1 l*S0LNb. 4 ) > -< LO 2_1_2 ... 1*S0LN6_5:> ) - <. L0_2_ 1 _3_ 1 * 
S0LNb_b ) .• 

UDOT_ 1 = S0LNb_4 
UD0Tl2= S0L.Nb_5 
UDOT _3 = SOLNb_b 
UDOT ~8 = S0LNb_ 1 
UD0T_9= S0LN6_2 
UD0T_10= SOLNb _3 

ALPH2 1= ( (UDOT 1*C4.) + ' UD0T 2*S4) ) 

ALPH2 2= ((UDOT 2*C4) - ‘ UDC)T_1 *S4.) .) 
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ACC2 1 = < <S4*> UDO T ? l . 3D0*UD0T 1 ) ) ' + <. L4 * ( UDu T d - < 1 . 5D0* 

UDOT 2 .) > .1 > 

ACC2 2= C (C4*CUD0T._9+C 1 . 5D0*UD0T 1 )).< LIDO T B-U.5D0* 

UDOT „2 :> > ) ) 

UDGT_4=( taup4- C C udot _3+ ( < a 1 Dh 2 _ 1 *v im4 _2._ 1 ) + ' alor.2 _2*v im4 2 2) > > + 
C ( udot _ 10*v im4_l_3) + ( (ac c2_ 1 *v im4_ 1 _1 .> + < a c c2_.2*v im4_ l _2 .> ) > ) ) 
ALPH3_2= ( : (AIPH2_2*C5> -t 55*< UDOT _3+UD0T_4> > } 

ALPH3 3= < < ALPH2 2*S5 > + C C5* < UDOT " i-^UDOT _4 ) ') ) 

ACC3._ 1 = (ACC2_1+C0. 1 DO* (UDOT _2+UD0T _4.) ) ) 

ACC3_2= <: (. ACC 2 2*C5 > - CSS* (UDOT lO-CO. 1D0*ALPH2. 1 ) ') > ) 

ACC3 _3 = C C ACC2 _2*S5 ') + < C5* C UDOT 1 0 - C 0 . 1 D0*ALPH2._ 1 ) ) ) ) 

UDOT 5 = ( tauo5-C ( CO. 0070472 1 6344542 JO** 1 oa3 3) -Caloh2..1 + C 

0. 04127655290446d0*aloh3 2) > ') -C C 1 . 1 5775697 1 7 1 046d0*a c c3 .2 ) + 

C 1 . 273532668881 5 ldO*ac c3 3) ') ) J 
ALPH4 _ 1 = (• (UDOT . 1*C6:> + < UDOT 2*Sfa ) ) 

ALPH4 2= ( i UQG I 2*C6 ) - * UDU T 1 *96 .> > 

ACC4 _ 1 = i (.UDOT 9*S6» + (Cfo*(UDQ7 S+ «. 1 . 2 bO*uOOT 3) > > ') 

AUC4 2- t < UDU i -?*U6 i - ■ — <6* c i mi! i o ' a < l . JUO*ULOT 
AUC4 3- , UDOf iO-t. l. 2l?0*UDUT l * > 

UDO r 3- . AijPS - (. 1 AUC4 l * V i .'16 i i .) + < UDOT 3-M.ALPH4 2*VIM6 2 2.> < ' * 

ALPH5 2= ( i ALHH4 2*0 /.>*■' 37 #■ < Ui,'U 7 H JDU i o.i .5 .i 

AL Hho .j z ( i w / *■ l uL'i j . jTljJl.. • r.* i .< - c i”t A ,2 * 7 / t .) 

ACCo b = * < . ALrC4 7 > *+■ (. ACU A LJ* T b / ' j .* 

ALLj b — (. i. AULA 3*1./ .) — ( AUU4 u. A b / .) ' 

UDOT 7= (.1 AUP7 - (ALPHA 1 - CO. 24568092 7a96o22DG*ACC5 .3) ') ) 


! USED 27.06 SECONDS CPU TINE. 

28805 BYTES OF EXPRESSION STORAGE. 

EQUATIONS CONlAlN 514 ADDS/6US TRAC fs/ NEGATES 

733 NULTI PLIES 
26 DIVIDES 
577 ASSIGNMENTS 


END REGION 3d_dvnamics 
DYNAMIC sd . states 
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= 

QDOT 
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□ DOT 
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QDOT 
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7 f 
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QDOT 
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□DOT 
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a 
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= 

Li DOT 



Q I O' = UDU T 10 

U 1 ' = UDOT 1 

U 2' = UDOT 2 

U 3 ’ = UDOT .3 

u 4 ' = udot _4 

u_5 ' = udot._5 

U 6' = UDOT .6 

u” 7 ’ = UDOT _7 

U_S ' = UD0T..S 
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U = IJDUT 3 

U 10’ •= •-•Du I 1 

end. nan = < system time.at.J0> 
END DYNAMIC sd states 
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DADS: Dynamic Analysis and Design System 


PRODUCT DESCRIPTION 

DADS is a mechanical computer- 
aided engineering (MCAE) software 
package that performs nonlinear large- 
displacement transient analysis and 
simulation. DADS is a cost-effective 
tool that allows you to model real- 
world behavior of complex mecha- 
nisms and systems before you build 
expensive prototypes. 

Unique features of DADS include: 

■ Feedback, control and hydraulics: 
DADS has a library of control and 
hydraulic components that enable you 
to model dynamics of feedback, control 
and hydraulic subsystems, their cou- 
pling with a mechanical system, and 
forces and torques that are fed back 
and act on components in the mechan- 
ical system. 

■ Flexible bodies: DADS combines 
flexible body analysis capabilities with 
traditional rigid body dynamics, using 
natural modes of vibration and static 
correction modes that are associated 
with kinematic constraint reaction 
forces. It includes an interface proces- 
sor for use with finite element struc- 
tural analysis codes to calculate 
stiffness, mass and deformation 
properties. 

Applications of DADS software are 
used in many industries, including the 
following: automotive, aerospace, air- 
craft, heavy equipment, materials han- 
dling equipment, military vehicles and 
robofcics. ^ 

I *Vy 1 / 

PROGRAM MODULES 

DADS consists of several modules 
that interact to perform an analysis. 
Commonly used modules include: 

■ PREPROCESSOR - The Preproces- 
sor allows you to define a model for 
analysis. A command-driven menu 
facilitates the input and editing of 
model data, while data entry errors are 
automatically screened. Extensive on- 
line help and documentation for begin- 
ners and an "expert" model of data 
entry make the program easy to use. 

■ ANALYZER - The Analyzer is the 
"heart" of DADS. A model is analyzed 
based on initial conditions and driving 
functions. Once a mechanism has been 
defined from the library of 2D or 3D 
model elements, the data set is pro- 
cessed and the mechanism is mathe- 
matically assembled and analyzed. 

DADS solves for joint position, dis- 
placement, velocity, acceleration, total 
and potential energy, and internal 
reaction forces of models. It performs 


assembly, static, kinematic, inverse 
dynamic and dynamic analysis, 
c Assembly analysis - determines 
whether all parts of a model can be 
successfully connected by the defined 
joints. 

g Static analysis - calculates the static 
equilibrium position and potential 
energy of the modeled system and 
components. 

□ Kinematic analysis - calculates the 
relative motions of bodies in a mecha- 
nism without regard for the mass 
effects of the components or any for- 
ces in the system. 

□ Inverse dynamic analysis - calculates 
the forces necessary to move a model 
through a predetermined path or 
given motion. 

□ Dynamic analysis - calculates the rel- 
ative motion of bodies in a mechanism, 
including mass effects of the compo- 
nents and any forces in the system. 

■ POSTPROCESSOR - The Postpro- 
cessor presents the results of an analy- 
sis in tabular or graphical form. Plots 
of any variable as a function of time, 
body position, joint forces, accelera- 
tion, velocity, total potential energy, 
and all internal reaction forces can be 
quickly examined in various formats. 


■ DADS GRAPHIC ENVIRON- 
MENT - This module allows anima- 
tion of analysis results. With simple 
assembly geometry commands, DADS 
Graphic Environment automatically 
links model geometry with position 
data from the Postprocessor and 
creates a realistic animation of a simu- 
lation. These animations allow you to 
quickly visualize the performance of a 
model, and provide an excellent tool 
for demonstrations and presentations. 
Options include wireframe and solid- 
shaded geometry with light sources 
and shadowing. 

INTERFACES 

DADS integrates with other MCAE 
applications such as FEA/FEM and 
CAD/CAE as part of a total design 
system. Interfaces are available for: 

■ AutoCAD • PDA/PATRAN ■ 
MSC/NASTRAN . COSMIC/NAS- 
TRAN ■ SDRC/I-DEAS ■ ANSYS. 

Other interfaces and custom inter- 
faces are also available. 

HARDWARE REQUIREMENTS 

DADS is available on all popular 
engineering workstations, main- 
frames, super minicomputers and PCs. 


STANDARD LIBRARY OF MODELING ELEMENTS 

DADS contains a standard library of modeling elements from which simula- 
tion models are constructed. 

Joints 

■ Revolute ■ Translational ■ Cylindrical ■ Spherical ■ Universal ■ Screw ■ Revolute-Revolute ■ 
Revolute-Translational ■ Revolute-Spherical ■ Revolute-Cylindrical ■ Spherical-Spherical ■ 
Bracket ■ Planar ■ Gear 

Forces 

■ Tire ■ Friction ■ Bushing ■ Leaf Spring ■ TSDA (Translational Spring-Damper-Actuator) ■ 
RSDA (Rotational Spring-Damper-Actuator) ■ User-Defined Force 

Constraints 

■ Fbint ■ Fbsition ■ Difference ■ Distance ■ Angle 

Controls 

■ Amplifier ■ Multiplier ■ Output ■ Dead Zone ■ Parameter ■ Delay ■ Sampled Delay ■ Sampled 
First Order ■ First Order ■ Sampled Second Order ■ Function ■ Second Order ■ General ■ 
Hysteresis ■ Summer ■ Input ■ Switch ■ Integrator ■ Limit ■ User-Defined Algebraic 

Hydraulics 

■ Accumulator ■ Hydraulic Motor ■ Check Valve ■ Double Actuator ■ Flow/Volume Update ■ 
Servo Value ■ Single Actuator ■ Valve 

Others 

■ Driver ■ Point of Interest ■ Initial Condition ■ Rigid Body ■ Curve Data ■ Header Text ■ 
Flexible Body ■ Impact ■ Reference Frame ■ Road ■ Track ■ Road Wheel 



COMPUTER AIDED DESIGN SOFTWARE. INC. 


P.Q Box 203 a Oakdale, IS 52319 □ (3l9) 337-8968 

Rev. 10/89 5M 
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CENTER AT A CROSSROADS 

After three years of operation, the Center for Simulation 
and Design Optimization has met many of its original three- 
year program goals, some with difficulty and others with 
a degree of success that has exceeded expectations. We 
have learned a number of lessons and have gained insights 
to guide the program in profitable directions for the future. 
Based on program restructuring agreed to at the Center 
Semi-Annual Program Review in April 1 990 and organiza- 
tional restructuring that has occurred during the summer 
of 1990, this newsletter (1) outlines Center progress in 
meeting original three-year goals, in terms of both soft- 
ware deliverables and basic technologies developed, (2) 
summarizes key lessons learned, (3) outlines a vision for 
the future of the Center, and (4) presents an organizational 
restructuring that dedicates management attention and re- 
sources to each of the major goals and provides a frame- 
work for efficient operation of the Center, 

Three-Year Progress: As indicated in the section that 
follows, most of the original three-year goals of the 
Center have been achieved. Spectacular advances have 
been made in real-time dynamic simulation, interactive 
dynamic simulation, operator-in-the-loop simulation, and 
continuum-based structural design sensitivity analysis. 
Significant progress has been made in computer network- 
ing, visualization, dynamic stress and life prediction, dy- 
namic system modeling, dynamic system linearization, 
Numerical integration of differential-algebraic equations 


(DAE), robotic system simulation, and track vehicle de- 
sign, Software integration for multidisciplinary simula- 
tion-based design has proven to be a more substantial chal- 
lenge then had been anticipated, particularly as regards 
communication of geometry data across different soft- 
ware platforms and creation of robust runstreams for 
multidisciplinary applications. While theoretical develop- 
ments in singular configuration and mechanism workspace 
analysis have been developed and published, creation of 
interactive tools for automated workspace analysis re- 
mains a difficult and illusive goal. Finally, extension of 
control and hydraulics capabilities in dynamic simulation 
proved to be overly ambitious undertakings, given the 
level of effort available in the Center program. 

Software deliverables, which are summarized in the fol- 
lowing section of the newsletter, are significant as indi- 
vidual capabilities. To realize their full potential, how- 
ever , they remain to be integrated into a broadly usable in- 
dustrial environment. For example, the new general-pur- 
pose dynamics simulation code, the dynamics workstation, 
visualization of dynamic system software, and the net- 
work computing system are now available individually. 
They now need to be integrated into a qualitatively new 
methodology for dynamic system simulation that can 
serve as the foundation for a next generation product to 
serve a broad cross section of US industry. 

Lessons Learned: The basic Center research program 
tenet that bal anoed- attention-shmild be placed on advancing 
the underlying theory and systematically implementing 
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methods in the form of usable software 
productive, it has been found tha evan , S proven ,0 be 
research software implementation A f ° r , expar 'mental 
planning is essential to accomplish d8S ' 9n and 

Object-oriented software desig! method ° d8d 9 ° a,s - 
mmg languages that were adopted m^ 3 ° d pr ° 9ram - 
three-year period of center operation ,he ,irst 

valuable constructs in numerous ! * proved ,0 be 

Panacea for resolving an challenne1 aS ’ bUt are not a 

offs involving the supportability !nd trada ' 

oriented design and the computational^ w ' ty ° f ° bject ’ 

crunching applications are ye, " 0 b . ° 1nUmb °' 


'ty across hardtJatl^ltfo^ms VaVS** S ° ,tware P° rta bil- 

than anticipated, provided tha, a i i n S i become ,ess difficult 

c °"»'«t, P or t^ilityof^appiic^^f?u l [ 0nment ’• used 
with geometric modelers has been t that must interface 
difficult than anticipated The lack of s ,an2 ,ar 
r,c m °deling makes it very difficult to on 3 990me ‘- 

metric information among app ica on, Unica,a 9a °' 
Finally, experience shows that / If pro 9rams. 
of methods developed in center research implemen,a tion 
m substantially different contexts an^ "’‘T b ® V ' ewed 
l.,.U of refinement and polish S"™ 11 at 
CAE software developed by the ™ h,erar chy of 

follows: <1),nd-us., ofi.nlladtoKwi!! "’“V vie » s0 •» 
the dynamics workstation and desion I£T*?-*®** stJch as 
workstation; (2) building block softi™ . T Y ana, y sia 
zation of dynamic systems, the new a8n « 3S v,suali ' 
namic simulation code and networ/ 9enara, - p urpose dy- 
(3) software frameworks sue” Ts tte Z'* 9 t0 °' S: and 

hfe prediction runstream that must bpL*^ 1 ® Stress and 
different combinations of CAE tools 9men,ed with 

plication environment. These and roi a » 8 !! h mdustr ial ap- 
outlined later in the newsletter have f. ®” 008 ,earn ®d 
agemen, in creatino a oraeSlf 9u,dad ca "‘er man- 

Center and in organizing the research oro" ,U,Ur9 ° f tha 
structive ends. 9 arch program towardcon- 


u «v.„, , y . using 

«"■“ f ">• llnssl d2S2 ITT* by ,h0 

Sim 8 ! y ,IVe years - until th « National* Ad US ,or a P prox ‘- 
Simulator (NADS) is constructed r Advancad Driving 

new facility and create a unique nation °i eXplol, ,his major 
factors research and mechamlT pr ° 9rani in human 
human operator, which is at the hearth" 1 d8S ' 9n ,0r ,ha 
neermg, opera, or-in-the-loon!!*./ co "current engi- 

— 


w -outers program 

While there is some 
tor-in -the-loop simulation goals' 1 of 9 ! neenn 9 and °P®ra- 
participants through the Center Adu danc 5 from Center 
past year suggests tha, it is i mnn J S ° ry 0oard over the 
function of the IDS and ultimateT/./ 18 /!.* 0 saparata the 
sic Center program in concurrent Inn •* NA0S ,rom ,he ba- 
us to restructure the manaoemlnf This hasled 
cated later in the newsletter * ° f ,be Can,ar - as indi- 

the basic thrusts and to better’manao V ' d8 Str ° n9 focus on 
to the Center. Professor Kyuno ChS h resourcas available 
sponsibility of Associate Director^ n ~ accapta d the re- 

neer.ng. Professor Jon Kuhlw.iul Conc urrent Engi- 

tor for Computing, and Professor j a as Assoc, ate Direc- 
as Associate Director for the lowL n S, ° ner wi " sarva 
restructuring has been carefullv/i° ° 9 Slmula, ° r - This 
merof iggo. with attention to mi/ 30089 durin 9 tha sum- 
ments for each of the Center's prai'/ft 3 ° d ,unc,ion s »ate- 
personnel. All concerned with 3 ° d mana 9 a ment 

are confident that the new struct!, 9e ^!° t ° f ,ha Can »er 
needs of participants and anh Ur8 Wl better serve the 
Center's program. "hance productivity of the 


MULTIBODY DYNAMICS 


PR°GRE ss in MEETING 
CENTER GOALS 


Vision for the Future • Th 

So " w,, ‘ 

Center, suggest that a refined vision for ! <ned b^ ‘he 
Center be established as the bas!i fo! .1 1" fU,Ure of tha 
and management. Two maior oonnrt 1 chnical Planning 

z*r.r~r.»« «*• «» 


coinplenienta,, b„, r~" «™ 

appear on the horizon for develoo/Jn! ^° ,Ca ' COn ‘ent, 
First, the emerging field of concurrent by tbe ^ an, er. 
the simulation-based design approach .® n9,nearin g. using 
gested by DoD industry study h3S been su 9* 

‘•cipants, represents an are 12 f h ° d ° Ur Ca "‘ a ' Pa!- 
a leadership role through computer aide^ C8 ° t8r Can play 
development and integration P to sudd!!, 809 !’ 88 " 09 * 001 
section of US industry The Center-!! 3 b r °ad cross 
strucuras. aad com™,. ^J.pfaa." '°°“ a »" 
adequately broad class of app|?ca,ion S ,n7 a " ab,e ' b ut 
tion for technology development and f °/ m the ,ou nda- 
‘ice of the profession in Ssl^Vry 8f l ° the prac ‘ 

x ce " ,e ' c °"- 

voh,c,a dn,,„ 9 Simula, io„^ C'ST ^ 
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Code(NGDC): The p- vers ion of slm ^latlo r] 

and has been released to several Dart ° C S availa b ,e now 
and examples manuals and an install / lcipan,s - with user's 
provides open- and closed-loo! !? 5 ?° 9U ' da - Tha NQ DC 
bihties. Five standard jo n ll , b ° dy ana ^ sis cap a- 
mants with linear and !o n i!ieTr 3 ° d tW ° forca e,a * 
actuator forces are provided Thr^'/u dampar - an d 
function generators fharmnn- ' brea different types of 
'<"*•) «. availab ,/ ,o a " d *P»"« 

A restarting capability is orn/i! °" llnaar ,or ce elements, 
“on. The NGDC provide! en!i!l ° f CPn,inuous simula- 
Position, velocity and a ! a o 9 8,s ,ha ,im e history of 

‘^ a coordinates' and each b^J 00 * e3Ch j ° iot ' in rala ' 

nates. Since the code is imolarr/’ ' n 9®neralized coordi- 
guage. the AT&T C ++ c SlTr 5 USi ° 9 ,he C++ lan ' 

Nr n!° 3re raquirad to install thi NG^C v U °' X opera,in 9 
NGDC will be available by the end rm 8rS '° n 10o “he 

version of the NGDC w,Vprov^e°d Ct ° b8 , r 1990 ‘ ™ s 

reaction force computation ca!!h!i!t r ^ 8r e,emants and 
documents, - ; ncta81n< rSPS t 'fnn y ^ om Piete software 

ficationf and SRS (S^ftwwL rL?*^* Parformanca Speci- 
will be available. Requirement Specification), 
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Dynamics Workstation (D\VS): .Version 10 of the 
DWS will be available by the end of October 1 990. The DWS 
is a modeling tool that assists engineers in defining me- 
chanical systems in a graphics-oriented, interactive way. 
Many useful modeling capabilities are included, such as 
topological layout, automatic assembly and loop closure, 
joint exercising, redundant joint checking, and initial 
condition definition. Geometric data that are needed to run 
VDS and modeling data that are needed to run the recursive 
dynamics code and DAOS have been produced. Using the 
DWS, journeyman engineers are able to easily create 
mechanical system models and evaluate their perform- 
ance, using advanced simulation and animation tools. Four 
technical reports and one software document have been 
written for the DWS. 

Basic Technologies 

Recursive Dynamics: The recursive formulation 

without acceleration elimination (Order N 2 ) has been em- 
ployed to develop the NGDC. An object-oriented software 
design has been developed to construct a software struc- 
ture that provides for extendibility. After extensive 
review of the {J -version code, an inter-object data access 
method has been created to improve efficiency. A tree- 
traversing method has been developed for closed-loop 
analysis and for system mass and force vector computa- 
tion. For numerical solution, alternate types of linear 
equation solvers will be employed to match the character- 
istics of matrices used in analysis. Recursive dynamics 
formulations, both with and without joint relative coordi- 
nate acceleration elimination, have been developed for 
rigid and flexible multibody dynamics. Extensive analysis 
and testing of alternate formulations have provided a 
thorough understanding of trade-offs among formulations, 
for generality and computational efficiency. Extensions in 
conceptual formulation of flexible multibody dynamics that 
define system topology in terms of frames and transfor- 
mations, rather than bodies and joints, have demonstrated 
potential for more than one order of magnitude speed-up 
over conventional formulations for flexible system dy- 
namics, to support broad classes of simulation and life pre- 
diction. These developments form the foundation for the 
next decade of multibody dynamic simulation, in both 
concurrent engineering and real-time operator-in-the- 
loop applications. 

Differential-Algebraic Equation (DAE) Solvers: 

Fundamental developments in methods for implicit numeri- 
cal integration of DAEs have been completed and imple- 
mented, with good computational efficiency and excellent 
numerical stability. Generalized coordinate partitioning 
and tangent space parameterizations of constraint mani- 
folds have been shown to be practical and computationally 
effective. Theoretical developments in error analysis 
have shown that the algorithms are rigorous, stable, and 
computationally effective. Algorithms based on these de- 
velopments are currently being implemented in the Carte- 
sian coordinate formulation by CADSI for the DADS code 
and in numerous applications within the Center research 
program. 

Workspace Analysis: Theoretical foundations for 

mechanism singular configuration and workspace analysis 
have been developed and demonstrated. Singular configu- 


ration manifold mapping methods developed by the Unive 
sity of Pittsburgh have been adapted and successfully use 
in Center research. While the feasibility of singul. 
configuration and workspace analysis has been estal 
lished, resources have not been available to implement 
general-purpose software capability for industrial appi 
cation. Continued thesis research is pursuing formulatior 
that are candidates for implementation with the ne. 
generation dynamics code. 

Validation Library: To carry out validation fr 

current flexible multibody simulation tools such as DAD! 
DISCOS, Order N DISCOS, and TREETOPS, a verificatic 
library system is being developed to store simulatio 
models in a standard format, translate the standard form, 
of the simulation model description into various simulatic 
inputs, launch various simulation tools, and manage simi 
lation and experimental data. A standard input dat 
definition has been adopted. A formulation for translatin 
standard input data into input data for DADS, DISCOS 
Order N DISCOS, and TREETOPS has also been developer 
Topology analyses of the mechanical system, recursiv 
position and velocity analysis, and initial assembly base 
on constraint error minimization are making use of trans 
lation of the standard data. Since most flexible dynamic 
simulation codes are based on the deformation mode ap 
proach. generating standard data for a flexible bod 
requires processing finite element output data to dynamic 
input data. 

OPERATOR-IN- THE- LOOP SIMULATION 
Deliverables 

Iowa Driving Simulator (IDS): With the aid of Dol 
and NASA participants in the Center, major graphics am 
motion base assets have been obtained to construct the IDS 
The University has provided $1.5 million to support com 
pletion of graphics equipment and construction of a build 
ing addition to house the IDS, and integration of the system 
The IDS will represent the finest driving simulator in thi 
US and will serve a broad range of Center participants an< 
other industrial and government organizations. An exten 
sive technical report that defines IDS capabilities and e 
business plan for sustained operation have been completer 
and distributed to participants and other potential users o 
this major facility. 

Basic Technologies 

Real-Time Dynamic Simulation: The feasibility o 
real-time dynamic simulation, through parallel processing 
on the Alliant FX/8, was demonstrated during the firs'i 
year of Center operation. Specialized parallel processing 
tools were developed to facilitate the development ol 
parallel dynamics applications. These tools can organize 
and manage parallel execution threads with very high ef- 
ficiency. The parallel processing tools expedite develop- 
ment of parallel applications and provide easy portability 
among various parallel computing platforms. These tools 
are currently being used in development of the new genera- 
tion parallel dynamics codes. 

During the second and third years of Center operation, the 
real-time dynamics focus shifted toward achieving real- 
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dynamic stress and life prediction 
Software Deliverables 
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Basic Technologies 
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Software Deliverables 
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database system has been developed to support the design 
process. The user interface and computational modules of 
the DSOW that support the user interface have been 
implemented in C++. Currently, the geometric modeler 
that is supported by the DSOW is PATRAN. The structural 
analysis codes ANSYS, MSC/NASTRAN, and ABAQUS are 
supported. For remote process execution and distribution 
of DSA computation, any Unix-operating machine can be 
used. The portability of the prototype DSOW has been 
verified on an Intergraph CAD workstation. Fourteen 
technical reports and three software documents have been 
written for the DSOW. 

Basic Technologies 

Shape Design Sensitivity and What-if 
Workstation: A methodology for integrating technolo- 
gies for planar and spatial structural component shape 
design has been developed. The feasibility of the technol- 
ogy has been demonstrated in a prototype Shape Design 
Sensitivity Analysis and What-if Analysis Workstation, 
which has been developed using Apollo Dialogue and FOR- 
TRAN. Capabilities developed in the Shape Design Worksta- 
tion include finite element error analysis and mesh adap- 
tation, design parameterization for curve and surface 
boundaries, automated boundary and domain velocity 
computations, automated design sensitivity computation 
using both direct and adjoint methods, visualization of 
design information, and automated geometric shape and 
finite element mesh updates that result from design changes. 

Configuration Design Sensitivity Analysis: For 
built-up structures, a new DSA theory for configuration 
design variables has been developed, to allow layout design 
of built-up structures that include truss, beam, mem- 
brane, and shell design components. Once implemented, it 
will provide a unique capability that can support layout 
design of built-up structures such as vehicle, aircraft, and 
space structures. 


DSA of Dynamic Frequency Response: Acontinuum 
design sensitivity analysis method for dynamic frequency 
response of structural systems has been developed. A 
variational approach for non-self adjoint operators, using 
complex variables, has been developed to obtain design 
sensitivities of frequency response with respect both 
sizing and shape design variables. This method has been 
tested using a large vehicle system finite element model 
that was developed by Ford Motor Company for NVH (noise, 
vibration, and harshness) design optimization. This new 
capability is being implemented under support from Ford 
Motor Company, so that it can function with a Ford 
proprietary NVH analysis code. 

DSA of Nonlinear Structural Systems with Criti- 
cal Loads: A new continuum formulation for DSA of 
critical loads with respect to sizing and shape design 
variables has been developed for nonlinear structural 
systems with geometric and material nonlinearities that 
are subjected to conservative loading. This new DSA ca- 
pability can be used to obtain an optimized vehicle body 
design that satisfies safety requirements for vehicles 
under crush conditions. 

DSA of Transient Structural Dynamic Response: 

A new continuum-based DSA theory has been developed for 
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transient dynamic response of built-up structures. Th 
DSA method does not require derivatives of eigenvector 
and Ritz vectors. Numerical experimentation with th 
method indicates that it is extremely efficient and accu 
rate. 


NETWORKING, VISUALIZATION 
COMPUTATION (NVPC) 


, AND PARALLEI 


Software Deliverables 


Visualization of Dynamic Systems (VDS): VDS 
was developed by e Center during its first year of opera 
tion. Since that tim^ has been extensively enhanced. VD5 
is currently supported on a wide range of graphics plat 
forms. Capabilities of VDS continue to be upgrade- tc 
enhance performance and extend functionality. An IGES 
translator that supports a number of entitles, including 
lines, planes, parametric surfaces, and finite element rep 
resentations, is currently undergoing internal testing anc 
will be included with the next release of VDS. Work is alsc 
continuing on a Motif-based graphical user interface foi 
VDS and on a number of enhancements to performance anc 
networking capabilities. 

Network Computing System (NCS): The NVPC 

project prove s the NCS - which is used by VDS - 
other Center deliverables to participants. The Center ; 
ported NCS to a variety of platforms and currently can 
provide distribution and support for seven different plat- 
forms. The NVPC project also furnishes internal technical 
support for other l/UCRC project areas that use these 
tools, X-windows, Motif, and parallel processing tools. 

Basic Technologies 

Para I lei Processing Tools: Specialized parallel 

processing tools have been developed to support the needs 
of high-speed dynamic simulation. Fully compatible ver- 
sions of these tools are currently supported on the Cen- 
ter's Alliant FX/8 and HP/Apoilo DN10000 multiproces- 
sors. These tools are currently being used by the m 
body dynamics and operator-in-the-loop simulation p |« 
ects to develop high-speed and real-time dynamic simula- 
tion codes. 


Distributed Computing and User Interfaces: Center 
staff are actively tracking developments in the Open 
Systems Foundation (OSF), particularly with respect to 
the selection of standards for distributed computing envi- 
ronments and user interfaces. Tools currently in use at the 
Center are compatible with selections made thus far by 
OSF . The Center will continue to assure that computing 
tools used in Center projects remain compatible with pre- 
vailing OSF standards. 

TOOL INTEGRATION FOR CONCURRENT 
ENGINEERING (TICE) 

Software Deliverables 

Integrated Analysis Capability (IAC): The Center 
has been successful in usinq an engineering database and 
module TTTHnagmTM»iii sysiern, called the integrated Analy- 
sis Capability (IAC), from Boeing. IAC is used in the 
Center's Dynamic Stress and Life Prediction project for 
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software Deliverables 
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yet modular, code that will provide detailed time and space 
dynamic tire loads. The pavement models use finite 
element and plate models that are capable of analyzing 
jointed PCC concrete pavements of varying thickness and 
slab length. 

LESSONS LEARNED 

In the development and implementation of simulation and 
design optimization software, engineers and programmers 
in the Center have learned that: 

(1 )Software planning is crucial for software develop- 
ment, even for research. 

(2) Object-oriented software design is valuable, but 
not a panacea. 

(3) Hardware platform portability is easier than an- 
ticipated, as long as Unix operating systems are 
used for all machines. 

(4) It is easier than anticipated to develop a unified 
interface to finite element analysis codes such as 
ANSYS, NASTRAN, and ABAQUS. 

(5) It is very difficult to develop a unified interface to 
geometric modelers, because there are no stan- 
dards. 

(6) Different categories of software deliverables have 
very different characteristics and levels of re- 
finement, as follows: (a) broadly usable and re- 
fined end user software, e.g., DSA workstation; 
(b) building block software, e.g., VDS and NGDC; 
and (c) software frameworks, e.g., IAC and dy- 
namic stress and life prediction software. 

To illustrate some of these points, consider the OSOW, 
which is a large software system that can be character- 
ized by the following attributes: 

(1 ) Multi-platform Support: As mentioned in the July 
newsletter, due to requests from several industry par- 
ticipants, the Design Sensitivity Analysis and Optimization 
project switched its focus in the DSOW from capability ex- 
pansion to software porting across many different hard- - 
ware platforms. This is complicated because the required 
software packages may be distributed across several 
machines. Currently, we are targeting Apollo, SUN, 
Intergraph, and DEC workstations with Unix operating 
systems. Some of our participants have distributed access 
to finite element software packages such as ANSYS, 
NASTRAN, and ABAQUS, thereby requiring distributed 
computation on machines such as Cray and Alliant main- 
frames with Unix operating systems. 

(2) Large and Complex Software: The DSOW has a large 
database, sophisticated user interfaces, a large number of 
modules, sophisticated runstreams, many languages (e.g., 
FORTRAN, C, C++), and interfaces to many external 
software packages such as PATRAN, ANSYS, NASTRAN, 
and ABAQUS. Development of a unified interface to finite 
element analysis codes is easier than initially anticipated. 
However, development of a unified interface to geometric 
modelers is very difficult, because there are no standards 
in geometric modelers as there are in finite element 
analysis codes. Due to non-existence of standards in 
geometric modelers, it is necessary to design a different 
interface for each geometric modeler. Initially, the proj- 
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ect intended to develop the DSOW to work with PATRAN ar 
IDEAS modelers. However, major resources are require 
to design and implement a new interface between DSOW ar 
IDEAS; e.g., TACOM and DICE are supporting substanti 
DSA and Optimization project efforts to design and imp! 
ment a new interface between the DSOW and the Intergrar 
l/EMS modeler. 

VISION FOR THE FUTURE 

Technical progress in Center research and interaction wii 
both government and industrial organizations that ar 
pursuing significant national goals suggest that there ar 
two complementary but quite technically different area 
of future contribution from Center research. 

Concurrent Engineering: First, the emerging nation; 
thrust in concurrent engineering can profit from integn 
tion of computer-aided engineering tools that are being d i 
veloped and implemented by the Center. The vision fc 
evolution of concurrent engineering presented by the Dot 
Industry Working Group on Concurrent Engineering of 
chanicai Systems proposed that the evolution of concui 
rent engineering tools should first embody simulation 
based design, followed by evolution to a design optimize 
tion methodology. This conclusion, which was based in pai 
on interaction of the DoD/Industry Working Group with ou 
Center activity suggests that the Center can play a signifi 
cant role in development and implementation of concurrer 
engineering tools for mechanical system design. The rol 
and scope of our Center in this major development, how 
ever, must be carefully considered and designed in prac 
tical terms. 

The major DARPA Initiative in Concurrent Engineerint 
(DICE) is creating a framework for concurrent engineer*™ 
of broad classes of engineering systems and is including at 
emphasis on manufacturing. Our Center program, ii 
contrast, has focused on engineering tool integration tt 
support concurrent engineering of mechanical systems 
Experience gained to date in implementation of dynamii 
stress and life prediction software and in the TACOM pilo 
project on tracked vehicle design suggests that the scop* 
of our Center activity is adequately broad to make a majo 
contribution to the field of concurrent engineering o 
mechanical systems, it is also technically coherent enougf 
to permit the Center to make significant and practice 
contributions to the field. Joint projects being undertake* 
with the DICE program will provide the Center an opportu 
nity to contribute tools to a framework that our partici 
pants may ultimately expeetto use in concurrent engineer 

ing of their products. Developing such partnerships fur 
ther permits our Center to concentrate on its strengths 
work within its physical and manpower limitations, am 
take advantage of capabilities that are being developed ii 
large programs that are beyond the scope of our Center. 

Apart from the real-time operator- in- the-loop Simula 
tion emphasis in the Center program, virtually all of the 
basic Center program falls well within the broad goals o 
concurrent engineering. With minimal project restructur 
ing-agreeefio-at the April 1990 Semi-annual Prograrr 
Review, the Center program can be continued and strength 
ened to develop and transfer concurrent engineering desigr 
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approximately 56 million worth ol advanced o^„°h 
motion-base equipment and has allocated 3nd 

construction of a facility to house the i I,l0n for 

lator (IDS). The facility is now^nV '° Wa Drivin 9 Simu- 
be complete in February 1991 Thte C m a StrUC,IOn and wil1 
facility will b. th. host groonVyalmf. rt? IOr " eW rasaarel ’ 
the US from its completion in 100? * J" 9 s,mulat °rin 

NADS becomes the world's best advance dd? 6 ’ Whe " th ® 
tor. The IDS, with minor anhlll a<3vanc ®d driving simula- 

be comparable in capability to^he currem^t 9 ^ 88 ' Wi " 
art driving simulator that is operated bv Daim| atG D ° Mhe * 
Berlin. The University of Iowa ta ^ rnier-Benz in 

relations with the Daimler-Benz res^arcrim 11 ^'" 9 i0 ' nt 
a joint government/industry group in France th P » and Wlth 
to build an advanced driving simulator m hat ln,ends 
four years. The University of Iowa simulatn?t Pr ° Xlmately 
has been spawned by our Center is olavino e team ’ Wh,Ch 
evolving the new capability for vehiri. I 9 ma,0r ro,a in 
both in the US ^^elsewhera^n th^ simulation, 


tha^mam agamer!™ of* t h e ^ agreed 

should be separated from'*?.* 23? <?„£? Simu,a,ors 
concurrent engineering, for purposed of eff!^? ram m 
dedicated management The man effective and 

presented in thenext section T,. agemant ^structuring 
designed to dedicate aSh °o bothm"'*’*"*" has baa " 
neering and operator-in-the i 0ft „ b * h ^ concur '’ent engi- 
ne Center, maintaining a sTro nn > f S ' mula,ion objectives 0 f 
rent engineering toot develoDmant US °i n Can,0r c °ncur- 
Center participants with Ten?, To, h°. L’"* b ' b «'"3 
loop simulation activity. f the °Perator-in-the- 

CENTER ORGANIZATION 

and to improve th^ manlg?menTof P pro f , ft"**!. aC,ivities 
management structure has been re^?' The " Cen,er 
the new organizational structure is to 1*]? od) ® ct,v ® °f 

and main, 'a,?." 5 *'« id.ntillad 

::: .ss U .s 

among projects. 110 ' 0553 '^ ° verla ' > •' duphcatipn 
(5) Common comouter *ri* n *~ 

and provided for multiple projTST '* iden,i,ied 


< CC P AD > *- bB ®n 

l/UCRC. (2) Special ’Si LI ? 9 * V (1),h ® 

Support, and (4) the Iowa Drivinn e- C , ompu,ln 9 /Sof tware 
"«» das its iJ„ EaChbro »'>'" 

associated proaram anri 9 e,’ . e ^ ,or - who manages the 
dv.„ll managamant ^ f ,S ' h * 0ir,c ">' 01 =CAD in 
ministrativa and technical ssoclaIe Directors have ad- 
gram. P,ol."o, ' Ky„ «•’ pro- 

hay. accoptad r.sponsibihtla . ‘ 1 K “ h '' ““ Jim S “"«' 
Concurrent Engineerina Cnmr. *• Associale Directors for 
Simulator. Ed Haug will serva PUt i ng .’ and ,he ,owa Oriving 
Projects, an, it Si” Ac,,n ? Di ' 8 d"» Spscial 
l/UCRC, there are einht r,r P ® SOn can ba h,re d. Under the 

Haug): Oparator-innhe^Loop' Sim“!' i , b0,iV D, ' nam,es (Ed 

OSA and Optimization ,Ky° n g Ihoh Tmo", ! " 

Conc urrent Engineering (Kirk W ! ’ r ,nte 9 ra tion for 

Controls and Tele-Operat on /r^" 8 "* En 9 ineerin 9 and 
UCRC that were aoreaH new P ro f acts in the 1/ 

Review. ^ ,0 at the April 1990 Program 


^ ... mo-iuup simulation, the Cantor ^ - 

mented interactive simulations of telerohn., ' mple ' 
with the NASA/Goddard Space Flint!# ‘ , bo,lc systems 

("9 signilicant i 0 lrtt »rt^S^,lLS; n ': r * ndl **»P. 0 '- 

Uni,.,si,yolCalilo,nZo.W, in hmh a,C 9 '°“ p a ' Tl, » 

imn With high-spaed mallinroJas?n h , . ay ' epair “““a- 
in the Center laboratory and^he m a SUperworks *ations 
motion systems that are balnn i„7 ' ' VISUalizati °n and 
Driving Simulator, the Center and Un Sty h,! J h# ‘° Wa 
opportunity to play a leadership role in a broad qU ® 

pMcat t ions an ' PU * a,0r deSign and buma n facto^s^esearch^pl 

^n^r^“^^^^ a ;- p opsim u . 

scpport th. Iowa Oriying Simula, andmiTiS^! Ad" 
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SS53SSS uada ' 

Optimization, DICE Project S j| r ° S ’ o° rd NVH Desi 9n 

Suspension Design Project^To condflvt Proj8ct ' and »TC 
ware and computer seiin? a J consolidate computer hard- 
several projects the Cn su PP or ^ af 0 necessary for 
gram haIw^ a .5-°7S: 9 ' Sp, '» a ' a S ppp °rt Pro 
projects^ Graphics Facility r p ogram ' ,h « r s are three 
p artic,pan, ScSw ^.^^“rt^ 




F— 1 1 






















“«• ol cVmirTJsouTaT h Th e re |o Ulte n ia m0r8 efficiant 

Program has been developed to ma** 0nvin9 Simulator 
recently acquired Iowa Driving S i m u lalo ^ f5s 1° T^ 3 * 0 ^ 
lour projects in this proqranr i ma „ / ! DS, ‘ Ther ®are 
Database Management! ExperimelTaM? Genera,ion and 

and Systems. ana v . h|e ° 

«■-- 

Mission and Function Statements El has ' developed 
principal working relationships with other 38 d ® fined 
cipal working relationships identify com^T**’ Pri "' 
science support for each project and c °mputer 

that the projects provide P and receiv^to ' Yrh" 9 SUpport 
missions. The Center will n B „ ? 1 achieve their 

projects and Associate Directors ThYYl*? MFS ,or al1 

will be used for internal program LT ST 00 ^ 8 "' 
projects accomplish their missions aSSUr ® ,hat 





ls sues in the design and imnUmJ 3 d that tho im Portant 
"isms are different TrLSZT ° f ,h ® S ® macha ' 

support of loop level-parallelism An e V ' OUSly S,Udied ,or 

actual application is presented m . XT ’* 1 Uwl r et u 
the specific nature of run-tima 8 ,Udy mdlcat ®s that 

dramatically impact the performs^ 0 ^ ^ 8chanism s can 
gram. The study suooests tha, ° f a parall ®< pro- 

9'ven to the choice of an approphale sit a of en,l ° n mUSt b ® 
agement mechanisms to suoDort c 8, .° run ' ,lm ® man- 
•xplicil heterogeneous psrshelrsm Th!'" 3< ”’ a,a,ad " 
syntactically closed constructs '^suggested* S,mPl •• 

'NTEGRATO ccJKE 2 E 32 K& 
robotic systems u 

L "c:?,“rr" d u 

Center Technical Report R-69 


REAL ' T 21! DYNAMIC SIMULATION OF A RORnr 
system with rotor oynami^ OBOT 

U-Ping Chuang ano Sung-Soo Kim 
Center Technical Report R-67 

Rotor dynamics are important for robots that =, ^ • 

by motors and high aear ratio non, 00,s that ar ® dr| ven 
generic recurs!.. dynamic to ™ SS,dn *>*•"•■ Net, 
■ems with motor rote °' ,ob ° li ' : ***■ 
motors can be mounted on inboard P r ®sented, where 
actuated joints. An organic " du bdard todl « « the 
ms equations is obtained through adopttonof'a «« d "’ am ' 
representation in equation derivation ru a,e v ® c,or 
rotors on the overall system dvna' Th ® m ! ,uenc ® of 

Parallelism of the derived dvnam; y m,cs ,s identified. 

exploited and,, eahZ. oarall.^ '° ,mula,id " a » then 

is P-oposed. A seven d. 9 ?.. o a , C .°7 U 'f , I,nal al9 ° ri,hm 

to illustrate rotor dynamics elfects ,a, atobot is used 

model. A computational speed that is lasher m mP “' a "°" al 
simulation is demonstrated fhan real-time 

MANAGEMENT OF HETEROGENEOUS 
PARALLELISM ON SHARED MEMORY 

multiprocessors 

Athar B. Tayyab and Jon G. Kuhl 
Canter Technical Report R-68 

This report considers the problem oi men, 
heterogeneous parallelism on shared " r'*'" 01 
processing systems. Heterogeneous oareufi ' pa,allel 

^r^ ou ^^: pp s pa -'Sii,m^^s:d 

compiled g.n.rrdrmt d :,Tar.?, y ha C , 0dM a " P 


support 
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speed^Tmu?^ ,or ,h ® high- 

cal systems or independent control con,rol, ® d mechani- 
and implemented. A general control^*.* 8 " 18 ' S d#v ® lo P*d 
method is provided in order to no r yst ® m Preprocessing 
topological analysts andtoaSt. "T, th8 C ° ntro1 •*«*•<" 
computational algorithm The can^ V 88< ° P ,h * para,,el 
system and mechanical svsfnm °‘ 8ys,om ’ hydraulic 
analysis. A parallel comiutLon^i Jf* ln ‘ 99ra, ® d *n the 
tor the integral control/hll! /■' a ' 90r ' thm is developed 

analysis. A ■»«.-» 

so that the complete controK COd ® 18 9® nara ted 
tern can be simulated as a sinot » 8Chanical s y s ' 
recursive computational methn^t ' ntegrat ® d m °del. A 
robotic system analysis is alio f for In, ®9rated controlled 
formulations for robotic tranJfar* °?* d ’ New rec ursive 
dynamics of robots with rotnp ° rmatl0ns a dd for the 
analysis capability in robotic ® ff ® cts .P rovid ® a general 
and inverse kinematics are former ' a * ,ons ' Th ® forward 
approach. The so-called rob^iJ ’ u 9 a r ® cur sive 
Jacobian are also derived Two « Jacob,an and inverse 
,c * squations of robots with rotor « ° recursiv ® d ynam- 

th ® set is for rob^ts^wjtlT^otrire 18 3r8 ^ orm °f a **d : 
outboard bodies of the actfnn »«• niounted on the 

for robots with rotors mounted n^th’ th8 88COnd set is 
acting joints. Equations for a tvofcal n?° 3fd b ° di ® 8 of ,h ® 
a transmission model are atso^l °. C m ° ,0r mod ®' and 
planning method in Cartesian ,ormu lated. A trajectory 
-halhbd generates a aS N n^*, ' S deVa ' 0pad ' Th * 
specified points in Cartesianll betw «®n any two user- 
frajectory between fwo neiohbol 8 3 Smooth ,ran si tion 
ferent position and hybrid position/? ® S ® 9ments - Dil ' 
are studied and tested A resoivln f C ° ntr0 ‘ m ® ,hods 
silion/lorc. conlro, m«to* d " T, se d ™ T”" h,b,i “ pd ' 
verse kinematics and dynamic, r „ . h recursiv ® m- 

and implemented. Jointspace comT ^ atl0ns ’ ,s d ®signed 
ces are derived and usedln th« h ? ' 3nC8 selection matri- 
acceleration hybrid controller d ,vatl0n of this ^solved 
algorithm for tba iotegral_an a i^ P arallel com Putational 
underp^ioir/^on,^!^ 18 °! robo,ic s yst®ms 

Fine g,a.n parallelism is ful y exploitT/in !u d Imp, ® m ®nted. 

12 " ex P'oited in this method. The 
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efficiency of this algorithm is demonstrated through simu- 
lations on the Alliant FX/8 multiprocessor. 


INTEGRATION OF CONTROLLED MULTIBODY 
MECHANICAL SYSTEMS WITH NON-STIFF 
MECHANICAL SUBSYSTEMS AND 
STIFF CONTROL SUBSYSTEMS 

Shih-tin Lin and George M. Lance 
Center Technical Report R-70 

The simulation of controlled multibody mechanical sys- 
tems, with nonstiff mechanical subsystems and stiff con- 
trol subsystems, usually suffer from the long run times. 
In this report, a hybrid integration scheme that uses 
separate algorithms for the mechanical and control sub- 
systems is presented. A double pendulum driven by two AC 
motors and a backhoe mechanism driven by three hydraulic 
actuator were used to demonstrate the computational 
efficiency of the proposed integration scheme. 


COMPUTATIONAL LIFE PREDICTION METHODOLOGY 
FOR MECHANICAL SYSTEMS USING DYNAMIC SIMULA- 
TION, FINITE ELEMENT ANALYSIS, 

AND FATIGUE LIFE PREDICTION METHODS 

Woon K. Baek and Ralph I. Stephens 
Center Technical Report R-71 

This study deals with an integrated life-prediction rr n- 
odology using dynamic simulation, finite element analysis, 
the fatigue life-prediction method, and experimental vali- 
dation for the finite life evaluation of a mechanical system. 
As a practical example, a multi-body dynamic model of an 
existing ground vehicle was developed using kinematic 
joints and components. The vehicle model was hypotheti- 
cally run over a measured road profile at constant speed. 
From this dynamic analysis, load histories were obtained 
for each component. A lower control arm, which is a 
critical suspension component, was selected for the co- 
mopnent fatigue life prediction, and a finite element model 
of this component was developed. Several high-stress 
regions were identified from the finite element- stress 
analysis. Local notch stresses”at each high-stress region 
were then obtained at potential fatigue crack ‘initiation" 
points, which can be called potential fatigue-critical loca- 
tions. Dynamic stress histories at potential fatigue-criti- 
cal locations were produced by the quasi-static approach. 
The local strain-life method was used to predict fatigue life 
of each potential fatigue-critical location. The fatigue life 
was defined as the typical crack "initiation" life of a crack 
about 2 mm in length. The fatigue life of the component was 
defined as the shortestfatigue life among several potential 
fatigue-critical locations. To validate this computerized 
■>rocedure, the lower control arm was experimentally 
.ested for stress and fatigue durability. The brittle coating 
method was used to identify high-stress regions. Experi- 
mental stress analysis was carried out using strain gauges. 

F- 


The experimental results and predicted results based upc 
finite element analysis were very close: every compai 
son showed a difference of less than 5 percent. Als< 
fatigue durability tests were done for the component t 
repeatedly applying the same load history that was applic 
to the finite element model, until a 2 mm long fatigue crac 
formed. The breakthrough in this integration is th; 
dynamic stress and fatigue life can be predicted in ear 
des stages without experimental measurement of eitht 
loads it stresses. This methodology can be extended 1 
design optimization, based upon durability. 


REAL-TIME OPERATOR-IN-THE-LOOP 
SIMULATION OF MULTIBODY SYSTEMS 

Joe Lan Chang, Sang-Sup Kim, 
and Edward J. Haug 

Center Technical Report R-72 

This thesis presents a general approach to achieving real 
time operator-in-the-loop simulation for multibody dy 
namic systems. Emerging real-time dynamic simuiatioi 
methods are used to demonstrate the potential for creatim 
interactive design workstations and for teleoperatim 
space robots, with a human operator in the control loop 
The recursive formulation of multibody system dynamic 
with relative coordinates is employed for efficient nu 
merical analysis and implementation on a parallel com 
puter. High-speed computer graphic techniques are evn 
ployed to create realistic visual cues for the simulator. 
A simulator is developed in this research by integrating thi 
real-time dynamics program, a realistic graphics display 
and the operator's control interface. Real-time operator 
in-the-loop simulation is analyzed, as regards the goal o 
real clock time, not only with respect to dynamic Simula 
tion but also with respect to graphics display and the 
operator interface. Synchronization of the simulation ii 
found to be most important for realism of the simulator. 
A backhoe simulation is implemented to demo nstr ate the 
capability for operator-in-the-loop simulation. The simu 
lator is developed by modeling backhoe dynamics anc 
hydraulic systems with the recursive formulation tc 
achieve real-time simulation, developing an interactive 
graphics program for visual cues, and interfacing the op 
erator’s control action with the dynamic simulation througt 
a pair of joysticks. The simulator is also implemented fo 
teleoperated simulation of space robots. 


CONTINUUM DESIGN SENSITIVITY 
ANALYSIS OF STRUCTURAL DYNAMIC 
RESPONSE USING RrTZ SEQUENCE 

Kyung K. Choi and Semyung Wang 
Center Technical Report R-73 

In this report, a unified continuum-based sizing Design 
Sensitivity Anal ysis (DSA) method is developed for the 
transient- dynamte “response of built-up structures by 
taking design derivatives of the variational equation of the 
built-up structure to obtain a variational equation of the 



The direct 

built-up structure is used to obtain an aDn?i eqUa “° n ° f ,he 
of the variational equation of the rie pprox| mate solution 
large size built-up structures the 9 ° S8nsi,ivit V- For 
method that is used to reduce then S3m0 super P°S't<on 
aquation of the bull up s^c u r s 1 ™"*'™ ° f ,he 
dimens, on of the matrix e S L f8duc ® tha 

For accuracy of analysis results and effime 9 ° S8ns,tivit y- 
tions of eigenvectors and Ritr v* ®ncy, combina- 

can ba implemented outside established Ffa^a method 
postprocessing data only s/nro *h« d f codes using 
quire derivatives „? J” ^ V method p oes no, r e 9 
matrices. Moreover the method it' Jr ‘” P "’ 3 ' a " d mass 
not require rterivatiei. o™ at, “c/om* T ^ ^ 
are presented to demonstrate * ; Two exam ples 

amples treated in this paper indicate^a?!^® me,hod - Ex ‘ 
of basis vectors that are used for a f th ® sam ® number 
-cure is Ohouqh ,0, IXtXl Z 


linearization in the 

RECURSIVE DYNAMICS FORMULATION 

Tsung-Chieh Lin and K. Harold Yae 
Center Technical Report R-74 

pose a d if He u 1 1 ^p r o b^e m ° In 71 linear 0 multi '[ )od '' dynamics 
therefore deslrabl. to hav , TnlJrZ d0Si9 "- “ is 
conjunction with a general -™ *",?■?. capabili, y in 
modeiing technique A new commit**™" 1 tl !* 0(i Y dynamics 
aarization is obtained by apoS f °°f ' m8,hod ,or lin ' 
analytical approximations to the fecuiIhrV ki flrS, ‘° rder 
lationships. The method has proved to b*L™ T** r8 ‘ 
more efficient. It has also turned out to 8 * computationally 
because the analytical Derturhat; ' ° b ® m ° re acc ura,e 
vector operations by circumvent no r ® qU,r8s matrix a "d 
tiation and other associated mim • umer,ca ^ dWeren- 
accumulate computational error opera,i ons that may 


RECURSIVE LINEARIZATION OF MULTIBODV 

dynamics and application 

TO CONTROL DESIGN 

Tsung-Chieh Lin and K. Harold Yae 
Center Technical Report R-75 

The non-linear equations of motion in multi b n H u w 

pose a difficult problem in ii„ a " mum-body dynamics 
solution is ,0 haveThe non linla T r °' desi 9 n ' One 

« dd -SitSTSS '" ,ui - 

from a non-linear model by omittino non , 8X,rac,ion 
such as Coriolis force, centrifugal force andoth f ^ ,,8Cfs 
interacting between bodies S.w, ° 0,her ,orces 

anadhocapproacrdependmo oo ^ >,mp, ' f ! cation >• rather 
therefore desirable to have ZarLZ ^ ^ " is 

i-«Oh »„h a g .„^~ a ^^h S 
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modeling technique. The linearization * 
research starts with non-linoar 1 dev ®'oped in this 
ten in the Newton-Euler form whirh^' 005 0< rT10, ' on writ- 
then, they are transformed into joint coo!?** *° construct : 
tation in two steps* Cartesian ua t-i rdma,e represen- 
ables, and state variahilt . ,ables into state vari- 
transformation the kinematic rela,'"* Variables - At ® a ch 
their linearized relations A newcomn S » V* replaced with 
linearization is obtained bv aonlvin P at '° nal me,hod for 
analytical approximation J X 9 3 S0r ' 8S 0t border 
lationships. The method has proved^b7co kln ® ma,iC f8 ' 
more efficient. It has also turned ou B 30 COmpu,ationai| y 
because the analytical oerturha.i * 1 b fT10ra acc urate 
vector operations by ci?cumvSminn r0quir8s ma,rix and 
tiation and other associated numor umerica ^ differen- 
accumulate computational error Tha 9 ®' 3110 " 8 ,hat m a f 
posed linearization alaorithm i« ^ 8 power °f the pro- 

son to a n um^Hca? pe^/ur^aho n ^ afh f ^ e d ’ * COmpari ‘ 
robotic manipulators Also damn b ° d ' tbrou 9h four 
to control desi g r n'add Jion a 2! ^ e , d * i,S app 'ioation 
linearization is also disced ' 3 ' 9 ° rithm ,or tha 

implemented on Alliant/FX8, a shared G mo th8 a,g a ori,hm is 

sor machine, it takes 16 5 ' m ;ir memory 8-proces- 

dynamic model of a 7 dearee ' 88C0nds ,0 linearize the 
tor. degree-of-freedom robot manipula- 
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